Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,077,212 A→T 88.4% 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,2120AT88.4% 593.3 / 58.7 206N326Y (AAT→TAT) sbcBexodeoxyribonuclease I; exonuclease I
Reads supporting (aligned to +/- strand):  ref base A (5/19);  new base T (84/98);  total (89/117)
Fisher's exact test for biased strand distribution p-value = 2.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 206 total aligned reads displayed.

CACATTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTG  >  CP009273/2077067‑2077356
                                                                                                                                                 |                                                                                                                                                
   aTTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGt                                                                                                                                             >  1:3178532/1‑148 (MQ=255)
    ttGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTc                                                                                                                                            >  1:2403657/1‑148 (MQ=255)
    ttGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTc                                                                                                                                            >  1:2405705/1‑148 (MQ=255)
        gcgAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCa                                                                                                                                        <  2:2298091/148‑1 (MQ=255)
           aGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTg                                                                                                                                     <  1:1790968/148‑1 (MQ=255)
             cGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGcc                                                                                                                                   <  1:1600318/148‑1 (MQ=255)
               tttATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTc                                                                                                                                 >  1:250196/1‑148 (MQ=255)
               tttATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTc                                                                                                                                 >  1:2587018/1‑148 (MQ=255)
               tttATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTc                                                                                                                                 >  1:1918875/1‑148 (MQ=255)
                ttATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCg                                                                                                                                >  2:3296238/1‑148 (MQ=255)
                ttATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGATGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCg                                                                                                                                <  1:7770/148‑1 (MQ=255)
                   tataCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATa                                                                                                                             <  1:1098046/148‑1 (MQ=255)
                   tataCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATa                                                                                                                             <  1:639701/148‑1 (MQ=255)
                   tataCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATa                                                                                                                             <  1:1257539/148‑1 (MQ=255)
                    ataCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATaa                                                                                                                            >  1:819645/1‑148 (MQ=255)
                    ataCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATaa                                                                                                                            >  2:3138685/1‑148 (MQ=255)
                    ataCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATaa                                                                                                                            >  2:419975/1‑148 (MQ=255)
                     taCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAAc                                                                                                                           <  1:2939809/148‑1 (MQ=255)
                           aaaaaCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGaaa                                                                                                                     <  2:1228668/148‑1 (MQ=255)
                           aaaaaCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGaaa                                                                                                                     <  1:1294974/148‑1 (MQ=255)
                            aaaaCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGaaaa                                                                                                                    <  2:188080/148‑1 (MQ=255)
                                      ttGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGa                                                                                                          <  2:1837552/148‑1 (MQ=255)
                                      ttGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGa                                                                                                          <  2:2166043/148‑1 (MQ=255)
                                      ttGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGa                                                                                                          <  2:891619/148‑1 (MQ=255)
                                       tGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGaa                                                                                                         >  2:2031578/1‑148 (MQ=255)
                                       tGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGaa                                                                                                         >  1:3069803/1‑148 (MQ=255)
                                        ggCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGaaa                                                                                                        <  2:3305768/148‑1 (MQ=255)
                                         gCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGaaaa                                                                                                       >  2:1089855/1‑148 (MQ=255)
                                         gCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGaaaa                                                                                                       <  1:278406/148‑1 (MQ=255)
                                         gCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGaaaa                                                                                                       <  1:462796/148‑1 (MQ=255)
                                          cGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAAt                                                                                                      <  1:3109626/148‑1 (MQ=255)
                                           gATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGTCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATc                                                                                                     >  1:725439/1‑148 (MQ=255)
                                            aTAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATcc                                                                                                    >  1:1854898/1‑148 (MQ=255)
                                            aTAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATcc                                                                                                    >  2:2075142/1‑148 (MQ=255)
                                             tAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCg                                                                                                   >  1:1368687/1‑148 (MQ=255)
                                               aCGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCa                                                                                                 <  1:55903/148‑1 (MQ=255)
                                                     ccgTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTgcg                                                                                           <  2:3183435/148‑1 (MQ=255)
                                                      cgTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTgcgc                                                                                          >  2:2301301/1‑148 (MQ=255)
                                                         tCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGaa                                                                                       >  2:1030070/1‑148 (MQ=255)
                                                           cGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGaaaa                                                                                     <  2:1923507/148‑1 (MQ=255)
                                                           cGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGaaaa                                                                                     <  1:453313/148‑1 (MQ=255)
                                                           cGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGaaaa                                                                                     <  1:279506/148‑1 (MQ=255)
                                                           cGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGaaaa                                                                                     <  1:2301301/148‑1 (MQ=255)
                                                             gTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAg                                                                                   >  2:2153230/1‑148 (MQ=255)
                                                             gTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAg                                                                                   <  1:1538141/148‑1 (MQ=255)
                                                             gTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAg                                                                                   <  2:2920248/148‑1 (MQ=255)
                                                              ttAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAgt                                                                                  >  1:57136/1‑148 (MQ=255)
                                                                aaGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAgtgg                                                                                >  1:2308442/1‑148 (MQ=255)
                                                                   cTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGtggtgg                                                                             <  2:1199579/148‑1 (MQ=255)
                                                                     ggTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCg                                                                           >  2:225555/1‑148 (MQ=255)
                                                                        gCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGata                                                                        >  1:414058/1‑148 (MQ=255)
                                                                        gCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGata                                                                        <  2:1271321/148‑1 (MQ=255)
                                                                        gCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGata                                                                        <  2:3226215/148‑1 (MQ=255)
                                                                         cATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGatat                                                                       >  2:3174428/1‑148 (MQ=255)
                                                                          atatCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAACAGTGGTGGCGATAtt                                                                      >  2:2352370/1‑148 (MQ=255)
                                                                           tatCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTc                                                                     >  2:2594350/1‑148 (MQ=255)
                                                                             tCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTcgc                                                                   <  2:1954567/148‑1 (MQ=255)
                                                                             tCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTcgc                                                                   >  1:976334/1‑148 (MQ=255)
                                                                              cAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTcgcg                                                                  >  1:2596011/1‑148 (MQ=255)
                                                                               aaTAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCgg                                                                 <  2:1475760/148‑1 (MQ=255)
                                                                               aaTAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCgg                                                                 >  2:353923/1‑148 (MQ=255)
                                                                               aaTAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCgg                                                                 >  1:1451352/1‑148 (MQ=255)
                                                                                aTAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGa                                                                >  1:1911061/1‑148 (MQ=255)
                                                                                aTAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGa                                                                <  2:3306833/148‑1 (MQ=255)
                                                                                 tAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGaa                                                               <  1:2293376/148‑1 (MQ=255)
                                                                                    aTGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGcc                                                            <  1:2841074/148‑1 (MQ=255)
                                                                                    aTGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGcc                                                            <  1:2841696/148‑1 (MQ=255)
                                                                                      gTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGa                                                          <  1:953849/148‑1 (MQ=255)
                                                                                             gCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGttt                                                   <  2:2214998/148‑1 (MQ=255)
                                                                                             gCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGttt                                                   <  2:284142/148‑1 (MQ=255)
                                                                                                 gCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGc                                               >  1:1708051/1‑148 (MQ=255)
                                                                                                  cccAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGcc                                              <  2:1206787/148‑1 (MQ=255)
                                                                                                       gCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAg                                         >  1:1191046/1‑148 (MQ=255)
                                                                                                       gCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAg                                         >  2:2500632/1‑148 (MQ=255)
                                                                                                          aaTACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATa                                      <  2:625265/148‑1 (MQ=255)
                                                                                                          aaTACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATa                                      <  2:2467213/148‑1 (MQ=255)
                                                                                                          aaTACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATa                                      >  2:2326648/1‑148 (MQ=255)
                                                                                                             acgctacgcCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGACCCGTTTACGCCTTCAGATAACg                                   >  1:336196/1‑148 (MQ=255)
                                                                                                              cgctacgcCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGt                                  <  1:2031578/148‑1 (MQ=255)
                                                                                                                  acgcCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGAt                              <  1:1849574/148‑1 (MQ=255)
                                                                                                                    gcCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGc                            <  2:2586906/148‑1 (MQ=255)
                                                                                                                       cGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGcaca                         <  1:2230365/148‑1 (MQ=255)
                                                                                                                           aGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCtt                     <  1:2831804/148‑1 (MQ=255)
                                                                                                                            gATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCttt                    >  1:904007/1‑148 (MQ=255)
                                                                                                                               gCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTtata                 >  1:1510869/1‑148 (MQ=255)
                                                                                                                               gCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTtata                 <  1:1737514/148‑1 (MQ=255)
                                                                                                                                ccgaccgaCTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATaa                <  2:1605958/148‑1 (MQ=255)
                                                                                                                                 cgaccgacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAAc               <  2:1850885/148‑1 (MQ=255)
                                                                                                                                 cgaccgacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAAc               <  1:211370/148‑1 (MQ=255)
                                                                                                                                    ccgacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGc            >  2:1100277/1‑148 (MQ=255)
                                                                                                                                    ccgacTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGc            <  1:1765547/148‑1 (MQ=255)
                                                                                                                                     cgacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCt           <  2:27718/148‑1 (MQ=255)
                                                                                                                                      gacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCtt          <  2:819333/148‑1 (MQ=255)
                                                                                                                                      gacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCtt          <  2:819656/148‑1 (MQ=255)
                                                                                                                                       acTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCttt         <  2:3150098/148‑1 (MQ=255)
                                                                                                                                       acTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCttt         <  2:2707803/148‑1 (MQ=255)
                                                                                                                                       acTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCttt         <  2:1597432/148‑1 (MQ=255)
                                                                                                                                             aaTTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGt   <  2:2521300/148‑1 (MQ=255)
                                                                                                                                              aTTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTg  <  2:57136/148‑1 (MQ=255)
                                                                                                                                              aTTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTg  >  1:1062157/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: