Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,133,311 T→A 7.7% L162Q (CTG→CAG)  flgJ → flagellar rod assembly protein and murein hydrolase; flagellum‑specific muramidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,133,3110TA 7.7% 538.3 / 26.1 180L162Q (CTG→CAG) flgJflagellar rod assembly protein and murein hydrolase; flagellum‑specific muramidase
Reads supporting (aligned to +/- strand):  ref base T (81/84);  new base A (8/6);  total (89/91)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

ATGAAATTCCCGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAG  >  CP009273/1133166‑1133457
                                                                                                                                                 |                                                                                                                                                  
 tGAAATTCCCGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGc                                                                                                                                                 <  2:417453/148‑1 (MQ=255)
  gAAATTCCCGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCa                                                                                                                                                <  1:2394910/148‑1 (MQ=255)
  gAAATTCCCGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCa                                                                                                                                                <  1:2395089/148‑1 (MQ=255)
  gAAATTCCCGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCa                                                                                                                                                >  1:1393126/1‑148 (MQ=255)
    aaTTCCCGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAg                                                                                                                                              >  1:2180216/1‑148 (MQ=255)
    aaTTCCCGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAg                                                                                                                                              <  1:1345182/148‑1 (MQ=255)
     aTTCCCGATCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCAAACTGGCAAGc                                                                                                                                             >  2:2871939/1‑148 (MQ=255)
          cGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCa                                                                                                                                        <  1:3006582/148‑1 (MQ=255)
            cTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCaaa                                                                                                                                      >  2:1150919/1‑148 (MQ=255)
            cTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCaaa                                                                                                                                      >  1:2626416/1‑148 (MQ=255)
             tCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAg                                                                                                                                     <  1:1507585/148‑1 (MQ=255)
             tCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAg                                                                                                                                     <  1:1507836/148‑1 (MQ=255)
             tCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAg                                                                                                                                     <  2:2359264/148‑1 (MQ=255)
             tCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAg                                                                                                                                     >  1:516018/1‑148 (MQ=255)
             tCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAg                                                                                                                                     <  1:1265710/148‑1 (MQ=255)
               gAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCg                                                                                                                                   <  2:1925224/148‑1 (MQ=255)
               gAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCg                                                                                                                                   <  2:2598857/148‑1 (MQ=255)
                aaaCTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCgg                                                                                                                                  >  1:2787256/1‑148 (MQ=255)
                aaaCTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCgg                                                                                                                                  <  1:643919/148‑1 (MQ=255)
                 aaCTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGgt                                                                                                                                 <  2:97426/148‑1 (MQ=255)
                 aaCTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGgt                                                                                                                                 <  2:3461214/148‑1 (MQ=255)
                 aaCTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGgt                                                                                                                                 <  2:1701808/148‑1 (MQ=255)
                    tgtgGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGtgtg                                                                                                                              <  1:1040836/148‑1 (MQ=255)
                         tgCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCcaca                                                                                                                         <  2:3287230/148‑1 (MQ=255)
                         tgCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCcaca                                                                                                                         <  2:1523429/148‑1 (MQ=255)
                          gCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCAcat                                                                                                                        >  1:197329/1‑148 (MQ=255)
                            gTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCAcatca                                                                                                                      >  1:3346006/1‑148 (MQ=255)
                             ttATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCAcatcat                                                                                                                     >  2:2503345/1‑148 (MQ=255)
                              tATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCAtt                                                                                                                    >  2:1014242/1‑148 (MQ=255)
                              tATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCAtt                                                                                                                    <  2:3044351/148‑1 (MQ=255)
                               aTCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCAttt                                                                                                                   >  2:932964/1‑148 (MQ=255)
                               aTCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCAttt                                                                                                                   <  1:1618372/148‑1 (MQ=255)
                               aTCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCAttt                                                                                                                   >  1:2746650/1‑148 (MQ=255)
                                tCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTg                                                                                                                  <  2:3075933/148‑1 (MQ=255)
                                tCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTg                                                                                                                  <  2:2180216/148‑1 (MQ=255)
                                   aaaTCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATc                                                                                                               >  1:2112077/1‑148 (MQ=255)
                                   aaaTCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATc                                                                                                               >  1:2574996/1‑148 (MQ=255)
                                   aaaTCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATc                                                                                                               >  1:2112074/1‑148 (MQ=255)
                                     aTCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCt                                                                                                             >  2:1277484/1‑148 (MQ=255)
                                          gcgcTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTc                                                                                                        <  1:2352640/148‑1 (MQ=255)
                                          gcgcTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTc                                                                                                        <  1:3378873/148‑1 (MQ=255)
                                             cTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAgg                                                                                                     >  2:2994294/1‑148 (MQ=255)
                                               ttCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAggcg                                                                                                   >  2:2404983/1‑148 (MQ=255)
                                                     gCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTg                                                                                             >  1:1840773/1‑148 (MQ=255)
                                                     gCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTg                                                                                             >  1:1252830/1‑148 (MQ=255)
                                                         gTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAAt                                                                                         >  1:1402866/1‑148 (MQ=255)
                                                          tGCAAATGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATc                                                                                        <  1:3268498/148‑1 (MQ=255)
                                                          tGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATc                                                                                        >  1:3011816/1‑148 (MQ=255)
                                                          tGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATc                                                                                        >  2:2785361/1‑148 (MQ=255)
                                                             aaaaGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTgg                                                                                     >  2:1247949/1‑148 (MQ=255)
                                                             aaaaGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGAGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTgg                                                                                     >  2:3479690/1‑148 (MQ=255)
                                                               aaGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGtt                                                                                   >  2:761573/1‑148 (MQ=255)
                                                               aaGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGtt                                                                                   >  2:1957640/1‑148 (MQ=255)
                                                               aaGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGtt                                                                                   >  2:2545387/1‑148 (MQ=255)
                                                                aGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTg                                                                                  <  1:539736/148‑1 (MQ=255)
                                                                  gCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTggg                                                                                >  1:1658047/1‑148 (MQ=255)
                                                                   ccGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTgggg                                                                               <  2:2059245/148‑1 (MQ=255)
                                                                    cGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTggggg                                                                              <  1:1237943/148‑1 (MQ=255)
                                                                       gCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCaa                                                                           <  1:1741899/148‑1 (MQ=255)
                                                                        ccACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAAc                                                                          >  1:986321/1‑148 (MQ=255)
                                                                          aCGTAACTACGATAATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACgg                                                                        >  2:1514910/1‑148 (MQ=255)
                                                                           cGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGc                                                                       >  2:1286159/1‑148 (MQ=255)
                                                                               aCTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAAt                                                                   <  1:3067532/148‑1 (MQ=255)
                                                                                cTACGATGAGTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATc                                                                  <  2:1831939/148‑1 (MQ=255)
                                                                                 tACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATcc                                                                 >  2:1875127/1‑148 (MQ=255)
                                                                                   cGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATccgc                                                               <  2:516018/148‑1 (MQ=255)
                                                                                   cGATGAGTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATccgc                                                               <  1:945291/148‑1 (MQ=255)
                                                                                     atgatTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATccgccg                                                             >  1:1001949/1‑148 (MQ=255)
                                                                                     atgatTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATccgccg                                                             <  1:907131/148‑1 (MQ=255)
                                                                                     atgatTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATccgccg                                                             >  1:1001942/1‑148 (MQ=255)
                                                                                            gCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAAc                                                      <  1:517763/148‑1 (MQ=255)
                                                                                             cTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACg                                                     <  2:673098/148‑1 (MQ=255)
                                                                                              tGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACgg                                                    <  2:197329/148‑1 (MQ=255)
                                                                                               gCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGc                                                   >  1:1556525/1‑148 (MQ=255)
                                                                                                      gACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGcgagccga                                            <  2:1728540/148‑1 (MQ=255)
                                                                                                       aCAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGcgagccgag                                           <  2:238331/148‑1 (MQ=255)
                                                                                                             aaGCATTCCTCGCGCAACTCGCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCTGGCGGCACTGGAATCTGGTTGGGGGCAACGGAAAATCCGCCGCGAAAACGGCGAGCCGAACTATaa                                     <  1:967408/148‑1 (MQ=255)
                                                                                                             aaGCAGTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATaa                                     <  1:582724/148‑1 (MQ=255)
                                                                                                              aGCATTCCTCGCGCAACTCTCGCTTCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTGGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAAc                                    >  1:2039684/1‑148 (MQ=255)
                                                                                                              aGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAAc                                    <  1:71892/148‑1 (MQ=255)
                                                                                                               gCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAAcc                                   >  1:1836021/1‑148 (MQ=255)
                                                                                                                cATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCt                                  <  2:537162/148‑1 (MQ=255)
                                                                                                                cATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCt                                  <  1:2602545/148‑1 (MQ=255)
                                                                                                                 aTTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTg                                 <  1:1014096/148‑1 (MQ=255)
                                                                                                                  ttCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGt                                <  2:1852548/148‑1 (MQ=255)
                                                                                                                   tCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGtt                               >  1:625169/1‑148 (MQ=255)
                                                                                                                         cgcAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGgtgt                         <  2:554362/148‑1 (MQ=255)
                                                                                                                               tctcGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGc                   >  2:1084775/1‑148 (MQ=255)
                                                                                                                                ctcGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGcc                  >  1:861043/1‑148 (MQ=255)
                                                                                                                                 tcGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCct                 <  2:2619158/148‑1 (MQ=255)
                                                                                                                                  cGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCctc                >  2:460913/1‑148 (MQ=255)
                                                                                                                                  cGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCctc                <  2:1403002/148‑1 (MQ=255)
                                                                                                                                   gCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCctct               >  1:670176/1‑148 (MQ=255)
                                                                                                                                   gCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCctct               >  1:670177/1‑148 (MQ=255)
                                                                                                                                        ccgcccAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCaa          <  2:2159358/148‑1 (MQ=255)
                                                                                                                                            ccAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTgg      <  1:2265126/148‑1 (MQ=255)
                                                                                                                                              aaCTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGaa    <  1:308391/148‑1 (MQ=255)
                                                                                                                                               aCTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGaaa   >  2:597827/1‑148 (MQ=255)
                                                                                                                                               aCTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGaaa   >  2:1551082/1‑148 (MQ=255)
                                                                                                                                                cTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAg  >  2:244336/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ATGAAATTCCCGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAG  >  CP009273/1133166‑1133457

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