Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,788,346 C→A 29.1% T7T (ACG→ACT)  btuE ← glutathione peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,788,3460CA29.1% 601.0 / 204.1 306T7T (ACG→ACT) btuEglutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base C (106/110);  new base A (45/45);  total (151/155)
Fisher's exact test for biased strand distribution p-value = 9.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 317 total aligned reads displayed.

CCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCA  >  CP009273/1788197‑1788495
                                                                                                                                                     |                                                                                                                                                     
   cgatcgaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGgtcgtc                                                                                                                                                    <  1:1430226/150‑1 (MQ=255)
   cgatcgaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGgtcgtc                                                                                                                                                    >  1:928051/1‑150 (MQ=255)
    gatcgaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCa                                                                                                                                                   >  2:4113956/1‑150 (MQ=255)
    gatcgaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCa                                                                                                                                                   <  2:4221919/150‑1 (MQ=255)
    gatcgaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCa                                                                                                                                                   >  1:2782840/1‑150 (MQ=255)
     atcgaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAg                                                                                                                                                  <  1:3618524/150‑1 (MQ=255)
      tcgaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGa                                                                                                                                                 <  1:2394937/150‑1 (MQ=255)
      tcgaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGa                                                                                                                                                 >  2:1541006/1‑150 (MQ=255)
        gaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAAt                                                                                                                                               <  1:3976965/150‑1 (MQ=255)
        gaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGAGATCTTTCACTACGGTCGTCAGaaa                                                                                                                                               >  1:4950823/1‑149 (MQ=255)
           ccAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGa                                                                                                                                            >  1:1686308/1‑150 (MQ=255)
            cAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGaa                                                                                                                                           >  1:4741022/1‑150 (MQ=255)
            cAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGaa                                                                                                                                           >  2:5243717/1‑150 (MQ=255)
            cAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGaa                                                                                                                                           >  2:2623503/1‑150 (MQ=255)
               gCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCt                                                                                                                                        >  1:3766143/1‑150 (MQ=255)
                 tttCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTg                                                                                                                                      >  2:1984943/1‑150 (MQ=255)
                  ttCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGc                                                                                                                                     >  2:2540242/1‑150 (MQ=255)
                    cTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCAt                                                                                                                                   >  2:214624/1‑150 (MQ=255)
                     tGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATc                                                                                                                                  >  2:430929/1‑150 (MQ=255)
                         tatTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGttt                                                                                                                              <  1:3187200/150‑1 (MQ=255)
                           ttCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTtct                                                                                                                            <  1:1539999/150‑1 (MQ=255)
                               ccAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTg                                                                                                                        <  1:1537306/150‑1 (MQ=255)
                               ccAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTg                                                                                                                        <  1:2478380/150‑1 (MQ=255)
                                  aCTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGgt                                                                                                                     <  2:1924699/150‑1 (MQ=255)
                                   cTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGgtg                                                                                                                    >  2:4980987/1‑150 (MQ=255)
                                     gCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTACTGGGtgtg                                                                                                                  >  1:2273227/1‑150 (MQ=255)
                                      cTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTgg                                                                                                                 >  1:4647446/1‑150 (MQ=255)
                                         atatTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCa                                                                                                              >  2:3996162/1‑150 (MQ=255)
                                           atTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGt                                                                                                            <  1:2123762/150‑1 (MQ=255)
                                               cggcgTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTaaaa                                                                                                        <  2:1087349/150‑1 (MQ=255)
                                                 gcgTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAAt                                                                                                      <  2:2044415/150‑1 (MQ=255)
                                                   gTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCt                                                                                                    <  1:5281747/150‑1 (MQ=255)
                                                    ttAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCtt                                                                                                   >  2:319961/1‑150 (MQ=255)
                                                     tAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACATCACCGTCGATATCTTTCACTACGGTCGTCAGAAAGGAAACTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCtaa                                                                                                  >  1:3844102/1‑148 (MQ=255)
                                                        gCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCt                                                                                               <  2:2722419/150‑1 (MQ=255)
                                                               tttGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATc                                                                                        <  1:1263365/150‑1 (MQ=255)
                                                               tttGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATc                                                                                        <  1:3310351/150‑1 (MQ=255)
                                                                  gAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATa                                                                                     >  2:3445780/1‑150 (MQ=255)
                                                                  gAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCAATACGGTAGTCAGAATGGAATCTTGCATCGTTTCTACTGGGTTTGGTCAGTAAAAAACTTATCTTTTAATCATa                                                                                     >  2:3737811/1‑150 (MQ=255)
                                                                   aGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAg                                                                                    <  2:244671/150‑1 (MQ=255)
                                                                   aGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAg                                                                                    <  2:1505454/150‑1 (MQ=255)
                                                                   aGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAg                                                                                    <  2:1405824/150‑1 (MQ=255)
                                                                   aGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAg                                                                                    <  2:1974702/150‑1 (MQ=255)
                                                                     gCGACAGTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGAc                                                                                  <  1:180504/150‑1 (MQ=255)
                                                                       gACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCg                                                                                <  1:4431222/150‑1 (MQ=255)
                                                                         cATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTc                                                                              >  2:2324665/1‑150 (MQ=255)
                                                                          aTTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCt                                                                             >  1:4012347/1‑150 (MQ=255)
                                                                          aTTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCt                                                                             >  1:4848011/1‑150 (MQ=255)
                                                                            tGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCttt                                                                           <  2:783245/150‑1 (MQ=255)
                                                                             gACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCtttt                                                                          >  1:1752366/1‑150 (MQ=255)
                                                                             gACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCtttt                                                                          <  1:2975727/150‑1 (MQ=255)
                                                                              aCAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAATAATCTTAGCTTTTAATCATAGACCGTCttttt                                                                         >  2:1199759/1‑150 (MQ=255)
                                                                              aCAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCttttt                                                                         >  2:905500/1‑150 (MQ=255)
                                                                              aCAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCttttt                                                                         <  1:4988800/150‑1 (MQ=255)
                                                                               cAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTg                                                                        <  2:2926008/150‑1 (MQ=255)
                                                                               cAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTg                                                                        <  1:3272741/150‑1 (MQ=255)
                                                                                aaTCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATATTATCTTTTAATCATAGACCGTCTTTTTGc                                                                       >  1:5131102/1‑150 (MQ=255)
                                                                                aaTCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGc                                                                       >  2:3955580/1‑150 (MQ=255)
                                                                                 aTCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCg                                                                      <  1:1854612/150‑1 (MQ=255)
                                                                                 aTCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCg                                                                      <  1:577758/150‑1 (MQ=255)
                                                                                     aCAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTa                                                                  >  2:655359/1‑150 (MQ=255)
                                                                                       aGCACATTACCGGCGAACTTCTCCAGCGTTGTCACTTCACCGTCGATATCTTTCACGACGGTCGTCAGAATTGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAAc                                                                <  1:4481001/150‑1 (MQ=255)
                                                                                       aGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATTGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAAc                                                                <  1:2310189/150‑1 (MQ=255)
                                                                                       aGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAAc                                                                >  2:3463896/1‑150 (MQ=255)
                                                                                       aGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAAc                                                                >  2:68236/1‑150 (MQ=255)
                                                                                       aGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAAc                                                                >  1:1150179/1‑150 (MQ=255)
                                                                                       aGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAAc                                                                <  1:1401421/150‑1 (MQ=255)
                                                                                           caTTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTcc                                                            <  2:3922082/150‑1 (MQ=255)
                                                                                           caTTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTcc                                                            >  2:1327781/1‑150 (MQ=255)
                                                                                            aTTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCt                                                           <  1:4354504/150‑1 (MQ=255)
                                                                                            aTTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCt                                                           <  2:4016261/150‑1 (MQ=255)
                                                                                             ttACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTg                                                          <  2:590667/150‑1 (MQ=255)
                                                                                                 cGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCttt                                                      <  2:17754/150‑1 (MQ=255)
                                                                                                  ggCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCtttt                                                     <  1:3860552/150‑1 (MQ=255)
                                                                                                     gAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAAc                                                  >  2:4942414/1‑150 (MQ=255)
                                                                                                         ttCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATa                                              <  1:4559994/150‑1 (MQ=255)
                                                                                                           ctcCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAAc                                            >  2:2361103/1‑150 (MQ=255)
                                                                                                            tcCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAAcc                                           >  1:199904/1‑150 (MQ=255)
                                                                                                            tcCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAAcc                                           <  2:1589280/150‑1 (MQ=255)
                                                                                                              cAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAg                                         >  1:3342076/1‑150 (MQ=255)
                                                                                                                gCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGAt                                       >  1:724689/1‑150 (MQ=255)
                                                                                                                  gTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATaa                                     >  2:3099547/1‑150 (MQ=255)
                                                                                                                   tGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATaaa                                    <  2:1657823/150‑1 (MQ=255)
                                                                                                                    ggTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAac                                   <  1:1909872/150‑1 (MQ=255)
                                                                                                                     gTCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAaca                                  <  1:4055387/150‑1 (MQ=255)
                                                                                                                      tCACTTCACCGTCGATATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTATCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAacac                                 >  1:3584513/1‑150 (MQ=255)
                                                                                                                      tCACTTCACCGTCGAGATCTTTCACTACGGTAGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAacac                                 <  2:5075873/150‑1 (MQ=255)
                                                                                                                      tCACTTCACCATCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAacac                                 >  2:4519927/1‑150 (MQ=255)
                                                                                                                         cTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCgg                              <  1:2368686/150‑1 (MQ=255)
                                                                                                                          ttCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                             <  1:2730893/150‑1 (MQ=255)
                                                                                                                            cACCGTCGATATCTTTCAATACGGTAGTCAGAATTGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGTCTAACTTCCTGCTTTTAACAATAACCAGATAAACACCGgcgc                           <  2:3764570/150‑1 (MQ=255)
                                                                                                                                  cGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTaa                     >  2:21752/1‑150 (MQ=255)
                                                                                                                                      atCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGtt                 >  2:4406941/1‑150 (MQ=255)
                                                                                                                                      atCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGtt                 >  1:5200058/1‑150 (MQ=255)
                                                                                                                                            cACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTg           >  2:4662664/1‑150 (MQ=255)
                                                                                                                                               tACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGAcc        >  2:5125534/1‑150 (MQ=255)
                                                                                                                                               tACGGTCGTCAGAATGGAATCTTGCATAGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGAcc        <  2:5160246/150‑1 (MQ=255)
                                                                                                                                                aCGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCa       >  2:2483372/1‑150 (MQ=255)
                                                                                                                                                  ggTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACg     >  1:3556086/1‑150 (MQ=255)
                                                                                                                                                     cgtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCa  >  1:3660355/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCA  >  CP009273/1788197‑1788495

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: ‑