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,2120AT97.2% 1067.9 / 12.7 330N326Y (AAT→TAT) sbcBexodeoxyribonuclease I; exonuclease I
Reads supporting (aligned to +/- strand):  ref base A (6/4);  new base T (145/175);  total (151/179)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 330 total aligned reads displayed.

GACACATTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGAC  >  CP009273/2077065‑2077358
                                                                                                                                                   |                                                                                                                                                  
   acaTTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGt                                                                                                                                               <  1:1072807/150‑1 (MQ=255)
   acaTTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGt                                                                                                                                               <  2:3200500/150‑1 (MQ=255)
       tGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGc                                                                                                                                           >  2:3827422/1‑150 (MQ=255)
          gcgAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCAtt                                                                                                                                        <  2:1040497/150‑1 (MQ=255)
             aGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGcc                                                                                                                                     <  2:1378798/150‑1 (MQ=255)
             aGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGcc                                                                                                                                     <  1:3257579/150‑1 (MQ=255)
               cGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTc                                                                                                                                   >  1:1393122/1‑150 (MQ=255)
                gTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCg                                                                                                                                  <  1:3857984/150‑1 (MQ=255)
                 tttATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGa                                                                                                                                 >  2:4450478/1‑150 (MQ=255)
                 ttgatataCCGCAAAACCCGATCTGGGCGATAACGCCGCCGTGCCGGTTAAGCTGGTGCATATCAAAAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGa                                                                                                                                 <  2:1953220/147‑1 (MQ=255)
                     tataCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAAc                                                                                                                             >  2:550591/1‑150 (MQ=255)
                     tataCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAAc                                                                                                                             >  2:4585249/1‑150 (MQ=255)
                      ataCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAAcc                                                                                                                            <  1:3003313/150‑1 (MQ=255)
                        aCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTg                                                                                                                          <  2:3883071/150‑1 (MQ=255)
                        aCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTg                                                                                                                          >  2:3392486/1‑150 (MQ=255)
                         ccGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGa                                                                                                                         >  2:615634/1‑150 (MQ=255)
                          cGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGaa                                                                                                                        >  2:4252731/1‑150 (MQ=255)
                          cGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGaa                                                                                                                        <  2:1452724/150‑1 (MQ=255)
                          cGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGaa                                                                                                                        <  2:267086/150‑1 (MQ=255)
                             aaaaaCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAAt                                                                                                                     <  1:648306/150‑1 (MQ=255)
                             aaaaaCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAAt                                                                                                                     >  1:40623/1‑150 (MQ=255)
                                 aCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTg                                                                                                                 >  2:1568857/1‑150 (MQ=255)
                                 aCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTg                                                                                                                 >  2:1096570/1‑150 (MQ=255)
                                  ccGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGc                                                                                                                <  2:4677244/150‑1 (MQ=255)
                                   cGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCg                                                                                                               >  2:2752528/1‑150 (MQ=255)
                                   cGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATAAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCg                                                                                                               >  2:325240/1‑150 (MQ=255)
                                    gctCTTGGCGAGAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCGACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGGCAGCATTGCCTCGATAACCTGAAAATTCTGCGt                                                                                                              <  2:2108803/148‑1 (MQ=255)
                                    gATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGt                                                                                                              >  1:2378045/1‑150 (MQ=255)
                                       cTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGaa                                                                                                           >  1:199440/1‑150 (MQ=255)
                                          ggCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAAt                                                                                                        <  1:4152621/150‑1 (MQ=255)
                                            cGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGACAATTCTGCGTGAAAATcc                                                                                                      >  2:4698155/1‑150 (MQ=255)
                                             gATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGACTTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCg                                                                                                     >  2:450745/1‑150 (MQ=255)
                                                 aCGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAg                                                                                                 <  1:4385975/150‑1 (MQ=255)
                                                  cgccgcCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGt                                                                                                <  2:5103961/150‑1 (MQ=255)
                                                   gccgccGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTg                                                                                               >  1:2868675/1‑150 (MQ=255)
                                                   gccgccGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTg                                                                                               >  2:457540/1‑150 (MQ=255)
                                                     cgccgTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTgcg                                                                                             >  2:1933597/1‑150 (MQ=255)
                                                     cgccgTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTgcg                                                                                             <  2:2889062/150‑1 (MQ=255)
                                                     cgccgTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTgcg                                                                                             >  1:2579554/1‑150 (MQ=255)
                                                       ccgTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGcgcg                                                                                           >  2:403950/1‑150 (MQ=255)
                                                        cgTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGa                                                                                          >  1:4604758/1‑150 (MQ=255)
                                                         gTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGaa                                                                                         >  1:617901/1‑150 (MQ=255)
                                                         gTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAAACTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGaa                                                                                         >  1:4834432/1‑150 (MQ=255)
                                                             cGGTTAAGCTGGTGCATATCAATAAATGTCCGTTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAg                                                                                     <  1:2249311/150‑1 (MQ=255)
                                                                  aaGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAgtggtg                                                                                <  1:4532807/150‑1 (MQ=255)
                                                                     cTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCg                                                                             >  2:4907106/1‑150 (MQ=255)
                                                                       ggTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGat                                                                           <  1:5162643/150‑1 (MQ=255)
                                                                        gTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGata                                                                          <  2:171490/150‑1 (MQ=255)
                                                                               tCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCgg                                                                   <  1:1930751/150‑1 (MQ=255)
                                                                                cAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGa                                                                  <  2:1545832/150‑1 (MQ=255)
                                                                                 taTAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGaa                                                                 >  2:1730051/2‑150 (MQ=255)
                                                                                 aaTAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGaa                                                                 >  1:1851108/1‑150 (MQ=255)
                                                                                 aaTAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGaa                                                                 >  2:714210/1‑150 (MQ=255)
                                                                                  aTAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAg                                                                >  2:4989553/1‑150 (MQ=255)
                                                                                    aaaTGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGcc                                                              <  1:1623151/150‑1 (MQ=255)
                                                                                    aaaTGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGcc                                                              >  2:2574975/1‑150 (MQ=255)
                                                                                     aaTGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCg                                                             <  2:3297799/150‑1 (MQ=255)
                                                                                      aTGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGa                                                            >  1:3536263/1‑150 (MQ=255)
                                                                                      aTGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGa                                                            <  1:3659897/150‑1 (MQ=255)
                                                                                       tGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGaa                                                           >  2:1601561/1‑150 (MQ=255)
                                                                                              tGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTa                                                    <  1:1333963/150‑1 (MQ=255)
                                                                                                cTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACg                                                  <  1:4932373/150‑1 (MQ=255)
                                                                                                 tGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGc                                                 >  1:3042547/1‑150 (MQ=255)
                                                                                                  ggCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGcc                                                <  2:2065094/150‑1 (MQ=255)
                                                                                                     ccAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTc                                             <  1:3759032/150‑1 (MQ=255)
                                                                                                         gCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGAt                                         <  2:4308618/150‑1 (MQ=255)
                                                                                                          cGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATa                                        <  1:1642096/150‑1 (MQ=255)
                                                                                                          cGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATa                                        <  1:3992964/150‑1 (MQ=255)
                                                                                                          cGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATa                                        <  1:3951354/150‑1 (MQ=255)
                                                                                                            aaTACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAAc                                      <  1:100522/150‑1 (MQ=255)
                                                                                                            aaTACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAAc                                      <  2:3042300/150‑1 (MQ=255)
                                                                                                            aaTACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAAc                                      <  2:3152553/150‑1 (MQ=255)
                                                                                                             aTACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACg                                     <  2:3959860/150‑1 (MQ=255)
                                                                                                                cgctacgcCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTgg                                  >  1:2930861/1‑150 (MQ=255)
                                                                                                                  ctacgcCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGAt                                >  2:4379253/1‑150 (MQ=255)
                                                                                                                   tacgcCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGTAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATg                               <  2:2264381/150‑1 (MQ=255)
                                                                                                                   tacgcCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGGTTACGCCTTCAGATAACGTGGATg                               >  2:2927718/1‑150 (MQ=255)
                                                                                                                        ccGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAg                          <  1:1159127/150‑1 (MQ=255)
                                                                                                                         cGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGc                         >  1:24878/1‑150 (MQ=255)
                                                                                                                         cGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGc                         >  1:941743/1‑150 (MQ=255)
                                                                                                                            aaGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCttt                      <  1:3147295/150‑1 (MQ=255)
                                                                                                                               aTGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTtata                   <  2:4110199/150‑1 (MQ=255)
                                                                                                                                tGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATaa                  >  2:2489312/1‑150 (MQ=255)
                                                                                                                                 gCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAAc                 <  2:4015986/150‑1 (MQ=255)
                                                                                                                                    gaccgacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGc              >  2:3638993/1‑150 (MQ=255)
                                                                                                                                     accgacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCt             >  1:209084/1‑150 (MQ=255)
                                                                                                                                     accgacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCt             <  2:4324329/150‑1 (MQ=255)
                                                                                                                                        gacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCtttt          <  1:703103/150‑1 (MQ=255)
                                                                                                                                        gacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCtttt          <  1:2760232/150‑1 (MQ=255)
                                                                                                                                         acTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCttttt         <  1:1577257/150‑1 (MQ=255)
                                                                                                                                          cTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTc        >  2:2393955/1‑150 (MQ=255)
                                                                                                                                           tGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCa       <  2:1652267/150‑1 (MQ=255)
                                                                                                                                           tGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCa       >  1:4706274/1‑150 (MQ=255)
                                                                                                                                            gggAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAg      <  1:3479996/150‑1 (MQ=255)
                                                                                                                                            gggAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAg      <  2:2306321/150‑1 (MQ=255)
                                                                                                                                            gggAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAg      <  1:1096533/150‑1 (MQ=255)
                                                                                                                                            gggAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAg      <  2:1252108/150‑1 (MQ=255)
                                                                                                                                             ggAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGt     >  1:3343558/1‑150 (MQ=255)
                                                                                                                                             ggAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGt     <  1:860255/150‑1 (MQ=255)
                                                                                                                                                aTTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGAc  <  2:345233/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GACACATTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGAC  >  CP009273/2077065‑2077358

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: