Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,798,112 C→A 8.6% pseudogene (338/1476 nt) arpB → pseudogene, ankyrin repeats

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,798,1120CA 8.6% 1060.4 / 51.8 290pseudogene (338/1476 nt)arpBpseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base C (125/140);  new base A (12/13);  total (137/153)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 298 total aligned reads displayed.

GCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCCACAGGACAGCACGGGGCTTTC  >  CP009273/1797963‑1798261
                                                                                                                                                     |                                                                                                                                                     
gcaTGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCcgc                                                                                                                                                       >  1:3933012/1‑150 (MQ=255)
   tGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGAt                                                                                                                                                    <  2:3908456/150‑1 (MQ=255)
     tCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGAttt                                                                                                                                                  <  1:4039506/150‑1 (MQ=255)
        ccTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTtat                                                                                                                                               <  1:1415685/150‑1 (MQ=255)
        ccTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGAGATTTtat                                                                                                                                               <  2:1840914/150‑1 (MQ=255)
         cTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTtata                                                                                                                                              <  1:4140789/150‑1 (MQ=255)
         cTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTtata                                                                                                                                              >  2:4051098/1‑150 (MQ=255)
          ttAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTtatat                                                                                                                                             >  1:2197324/1‑150 (MQ=255)
            aGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTCAGTATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTtatatct                                                                                                                                           >  2:3341724/1‑148 (MQ=255)
              tCCGGGACGCGGATTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAg                                                                                                                                         >  1:4667846/1‑150 (MQ=255)
                cGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGAc                                                                                                                                       <  1:188123/150‑1 (MQ=255)
                cGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGAc                                                                                                                                       >  1:265392/1‑150 (MQ=255)
                cGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGAc                                                                                                                                       <  1:2569943/150‑1 (MQ=255)
                 gggACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACa                                                                                                                                      >  1:4479884/1‑150 (MQ=255)
                 gggACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACa                                                                                                                                      >  1:3475536/1‑150 (MQ=255)
                 gggACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACa                                                                                                                                      >  1:2354155/1‑150 (MQ=255)
                 gggACGCGGAGTTCAGAGCATTACGTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACa                                                                                                                                      <  2:4848751/150‑1 (MQ=255)
                    aCGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGTATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAAATATGTCAGTCATATGTCCGCGATTTTATATATAGACAATc                                                                                                                                   <  2:1531894/150‑1 (MQ=255)
                    aCGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATc                                                                                                                                   <  2:1931450/150‑1 (MQ=255)
                        ggAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATAc                                                                                                                               >  1:2302896/1‑150 (MQ=255)
                        ggAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATAc                                                                                                                               <  2:2573011/150‑1 (MQ=255)
                        ggAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATAc                                                                                                                               >  2:2795721/1‑150 (MQ=255)
                        ggAGTTCAGAGCATTACTTCGCCATCCTCCGGCAGGGAGTCGCGTTCTCGCGAGTAAGGATTCTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATAc                                                                                                                               >  2:5144064/1‑150 (MQ=255)
                         gAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACt                                                                                                                              <  2:1334503/150‑1 (MQ=255)
                         gAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACt                                                                                                                              <  1:3399643/150‑1 (MQ=255)
                           gTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCTAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACttt                                                                                                                            <  1:183708/150‑1 (MQ=255)
                           gTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACttt                                                                                                                            <  1:5026980/150‑1 (MQ=255)
                           gTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACttt                                                                                                                            <  1:5027099/150‑1 (MQ=255)
                            ttCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGTATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTg                                                                                                                           >  1:2053554/1‑150 (MQ=255)
                              cAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGAGATTTTATATATAGACAATCATACTTTGTc                                                                                                                         <  2:4961093/150‑1 (MQ=255)
                              cAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGAGATTTTATATATAGACAATCATACTTTGTc                                                                                                                         >  2:2617349/1‑150 (MQ=255)
                                gagCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATAATTTGTCAg                                                                                                                       >  1:2905101/1‑150 (MQ=255)
                                  gCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTa                                                                                                                     <  2:3475433/150‑1 (MQ=255)
                                   cATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTaa                                                                                                                    <  1:2036072/150‑1 (MQ=255)
                                    aTTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATAATTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAg                                                                                                                   >  2:2591455/1‑150 (MQ=255)
                                     ttACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAgg                                                                                                                  <  1:2951558/150‑1 (MQ=255)
                                     ttACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATAATTTGTCAGTAAgg                                                                                                                  >  1:4890419/1‑150 (MQ=255)
                                      tACTTCGCCATCCTCCTCCAGGGAGTCGCGTTATCGCGAGTAAGGATATTGGCATCGCATTAGATATTTTCATTTGTCGAATGATGGCAAACAATGTCAGTCAAAAGTCCCCGAATTAATAAATAGACAAACATACTTAGTCAGTATGGc                                                                                                                 >  2:784972/1‑150 (MQ=255)
                                       aCTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCt                                                                                                                <  2:1790433/150‑1 (MQ=255)
                                         ttCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTAc                                                                                                              >  2:390260/1‑150 (MQ=255)
                                           cGCCATCCTCCTGCAGGGTGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACNa                                                                                                            <  1:3884633/150‑1 (MQ=255)
                                            gCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGaa                                                                                                           <  1:1858403/150‑1 (MQ=255)
                                             ccATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGTATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAAt                                                                                                          >  2:4778622/1‑150 (MQ=255)
                                                tcctccTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATaaa                                                                                                       <  2:4440040/150‑1 (MQ=255)
                                                tcctccTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATaaa                                                                                                       >  2:670283/1‑150 (MQ=255)
                                                tcctccTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATaaa                                                                                                       <  1:1181234/150‑1 (MQ=255)
                                                 cctcctGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAg                                                                                                      >  1:3749879/1‑150 (MQ=255)
                                                   tcctGCATGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGca                                                                                                    <  2:3322788/150‑1 (MQ=255)
                                                      tGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGcagtc                                                                                                 <  2:1140819/150‑1 (MQ=255)
                                                      tGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGcagtc                                                                                                 <  1:1328734/150‑1 (MQ=255)
                                                       gCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGcagtca                                                                                                <  1:3349454/150‑1 (MQ=255)
                                                         aGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGAGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGcagtcagt                                                                                              >  2:3068475/1‑150 (MQ=255)
                                                           ggAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAgtgt                                                                                            >  2:680953/1‑150 (MQ=255)
                                                             aGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTAt                                                                                          <  1:1160351/150‑1 (MQ=255)
                                                              gTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATg                                                                                         <  2:385534/150‑1 (MQ=255)
                                                                cgcgTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGa                                                                                       <  2:881165/150‑1 (MQ=255)
                                                                 gcgTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGAc                                                                                      <  1:437132/150‑1 (MQ=255)
                                                                  cgTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACa                                                                                     <  1:242333/150‑1 (MQ=255)
                                                                          cgAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACaattaaatt                                                                             >  2:1665302/1‑150 (MQ=255)
                                                                            aGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTAt                                                                           >  1:3034484/1‑150 (MQ=255)
                                                                             gTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTAtg                                                                          <  1:4562098/150‑1 (MQ=255)
                                                                                 ggATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTAtgttgt                                                                      >  2:2283805/1‑150 (MQ=255)
                                                                                  gATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTtgtg                                                                     >  1:4846635/1‑150 (MQ=255)
                                                                                     tttGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTcc                                                                  <  2:3265091/150‑1 (MQ=255)
                                                                                      ttGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCg                                                                 <  1:43444/150‑1 (MQ=255)
                                                                                           ttCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGAGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTAc                                                            <  2:4906011/150‑1 (MQ=255)
                                                                                               cTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGAcc                                                        <  1:728488/150‑1 (MQ=255)
                                                                                               cTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATAGAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGAcc                                                        <  2:2404498/150‑1 (MQ=255)
                                                                                                tttAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGAccc                                                       >  1:4234613/1‑150 (MQ=255)
                                                                                                    gATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAAc                                                   <  1:1233057/150‑1 (MQ=255)
                                                                                                       atTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAGCGAt                                                >  1:868183/1‑150 (MQ=255)
                                                                                                       atTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGAGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGAt                                                >  1:57419/1‑150 (MQ=255)
                                                                                                         tttCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATAc                                              <  1:3865231/150‑1 (MQ=255)
                                                                                                           tCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCa                                            >  1:1699312/1‑150 (MQ=255)
                                                                                                               ttGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGt                                        >  2:1432298/1‑150 (MQ=255)
                                                                                                                tGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGtt                                       >  2:3734280/1‑150 (MQ=255)
                                                                                                                 gTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTg                                      <  2:1435011/150‑1 (MQ=255)
                                                                                                                  tCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGc                                     >  2:4667672/1‑150 (MQ=255)
                                                                                                                   cGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGcc                                    <  2:2496675/150‑1 (MQ=255)
                                                                                                                   cGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGcc                                    <  1:1737945/150‑1 (MQ=255)
                                                                                                                         atgGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCaga                              >  2:1273058/1‑150 (MQ=255)
                                                                                                                         atgGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCaga                              <  2:2090512/150‑1 (MQ=255)
                                                                                                                              aaaCAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGacac                         >  1:1566460/1‑150 (MQ=255)
                                                                                                                                aCAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACAccc                       <  2:2325347/150‑1 (MQ=255)
                                                                                                                                  aaTGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCcac                     >  2:1190690/1‑150 (MQ=255)
                                                                                                                                  aaTGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCcac                     >  1:4306867/1‑150 (MQ=255)
                                                                                                                                  aaTGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCcac                     >  1:1947341/1‑150 (MQ=255)
                                                                                                                                  aaTGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCcac                     >  1:787709/1‑150 (MQ=255)
                                                                                                                                   aTGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCcaca                    <  1:652691/150‑1 (MQ=255)
                                                                                                                                   aTGTCAGTCATATGTCCGCGATTTTATATATAGAAAAAAATACTTTGTCAGAAAGGCAAATAATAAAACACTCAGTGTATGGACAAATAAATTATGTTGTGACCGTAAACGACCCGAACGATACCAACATTGCAGACAGAGAcacacaca                    >  2:594655/1‑150 (MQ=255)
                                                                                                                                          tcaTATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCCACAGGACAGc             <  2:1720007/150‑1 (MQ=255)
                                                                                                                                            atatGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCCACAGGACAGCAc           >  1:2722094/1‑150 (MQ=255)
                                                                                                                                            atatGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCCACAGGACAGCAc           <  1:989783/150‑1 (MQ=255)
                                                                                                                                             tatGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCCACAGGACAGCACg          >  2:84964/1‑150 (MQ=255)
                                                                                                                                              atGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCCACAGGACAGCACgg         >  1:3275147/1‑150 (MQ=255)
                                                                                                                                                gTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCCACAGGACAGCACgggg       <  2:4729708/150‑1 (MQ=255)
                                                                                                                                                 tCCGAGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCCACAGGACAGCACGGGGc      <  2:628130/150‑1 (MQ=255)
                                                                                                                                                   cgcgATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCCACAGGACAGCACGGGGctt    <  2:2666978/150‑1 (MQ=255)
                                                                                                                                                    gcgATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCCACAGGACAGCACGGGGcttt   <  2:426469/150‑1 (MQ=255)
                                                                                                                                                     cgATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCCACAGGACAGCACGGGGctttc  >  1:258951/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTCATATGTCCGCGATTTTATATATAGACAATCATACTTTGTCAGTAAGGCTACGAATAAAGCAGTCAGTGTATGGACAATTAAATTATGTTGTGTCCGTTTACGACCCGAACGATACCAACGTTGCCGTCAGAGACACCCACAGGACAGCACGGGGCTTTC  >  CP009273/1797963‑1798261

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: