Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,077,212 A→T 100% N326Y (AAT→TAT)  sbcB → exodeoxyribonuclease I; exonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,077,2120AT95.6% 714.2 / 17.3 225N326Y (AAT→TAT) sbcBexodeoxyribonuclease I; exonuclease I
Reads supporting (aligned to +/- strand):  ref base A (7/3);  new base T (109/105);  total (116/109)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 225 total aligned reads displayed.

CACATTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTG  >  CP009273/2077067‑2077356
                                                                                                                                                 |                                                                                                                                                
cacaTTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTAt                                                                                                                                                >  2:1059229/1‑148 (MQ=255)
cacaTTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTAt                                                                                                                                                <  1:2085035/148‑1 (MQ=255)
   aTTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGt                                                                                                                                             <  1:883915/148‑1 (MQ=255)
    ttGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTc                                                                                                                                            <  2:3142726/148‑1 (MQ=255)
     tGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCa                                                                                                                                           >  2:3119016/1‑148 (MQ=255)
           aGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTg                                                                                                                                     <  1:3480960/148‑1 (MQ=255)
           aGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTg                                                                                                                                     <  1:1617645/148‑1 (MQ=255)
            gCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGc                                                                                                                                    <  1:3300971/148‑1 (MQ=255)
            gCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGc                                                                                                                                    >  2:2178208/1‑148 (MQ=255)
            gCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGc                                                                                                                                    <  2:1804758/148‑1 (MQ=255)
              gTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCt                                                                                                                                  <  1:778096/148‑1 (MQ=255)
               tttATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGACTTTATCGTCAGCATTGCCTc                                                                                                                                 >  2:815666/1‑148 (MQ=255)
                ttATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCg                                                                                                                                >  1:1569271/1‑148 (MQ=255)
                ttATAGACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCg                                                                                                                                <  2:2024967/148‑1 (MQ=255)
                    ataCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATaa                                                                                                                            >  2:1905912/1‑148 (MQ=255)
                    ataCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATaa                                                                                                                            >  2:1976108/1‑148 (MQ=255)
                    ataCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATaa                                                                                                                            >  2:45428/1‑148 (MQ=255)
                      aCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAAcc                                                                                                                          >  2:1876576/1‑148 (MQ=255)
                      aCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAAcc                                                                                                                          >  2:1896668/1‑148 (MQ=255)
                            aaaaCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGaaaa                                                                                                                    <  1:567619/148‑1 (MQ=255)
                             aaaCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAAt                                                                                                                   >  1:1720723/1‑148 (MQ=255)
                             aaaCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAAAAAATGTCCGGTGCTGGCCCAGGCGAATACCCAAAGCCCGGAACATGCCGACAGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAAt                                                                                                                   >  1:1004143/1‑148 (MQ=255)
                               aCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTc                                                                                                                 >  1:3465702/1‑148 (MQ=255)
                                ccGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCt                                                                                                                >  1:3757148/1‑148 (MQ=255)
                                ccGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCt                                                                                                                >  1:1003967/1‑148 (MQ=255)
                                 cGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTg                                                                                                               <  2:1599784/148‑1 (MQ=255)
                                   aTCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCg                                                                                                             >  2:971486/1‑148 (MQ=255)
                                    tCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGt                                                                                                            >  1:3067490/1‑148 (MQ=255)
                                    tCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGt                                                                                                            >  1:1865695/1‑148 (MQ=255)
                                    tCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGt                                                                                                            >  1:2885440/1‑148 (MQ=255)
                                       tGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGaa                                                                                                         >  1:595605/1‑148 (MQ=255)
                                             tAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCg                                                                                                   <  1:1316059/148‑1 (MQ=255)
                                             tAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCg                                                                                                   >  2:1301487/1‑148 (MQ=255)
                                             tAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCg                                                                                                   >  1:2716616/1‑148 (MQ=255)
                                             tAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCg                                                                                                   >  1:3591977/1‑148 (MQ=255)
                                              aaCGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGc                                                                                                  >  2:2886327/1‑148 (MQ=255)
                                              aaCGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGc                                                                                                  <  2:2222116/148‑1 (MQ=255)
                                                cgccgcCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCaa                                                                                                >  2:904942/1‑148 (MQ=255)
                                                 gccgccGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAg                                                                                               <  2:2620085/148‑1 (MQ=255)
                                                  ccgccgTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGt                                                                                              <  2:2286492/148‑1 (MQ=255)
                                                   cgccgTGCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTg                                                                                             <  2:3723429/148‑1 (MQ=255)
                                                        ttCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGa                                                                                        >  2:2986144/1‑148 (MQ=255)
                                                         tCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGaa                                                                                       >  1:3491139/1‑148 (MQ=255)
                                                         tCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGACCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGaa                                                                                       >  1:2405916/1‑148 (MQ=255)
                                                          ccGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGaaa                                                                                      >  2:539968/1‑148 (MQ=255)
                                                          ccGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGaaa                                                                                      >  2:3377972/1‑148 (MQ=255)
                                                              ttAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAgt                                                                                  <  1:1571583/148‑1 (MQ=255)
                                                                 aGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAgtggt                                                                               >  2:604316/1‑148 (MQ=255)
                                                                  gCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAgtggtg                                                                              >  2:352459/1‑148 (MQ=255)
                                                                  gCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAgtggtg                                                                              >  1:3310087/1‑148 (MQ=255)
                                                                  gCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAgtggtg                                                                              >  2:290896/1‑148 (MQ=255)
                                                                     ggTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCg                                                                           >  2:2636327/1‑148 (MQ=255)
                                                                     ggTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCg                                                                           <  1:2038847/148‑1 (MQ=255)
                                                                      gTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGa                                                                          <  2:1929703/148‑1 (MQ=255)
                                                                        gCAGATCAATACATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGata                                                                        <  2:2078177/148‑1 (MQ=255)
                                                                         cATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGatat                                                                       <  1:506852/148‑1 (MQ=255)
                                                                          atatCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATAtt                                                                      <  1:1271196/148‑1 (MQ=255)
                                                                             tCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTcgc                                                                   <  2:2538931/148‑1 (MQ=255)
                                                                              cAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTcgcg                                                                  >  2:699218/1‑148 (MQ=255)
                                                                               aaTAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCgg                                                                 >  2:1938237/1‑148 (MQ=255)
                                                                                aTAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGa                                                                >  2:1055732/1‑148 (MQ=255)
                                                                                aTAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGa                                                                >  1:2556152/1‑148 (MQ=255)
                                                                                  aaaTGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAg                                                              <  1:229113/148‑1 (MQ=255)
                                                                                  aaaTGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAg                                                              >  2:3490856/1‑148 (MQ=255)
                                                                                   aaTGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGc                                                             <  2:1146980/148‑1 (MQ=255)
                                                                                    aTGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATGGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGcc                                                            >  1:2830266/1‑148 (MQ=255)
                                                                                     tGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCg                                                           >  1:3100543/1‑148 (MQ=255)
                                                                                      gTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGa                                                          <  2:2495968/148‑1 (MQ=255)
                                                                                       tCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGaa                                                         <  2:1003967/148‑1 (MQ=255)
                                                                                        ccGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAAc                                                        <  2:1530266/148‑1 (MQ=255)
                                                                                          ggTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCg                                                      <  2:2202670/148‑1 (MQ=255)
                                                                                          ggTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCg                                                      <  2:3646723/148‑1 (MQ=255)
                                                                                           gTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGt                                                     <  2:463706/148‑1 (MQ=255)
                                                                                            tGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGtt                                                    >  2:1575369/1‑148 (MQ=255)
                                                                                             gCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGttt                                                   <  1:1905912/148‑1 (MQ=255)
                                                                                                 gCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGc                                               >  1:2731928/1‑148 (MQ=255)
                                                                                                  cccAGGCGAATACGCTACGCCCGGAAGATGTCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGcc                                              <  1:3200048/148‑1 (MQ=255)
                                                                                                      ggCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCa                                          >  2:885201/1‑148 (MQ=255)
                                                                                                       gCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAg                                         <  1:304371/148‑1 (MQ=255)
                                                                                                          aaTACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATa                                      >  2:3436228/1‑148 (MQ=255)
                                                                                                                 tacgcCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGa                               <  1:45428/148‑1 (MQ=255)
                                                                                                                 tacgcCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGa                               <  1:1427051/148‑1 (MQ=255)
                                                                                                                 tacgcCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGa                               <  1:1976108/148‑1 (MQ=255)
                                                                                                                   cgcCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATg                             >  2:1278148/1‑148 (MQ=255)
                                                                                                                      ccGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGcac                          <  1:1715106/148‑1 (MQ=255)
                                                                                                                         gAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGc                       <  2:667098/148‑1 (MQ=255)
                                                                                                                            gATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCttt                    <  1:971486/148‑1 (MQ=255)
                                                                                                                             aTGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTta                   >  1:1912199/1‑148 (MQ=255)
                                                                                                                               gCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTtata                 <  2:3591977/148‑1 (MQ=255)
                                                                                                                                ccgaccgaCTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATaa                <  1:2009433/148‑1 (MQ=255)
                                                                                                                                ccgaccgaCTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATaa                <  2:2280247/148‑1 (MQ=255)
                                                                                                                                   accgacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACGGCTTTATAACgg             >  2:2793192/1‑148 (MQ=255)
                                                                                                                                    ccgacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGc            <  2:1367280/148‑1 (MQ=255)
                                                                                                                                    ccgacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGc            <  2:3285581/148‑1 (MQ=255)
                                                                                                                                    ccgacTGGGAATTGATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGc            <  2:1115461/148‑1 (MQ=255)
                                                                                                                                     cgacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCt           <  1:3409858/148‑1 (MQ=255)
                                                                                                                                      ggcTGGGAATTAATCGTCAGCAGTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCtt          <  2:1731741/146‑1 (MQ=255)
                                                                                                                                       acTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCttt         <  2:3757148/148‑1 (MQ=255)
                                                                                                                                            gAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAg    <  1:1107359/148‑1 (MQ=255)
                                                                                                                                            gAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAg    >  2:2665524/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CACATTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTG  >  CP009273/2077067‑2077356

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: