Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,809,935 G→T 6.3% A604A (GCG→GCT katE → catalase HPII, heme d‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,809,9350GT 6.3% 708.9 / 16.1 190A604A (GCG→GCTkatEcatalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base G (83/95);  new base T (3/9);  total (86/104)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 195 total aligned reads displayed.

CCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTACC  >  CP009273/1809788‑1810082
                                                                                                                                                   |                                                                                                                                                   
ccAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCg                                                                                                                                                     <  1:1485762/148‑1 (MQ=255)
 cAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGa                                                                                                                                                    >  2:172945/1‑148 (MQ=255)
 cAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGa                                                                                                                                                    <  1:1399064/148‑1 (MQ=255)
   ggCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAtt                                                                                                                                                  >  1:303420/1‑148 (MQ=255)
   ggCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAtt                                                                                                                                                  >  1:1617369/1‑148 (MQ=255)
    gCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAttt                                                                                                                                                 >  2:1234777/1‑148 (MQ=255)
     cGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAtttt                                                                                                                                                <  2:2224380/148‑1 (MQ=255)
     cGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAtttt                                                                                                                                                <  2:3592085/148‑1 (MQ=255)
       gTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTAc                                                                                                                                              <  1:503971/148‑1 (MQ=255)
       gTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTAc                                                                                                                                              >  2:2042817/1‑148 (MQ=255)
        tGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACt                                                                                                                                             >  1:1922327/1‑148 (MQ=255)
          gCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTa                                                                                                                                           <  1:1067467/148‑1 (MQ=255)
              aaaaTCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAatga                                                                                                                                       >  1:2709549/1‑148 (MQ=255)
                aaTCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAatgatg                                                                                                                                     >  1:1011650/1‑148 (MQ=255)
                 aTCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAAtgatga                                                                                                                                    <  2:3174143/148‑1 (MQ=255)
                  tctcGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGaa                                                                                                                                   <  2:364938/148‑1 (MQ=255)
                     cGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTg                                                                                                                                <  1:1677250/148‑1 (MQ=255)
                      ggTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTga                                                                                                                               <  2:1223945/148‑1 (MQ=255)
                       gTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTgag                                                                                                                              <  1:1881348/148‑1 (MQ=255)
                        tATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTgaga                                                                                                                             <  2:1537629/148‑1 (MQ=255)
                         aTCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGAt                                                                                                                            <  2:1190379/148‑1 (MQ=255)
                         aTCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGAt                                                                                                                            >  2:2052855/1‑148 (MQ=255)
                         aTCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGAt                                                                                                                            <  1:2825899/148‑1 (MQ=255)
                          tCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATc                                                                                                                           <  1:203450/148‑1 (MQ=255)
                           cGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCg                                                                                                                          >  2:3278664/1‑148 (MQ=255)
                             aaCTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGc                                                                                                                        >  2:454783/1‑148 (MQ=255)
                                tgactgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGa                                                                                                                     >  2:1565799/1‑148 (MQ=255)
                                tgactgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGa                                                                                                                     >  2:774605/1‑148 (MQ=255)
                                 gactgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGAc                                                                                                                    >  2:1018665/1‑148 (MQ=255)
                                  actgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGAcc                                                                                                                   <  1:394696/148‑1 (MQ=255)
                                    tgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGCAGACCtt                                                                                                                 <  1:3463233/148‑1 (MQ=255)
                                      acgacCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCt                                                                                                               <  1:1988965/148‑1 (MQ=255)
                                           cAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCa                                                                                                          <  2:2663331/148‑1 (MQ=255)
                                              cTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATtc                                                                                                       >  2:2647235/1‑148 (MQ=255)
                                                 aaTATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCa                                                                                                    <  1:2768875/148‑1 (MQ=255)
                                                      cACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGgcg                                                                                               >  2:1400428/1‑148 (MQ=255)
                                                       tccccACCTGCGGACGTCAACGTTCTGTTAATGGATACTTCCTTTAGTTTGGACGCCTTTCCTGTCGGTGATGTGTAAGGTCGCGTGGTAGCGTTTTTACTTAATGATGAAGTGAGATCGGTATACCTTCTTGCAATTCTAAAGgcgc                                                                                              <  2:1618434/147‑1 (MQ=255)
                                                          ccACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGa                                                                                           <  1:862523/148‑1 (MQ=255)
                                                          ccACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGa                                                                                           >  2:2707693/1‑148 (MQ=255)
                                                          ccACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGa                                                                                           <  2:1922327/148‑1 (MQ=255)
                                                           cACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGaa                                                                                          >  1:3638997/1‑148 (MQ=255)
                                                            aCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAg                                                                                         >  2:1107780/1‑148 (MQ=255)
                                                              cTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGc                                                                                       <  2:2233028/148‑1 (MQ=255)
                                                                ccGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCa                                                                                     >  1:456273/1‑148 (MQ=255)
                                                                 cGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCaa                                                                                    <  2:2322086/148‑1 (MQ=255)
                                                                  ggACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCaaa                                                                                   <  1:1003585/148‑1 (MQ=255)
                                                                  ggACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCaaa                                                                                   <  1:1003307/148‑1 (MQ=255)
                                                                   gACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAg                                                                                  >  2:851687/1‑148 (MQ=255)
                                                                      gTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCg                                                                               <  2:1447867/148‑1 (MQ=255)
                                                                        cAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGtt                                                                             <  1:3776461/148‑1 (MQ=255)
                                                                        cAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGtt                                                                             >  2:3525589/1‑148 (MQ=255)
                                                                          aCGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCa                                                                           <  2:233418/148‑1 (MQ=255)
                                                                            ggTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATg                                                                         <  1:2356279/148‑1 (MQ=255)
                                                                                 gAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCaaa                                                                    <  2:2242067/148‑1 (MQ=255)
                                                                                 gAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCaaa                                                                    <  2:3316503/148‑1 (MQ=255)
                                                                                 gAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCaaa                                                                    <  2:2598072/148‑1 (MQ=255)
                                                                                 gAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCaaa                                                                    <  1:2663/148‑1 (MQ=255)
                                                                                   aaaaGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACt                                                                  >  2:2082463/1‑148 (MQ=255)
                                                                                    aaaGGATCCATCCTTAAGTTTGTACTACATTCCTGACGTTGATGTGAAAGGTCTCTTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCTGTCATGCCCAACTg                                                                 >  1:312902/1‑148 (MQ=255)
                                                                                    aaaGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTg                                                                 >  1:2241431/1‑148 (MQ=255)
                                                                                    aaaGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTg                                                                 <  2:1963740/148‑1 (MQ=255)
                                                                                    aaaGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTg                                                                 >  1:2327017/1‑148 (MQ=255)
                                                                                      aGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGct                                                               >  1:714791/1‑148 (MQ=255)
                                                                                      aGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGct                                                               >  2:3451738/1‑148 (MQ=255)
                                                                                      aGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGct                                                               >  1:1204703/1‑148 (MQ=255)
                                                                                       ggATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGctc                                                              <  2:114177/148‑1 (MQ=255)
                                                                                        gATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGctct                                                             >  2:2957527/1‑148 (MQ=255)
                                                                                           ccatccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACt                                                          >  1:322752/1‑148 (MQ=255)
                                                                                            catccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGGCAAACTGCTCTACTc                                                         >  2:1984198/1‑148 (MQ=255)
                                                                                               ccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCg                                                      >  2:1443377/1‑148 (MQ=255)
                                                                                                 ttAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGaa                                                    >  2:528654/1‑148 (MQ=255)
                                                                                                  tAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAAt                                                   >  1:1586051/1‑148 (MQ=255)
                                                                                                  tAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAAt                                                   >  2:2507733/1‑148 (MQ=255)
                                                                                                     gTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATggg                                                <  1:3771113/148‑1 (MQ=255)
                                                                                                     gTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATggg                                                >  2:1606425/1‑148 (MQ=255)
                                                                                                      tttGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGt                                               >  1:690214/1‑148 (MQ=255)
                                                                                                        tGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGa                                             <  2:563529/148‑1 (MQ=255)
                                                                                                        tGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGa                                             <  2:1881355/148‑1 (MQ=255)
                                                                                                             gCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGa                                        <  1:2010709/148‑1 (MQ=255)
                                                                                                                  tCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCg                                   >  1:3437390/1‑148 (MQ=255)
                                                                                                                  tCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCg                                   >  1:1018741/1‑148 (MQ=255)
                                                                                                                  tCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCg                                   >  2:280220/1‑148 (MQ=255)
                                                                                                                  tCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCg                                   >  1:2970563/1‑148 (MQ=255)
                                                                                                                   ccTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCgg                                  <  1:3212561/148‑1 (MQ=255)
                                                                                                                    cTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGa                                 <  2:2903226/148‑1 (MQ=255)
                                                                                                                      gACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATg                               >  1:1775960/1‑148 (MQ=255)
                                                                                                                         ggTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGacg                            <  2:584893/148‑1 (MQ=255)
                                                                                                                          gTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGacgg                           <  1:2390462/148‑1 (MQ=255)
                                                                                                                           tGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGacggt                          <  1:1234331/148‑1 (MQ=255)
                                                                                                                           tGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGacggt                          <  1:1234879/148‑1 (MQ=255)
                                                                                                                                tgAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGacggtacggt                     <  2:1617369/148‑1 (MQ=255)
                                                                                                                                     ggTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGc                <  2:3756627/148‑1 (MQ=255)
                                                                                                                                      gTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGcc               <  1:1221277/148‑1 (MQ=255)
                                                                                                                                         gcgTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCtat            <  1:76627/148‑1 (MQ=255)
                                                                                                                                              gTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCg       <  2:3790848/148‑1 (MQ=255)
                                                                                                                                                aGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCt     >  2:704251/1‑148 (MQ=255)
                                                                                                                                                aGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCt     >  1:1710195/1‑148 (MQ=255)
                                                                                                                                                aGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCt     >  1:3603979/1‑148 (MQ=255)
                                                                                                                                                aGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCt     >  1:3169534/1‑148 (MQ=255)
                                                                                                                                                   gATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTAcc  <  1:1017978/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTACC  >  CP009273/1809788‑1810082

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: