Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 169,813 A→G 6.8% intergenic (+777/+2910) PP_RS08750 → / ← fhuA_1 tRNA‑Val/Ferrichrome outer membrane transporter/phagereceptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122169,8130AG 6.8% 253.5 / 6.9 88intergenic (+777/+2910)PP_RS08750/fhuA_1tRNA‑Val/Ferrichrome outer membrane transporter/phagereceptor
Reads supporting (aligned to +/- strand):  ref base A (43/39);  new base G (3/3);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTT  >  ACB122/169684‑169923
                                                                                                                                 |                                                                                                              
cATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAAGAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACtgtg                                                                                                           >  8:359668/1‑135 (MQ=17)
     tccAAGGCTAAATACTACTGACTGACCGATAGTGCCCCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTGGTTAgg                                                                                                      >  3:316851/1‑135 (MQ=12)
     tccAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAgg                                                                                                      >  8:266254/1‑135 (MQ=32)
     tccAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAgg                                                                                                      >  8:23288/1‑135 (MQ=32)
       cAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTg                                                                                                    >  8:18957/1‑135 (MQ=32)
           gCTAAATACTACTGCCTGACCGATAGTGAACCAGTACCGTGAGGGAAAGCCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTg                                                                                                <  7:338740/135‑1 (MQ=17)
                 tactacTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTAcc                                                                                          <  1:4917/135‑1 (MQ=32)
                  actactGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCt                                                                                         >  1:292960/1‑135 (MQ=32)
                  actactGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCt                                                                                         >  8:199182/1‑135 (MQ=32)
                       tgactgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGt                                                                                    <  7:199182/135‑1 (MQ=32)
                         actgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGta                                                                                   >  7:156261/1‑134 (MQ=32)
                         actgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGta                                                                                   >  5:255837/1‑134 (MQ=32)
                         actgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGta                                                                                   >  5:147961/1‑134 (MQ=32)
                          ctgacCGAAACTCAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGta                                                                                   <  4:76620/133‑1 (MQ=11)
                            gacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAAt                                                                               >  2:13742/1‑135 (MQ=32)
                            gacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAAt                                                                               >  2:201369/1‑135 (MQ=32)
                            gacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAAt                                                                               >  8:291079/1‑135 (MQ=32)
                                gATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATggg                                                                            >  7:76810/1‑134 (MQ=32)
                                 aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACt                                                                                                 <  7:192221/112‑1 (MQ=32)
                                 aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACt                                                                                                 >  8:192221/1‑112 (MQ=32)
                                 aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          <  1:201369/135‑1 (MQ=32)
                                 aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          <  6:370941/135‑1 (MQ=32)
                                 aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          <  7:86491/135‑1 (MQ=32)
                                 aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          <  6:285118/135‑1 (MQ=32)
                                 aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          <  8:156261/135‑1 (MQ=32)
                                 aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          <  2:328694/135‑1 (MQ=32)
                                      gAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGa                                                                     <  3:228797/135‑1 (MQ=32)
                                          cAGTACCGTGAGGGAAAGGCGATAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACtt                                                                  >  2:172525/1‑134 (MQ=11)
                                          cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCTAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          >  1:316810/1‑126 (MQ=11)
                                          cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCTAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          <  2:316810/126‑1 (MQ=11)
                                          cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACtt                                                                  >  3:366990/1‑134 (MQ=32)
                                               ccGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTatat                                                              >  6:147265/1‑133 (MQ=32)
                                               ccGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTatat                                                              >  1:269345/1‑133 (MQ=32)
                                               ccGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTatat                                                              >  3:328482/1‑133 (MQ=32)
                                                cGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCa                                                           <  5:374725/135‑1 (MQ=32)
                                                  tGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTCAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                         <  1:285992/135‑1 (MQ=17)
                                                     ggtaaaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          <  7:139938/112‑1 (MQ=17)
                                                     gggAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGtag                                                       >  7:306118/1‑132 (MQ=11)
                                                     gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          >  8:139938/1‑115 (MQ=32)
                                                     gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTgg                                                       >  7:346877/1‑134 (MQ=32)
                                                     gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTgg                                                       >  3:130965/1‑134 (MQ=32)
                                                     gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTgg                                                       >  1:6653/1‑134 (MQ=32)
                                                     gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                      <  4:14273/135‑1 (MQ=32)
                                                     gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                      >  1:362182/1‑135 (MQ=32)
                                                       gAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTCCCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGa                                                    >  8:137933/1‑135 (MQ=17)
                                                       gAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGa                                                    <  7:289684/135‑1 (MQ=32)
                                                        aaaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAg                                                   <  8:76810/135‑1 (MQ=32)
                                                          gggCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGc                                                  >  7:246742/2‑134 (MQ=17)
                                                          acgCGAACAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCt                                                 <  4:230959/133‑1 (MQ=17)
                                                          aGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACt                                                                                                 >  4:252221/1‑87 (MQ=33)
                                                          aGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACt                                                                                                 <  3:256968/87‑1 (MQ=33)
                                                          aGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACt                                                                                                 >  4:256968/1‑87 (MQ=33)
                                                          aGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACt                                                                                                 <  3:252221/87‑1 (MQ=33)
                                                          aGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGa                                                                     >  1:39564/1‑115 (MQ=32)
                                                          aGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGa                                                                     <  2:39564/115‑1 (MQ=32)
                                                          aGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCTAGCt                                                 <  1:55523/135‑1 (MQ=17)
                                                          aGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGc                                                  >  8:167329/1‑134 (MQ=32)
                                                            gCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTa                                               >  6:109294/1‑135 (MQ=32)
                                                             cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTaa                                              >  2:295540/1‑135 (MQ=32)
                                                              gAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAAc                                             >  8:111388/1‑135 (MQ=32)
                                                               aaaaGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAAcc                                            <  7:337552/135‑1 (MQ=32)
                                                                         cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGGACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGa                                  >  4:67305/1‑135 (MQ=17)
                                                                               ggggAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCCGTGGCGAGCTTAACCGAATAGGGGAGGCGTa                            <  5:262889/135‑1 (MQ=17)
                                                                                   aGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGTGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCg                         >  8:140237/1‑134 (MQ=17)
                                                                                   aGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTtat                                                                <  6:330937/95‑1 (MQ=33)
                                                                                   aGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTtat                                                                >  5:330937/1‑95 (MQ=33)
                                                                                   aGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGa                        <  2:94046/135‑1 (MQ=32)
                                                                                   aGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGa                        <  7:326950/135‑1 (MQ=32)
                                                                                   aGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGa                        <  4:130965/135‑1 (MQ=32)
                                                                                   aGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGa                        <  3:111410/135‑1 (MQ=32)
                                                                                   aGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGa                        <  1:172525/135‑1 (MQ=11)
                                                                                    gTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCg                         >  4:254435/1‑133 (MQ=32)
                                                                                     tGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGaaa                      <  8:236250/135‑1 (MQ=32)
                                                                                      gAAATAGAACCTGACACCGTGTGCGTTCAAGGAGTGGGAGCCTACTTTGTTATGTGACTGCGTACCTTTTGTATAATGGGTCAGCGTCTTATCTTCAGTGGCGAGGTTAACCGAATAGGGGAGGCGTAGCGaaa                      <  8:183130/134‑1 (MQ=11)
                                                                                      gAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAg                     >  7:274066/1‑135 (MQ=32)
                                                                                      gAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAg                     >  6:118037/1‑135 (MQ=32)
                                                                                       aaaTAGAACCTGAACCCGAGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGc                    <  2:103825/135‑1 (MQ=17)
                                                                                        aaTAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGc                    <  7:291079/134‑1 (MQ=32)
                                                                                             aaCCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTc               >  6:277583/1‑134 (MQ=32)
                                                                                                cTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTg                                                                                                <  5:97434/50‑1 (MQ=38)
                                                                                                cTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTg                                                                                                >  6:97434/1‑50 (MQ=38)
                                                                                                cTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCtt                                                                                        >  8:117284/1‑58 (MQ=35)
                                                                                                cTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCtt                                                                                        <  7:117284/58‑1 (MQ=35)
                                                                                                     aCCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTAATAggg      <  6:194083/135‑1 (MQ=32)
                                                                                                     aCCGTGTGCGTACAAGCAGGGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTGAATAggg      >  4:234365/1‑135 (MQ=17)
                                                                                                         tgtgCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGtt  <  8:346877/135‑1 (MQ=31)
                                                                                                         tgtgCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAAAGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGtt  <  6:147961/135‑1 (MQ=14)
                                                                                                         tgtgCGTACAAGCAGTGGGAGACTACTTTGTTCGGTGTCTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTAATAGCGCGtt  <  6:117353/135‑1 (MQ=255)
                                                                                                                                 |                                                                                                              
CATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTT  >  ACB122/169684‑169923

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: