Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,798,001 T→A 18.4% pseudogene (227/1476 nt) arpB → pseudogene, ankyrin repeats

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,798,0010TA18.4% 670.0 / 139.0 250pseudogene (227/1476 nt)arpBpseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base T (100/104);  new base A (19/27);  total (119/131)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 256 total aligned reads displayed.

CGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGG  >  CP009273/1797852‑1798149
                                                                                                                                                     |                                                                                                                                                    
cGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCAtt                                                                                                                                                      >  1:3229780/1‑150 (MQ=255)
  cTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTAc                                                                                                                                                    >  2:1245350/1‑150 (MQ=255)
   tATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACt                                                                                                                                                   >  1:2315817/1‑150 (MQ=255)
   tATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACt                                                                                                                                                   <  2:515649/150‑1 (MQ=255)
     ttttttGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTc                                                                                                                                                 <  2:3242267/150‑1 (MQ=255)
     ttttttGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATAACTTc                                                                                                                                                 <  2:2539660/150‑1 (MQ=255)
      tttttGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCg                                                                                                                                                >  2:3133708/1‑150 (MQ=255)
          tGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATAACTTCGCCAt                                                                                                                                            >  1:3570534/1‑150 (MQ=255)
           gTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCAtc                                                                                                                                           >  1:1756332/1‑150 (MQ=255)
            tCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCAtcc                                                                                                                                          >  1:4068713/1‑150 (MQ=255)
             cGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCAtcct                                                                                                                                         <  1:134746/150‑1 (MQ=255)
                cATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATcctcct                                                                                                                                      >  2:492138/1‑150 (MQ=255)
                cATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATcctcct                                                                                                                                      <  1:1260265/150‑1 (MQ=255)
                cATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCGTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATcctcct                                                                                                                                      <  2:3071086/150‑1 (MQ=255)
                   cTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCa                                                                                                                                   >  2:531561/1‑150 (MQ=255)
                    tGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAg                                                                                                                                  >  2:3346300/1‑150 (MQ=255)
                     gTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAgg                                                                                                                                 >  2:3105206/1‑150 (MQ=255)
                      tCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAggg                                                                                                                                >  2:2582871/1‑150 (MQ=255)
                      tCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAggg                                                                                                                                >  2:287538/1‑150 (MQ=255)
                       cATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGa                                                                                                                               <  2:3491232/150‑1 (MQ=255)
                        aTTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAg                                                                                                                              <  1:3406214/150‑1 (MQ=255)
                          tGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTc                                                                                                                            <  2:3361213/150‑1 (MQ=255)
                          tGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTc                                                                                                                            >  2:1706450/1‑150 (MQ=255)
                           gCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTcg                                                                                                                           <  1:3365423/150‑1 (MQ=255)
                               tACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGtt                                                                                                                       <  2:4068954/150‑1 (MQ=255)
                                 cGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATAACTTCGCCATCCTCCTGCAGGGAGTCGCGTTAt                                                                                                                     >  1:1184921/1‑150 (MQ=255)
                                     cTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATAACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATcgcg                                                                                                                 <  1:2550061/150‑1 (MQ=255)
                                      ttGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGa                                                                                                                >  1:2786153/1‑150 (MQ=255)
                                        gATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGt                                                                                                              >  2:3773567/1‑150 (MQ=255)
                                          ttCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTaa                                                                                                            <  2:2968198/150‑1 (MQ=255)
                                          ttCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTaa                                                                                                            <  2:3210685/150‑1 (MQ=255)
                                           tCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAg                                                                                                           >  1:3763954/1‑150 (MQ=255)
                                               tGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGAtt                                                                                                       <  2:4543520/150‑1 (MQ=255)
                                                 agaAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGAtttt                                                                                                     <  1:2588987/150‑1 (MQ=255)
                                                 agaAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGAtttt                                                                                                     <  1:536392/150‑1 (MQ=255)
                                                 agaAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATAACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGAtttt                                                                                                     <  2:3570534/150‑1 (MQ=255)
                                                   aaaTGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATAACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTgg                                                                                                   <  1:3989698/150‑1 (MQ=255)
                                                      tGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCtt                                                                                                <  2:774108/150‑1 (MQ=255)
                                                      tGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCtt                                                                                                >  1:4094145/1‑150 (MQ=255)
                                                       ggAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTc                                                                                               <  2:868400/150‑1 (MQ=255)
                                                       ggAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATAACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTc                                                                                               <  2:2650639/150‑1 (MQ=255)
                                                           aaGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCtt                                                                                           >  1:4674898/1‑150 (MQ=255)
                                                               cAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGa                                                                                       <  1:4713749/150‑1 (MQ=255)
                                                               cAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGa                                                                                       <  1:893053/150‑1 (MQ=255)
                                                                aaaGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGat                                                                                      <  2:4015606/150‑1 (MQ=255)
                                                                   gTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATAtt                                                                                   <  1:59882/150‑1 (MQ=255)
                                                                    tCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATAttt                                                                                  >  1:3399007/1‑150 (MQ=255)
                                                                      aaCCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTc                                                                                <  1:2793237/150‑1 (MQ=255)
                                                                      aaCCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTc                                                                                <  1:4018496/150‑1 (MQ=255)
                                                                      aaCCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATAACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTc                                                                                >  1:4784418/1‑150 (MQ=255)
                                                                         cTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATAACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCttt                                                                             <  1:753636/150‑1 (MQ=255)
                                                                          tGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTCTCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCtttt                                                                            >  1:4450629/1‑150 (MQ=255)
                                                                          tGATTAAACATTACTCCTCCGTTGAATACATACAGCAGCATGTCCCCTTCGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCtttt                                                                            <  1:4370956/150‑1 (MQ=255)
                                                                             ttAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATattttatttttta                                                                         >  1:4513666/1‑147 (MQ=255)
                                                                              tAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCg                                                                        >  1:2169768/1‑150 (MQ=255)
                                                                               aaaCATTACTCCTCCGTTGAATCCATACAGCAGCAAGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATAACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGa                                                                       >  2:37970/1‑150 (MQ=255)
                                                                                 aCATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAat                                                                     <  1:1706450/150‑1 (MQ=255)
                                                                                   aTTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAatga                                                                   <  1:344866/150‑1 (MQ=255)
                                                                                    ttACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAatgat                                                                  <  2:1049360/150‑1 (MQ=255)
                                                                                      aCTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATgg                                                                <  2:4815904/150‑1 (MQ=255)
                                                                                        tcctccGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCa                                                              >  2:4350823/1‑150 (MQ=255)
                                                                                         cctccGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCaa                                                             <  1:1873548/150‑1 (MQ=255)
                                                                                         cctccGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCaa                                                             <  2:2149826/150‑1 (MQ=255)
                                                                                         cctccGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATAACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTATTTTTTCGAATGATGGCaa                                                             >  1:2020010/1‑150 (MQ=255)
                                                                                              gTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAAt                                                        >  2:4468907/1‑150 (MQ=255)
                                                                                              gTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAAt                                                        >  2:2014268/1‑150 (MQ=255)
                                                                                                tGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATgt                                                      >  2:1340651/1‑150 (MQ=255)
                                                                                                   aTCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATgtcag                                                   <  1:3733982/150‑1 (MQ=255)
                                                                                                      cATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATgtcagtca                                                >  2:3659733/1‑150 (MQ=255)
                                                                                                      cATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATgtcagtca                                                >  2:3661751/1‑150 (MQ=255)
                                                                                                       aTACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATAACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCat                                               <  1:2984528/150‑1 (MQ=255)
                                                                                                             cagcaTGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATAACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTcc                                         <  2:3535400/150‑1 (MQ=255)
                                                                                                                 aTGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTATATATTTTATTTTTTCGAATGATGGCAAACAATGTCAGTCATATGTACGCGa                                     >  1:3283012/1‑150 (MQ=255)
                                                                                                                  tGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGAt                                    <  2:150157/150‑1 (MQ=255)
                                                                                                                   gTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGAtt                                   <  1:1946801/150‑1 (MQ=255)
                                                                                                                      cccTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTta                                <  2:1313872/150‑1 (MQ=255)
                                                                                                                       ccTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTATTTTGTCGAATGATGGCAAACAATGTCAGTCAAATGTGCGCGATTTtat                               >  1:2870389/1‑150 (MQ=255)
                                                                                                                        cTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTtata                              <  1:1322138/150‑1 (MQ=255)
                                                                                                                        cTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCGTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTtata                              >  1:1015172/1‑150 (MQ=255)
                                                                                                                         ttAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTtatat                             <  1:1689078/150‑1 (MQ=255)
                                                                                                                          tAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTtatata                            >  1:1914986/1‑150 (MQ=255)
                                                                                                                           aGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTtatatat                           >  2:2726526/1‑150 (MQ=255)
                                                                                                                             tCCGGGACGCGGAGTTCAGAGCATAACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAg                         <  1:2228466/150‑1 (MQ=255)
                                                                                                                              ccGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGa                        <  1:1856492/150‑1 (MQ=255)
                                                                                                                              ccGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGa                        >  2:4726293/1‑150 (MQ=255)
                                                                                                                               cGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGAc                       <  2:4068713/150‑1 (MQ=255)
                                                                                                                                 ggACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACaa                     >  1:2121254/1‑150 (MQ=255)
                                                                                                                                    cgcgGAGTTCAGAGCATAACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCa                  <  1:1479746/150‑1 (MQ=255)
                                                                                                                                       ggAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATAc               >  2:4439433/1‑150 (MQ=255)
                                                                                                                                         aGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACtt             >  2:2515786/1‑150 (MQ=255)
                                                                                                                                          gTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACttt            >  2:1010766/1‑150 (MQ=255)
                                                                                                                                           ttCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTg           >  2:1859076/1‑150 (MQ=255)
                                                                                                                                           ttCAGAGCATAACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTg           <  1:2926697/150‑1 (MQ=255)
                                                                                                                                            tCATAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGt          >  2:159054/1‑150 (MQ=255)
                                                                                                                                             cAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTc         >  1:3841805/1‑150 (MQ=255)
                                                                                                                                             cAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTc         <  2:1875192/150‑1 (MQ=255)
                                                                                                                                              agagCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCa        <  2:3573958/150‑1 (MQ=255)
                                                                                                                                                   aTTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCATTAAg   >  2:1447412/1‑150 (MQ=255)
                                                                                                                                                   aTTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAg   <  1:3133708/150‑1 (MQ=255)
                                                                                                                                                    ttACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAgg  <  2:78118/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGG  >  CP009273/1797852‑1798149

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: