Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,028,552 A→T 6.6% T320S (ACT→TCT)  hyaA → hydrogenase 1, small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,028,5520AT 6.6% 902.6 / 42.7 256T320S (ACT→TCT) hyaAhydrogenase 1, small subunit
Reads supporting (aligned to +/- strand):  ref base A (130/108);  new base T (6/11);  total (137/119)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

Only 100 of 257 total aligned reads displayed.

ACCTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAGGA  >  CP009273/1028405‑1028700
                                                                                                                                                   |                                                                                                                                                    
accTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTAc                                                                                                                                                     <  2:113266/149‑1 (MQ=255)
 ccTATAACGCCTGTTCCTCCACACGCTGTAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTTCt                                                                                                                                                    <  1:2700174/149‑1 (MQ=255)
 ccTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTTCt                                                                                                                                                    <  1:2699322/149‑1 (MQ=255)
      aaCGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTc                                                                                                                                               <  2:1112583/149‑1 (MQ=255)
        cGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCa                                                                                                                                             <  2:324040/149‑1 (MQ=255)
         gCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTTCTCATTCCAc                                                                                                                                            <  1:1093691/149‑1 (MQ=255)
             gTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCAccgcc                                                                                                                                        <  2:3020056/149‑1 (MQ=255)
               tcctccACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGa                                                                                                                                      >  1:3848024/1‑149 (MQ=255)
               tcctccACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGa                                                                                                                                      <  2:771648/149‑1 (MQ=255)
                  tccACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATAc                                                                                                                                   <  1:3728738/149‑1 (MQ=255)
                  tccACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATAc                                                                                                                                   >  1:3611941/1‑149 (MQ=255)
                   ccACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATAcc                                                                                                                                  >  1:4223474/1‑149 (MQ=255)
                     acacGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGt                                                                                                                                <  1:1205815/149‑1 (MQ=255)
                     acacGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTACACCGCCGATACCGt                                                                                                                                >  1:3949827/1‑149 (MQ=255)
                          cTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCAGTCCACCGCCGATACCGTCGGtt                                                                                                                           <  2:3904087/149‑1 (MQ=255)
                          cTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCGGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGGGGTCGATATTCCGCAAATGGGTACTCATTCCCCCGCCGATACCGTCgggt                                                                                                                           >  1:2007361/1‑147 (MQ=255)
                            ggAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTa                                                                                                                         <  2:2227580/149‑1 (MQ=255)
                              gctgatgGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCGGTGCGGAAAATGGTTTCGGGGATCGCGGGTCGGTGTACAGCCGCGGGGTCGATATTCCGCAACTGGGTCCTCATTCCACCGCCGCTACCGGCGGGTTAAc                                                                                                                       >  1:2453676/3‑149 (MQ=255)
                              aaTGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAAc                                                                                                                       >  1:4182837/1‑149 (MQ=255)
                                tgatgGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACcg                                                                                                                     >  1:1485646/1‑149 (MQ=255)
                                 gatgGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTTCTCATTCCACCGCCGATACCGTCGGTTTaaaccc                                                                                                                    >  1:3256746/1‑145 (MQ=255)
                                 gatgGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACcgc                                                                                                                    >  1:1532017/1‑149 (MQ=255)
                                   tgGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGAGACCGTCGGTTTAACCgcgc                                                                                                                  <  2:1243655/149‑1 (MQ=255)
                                       gTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTgg                                                                                                              <  1:1817271/149‑1 (MQ=255)
                                       gTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCCCGCATgg                                                                                                              >  2:1643673/1‑149 (MQ=255)
                                        tttctttcCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGc                                                                                                             >  1:3057667/1‑149 (MQ=255)
                                        tttctttcCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGc                                                                                                             >  2:2408508/1‑149 (MQ=255)
                                            tttcCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCgtgg                                                                                                         >  2:4438331/1‑149 (MQ=255)
                                            tttcCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCgtgg                                                                                                         >  1:2136932/1‑149 (MQ=255)
                                              tcCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTATCCGCGCTTGGCgtggtg                                                                                                       <  1:240476/149‑1 (MQ=255)
                                                ccAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGc                                                                                                     >  1:3374733/1‑149 (MQ=255)
                                                ccAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGGGGTGGc                                                                                                     >  1:1037217/1‑149 (MQ=255)
                                                   aTCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGc                                                                                                  >  1:4637043/1‑149 (MQ=255)
                                                       aGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCt                                                                                              >  1:3831235/1‑149 (MQ=255)
                                                       aGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCt                                                                                              >  2:2112451/1‑149 (MQ=255)
                                                       aGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCt                                                                                              >  1:434217/1‑149 (MQ=255)
                                                        gTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTg                                                                                             <  1:3106156/149‑1 (MQ=255)
                                                         tCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGt                                                                                            <  1:3576772/149‑1 (MQ=255)
                                                           tGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTg                                                                                          <  1:2344455/149‑1 (MQ=255)
                                                               cACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGgtgt                                                                                      <  2:1738886/149‑1 (MQ=255)
                                                                 cGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGc                                                                                    >  2:2561546/1‑149 (MQ=255)
                                                                   gCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCAc                                                                                  >  1:702323/1‑149 (MQ=255)
                                                                   gCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCAc                                                                                  <  1:493216/149‑1 (MQ=255)
                                                                     tGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGc                                                                                >  1:1918180/1‑149 (MQ=255)
                                                                     tGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGc                                                                                <  1:4132363/149‑1 (MQ=255)
                                                                       ccTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAg                                                                              <  1:169061/149‑1 (MQ=255)
                                                                       ccTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAg                                                                              <  1:3492576/149‑1 (MQ=255)
                                                                        cTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGt                                                                             <  2:1587634/149‑1 (MQ=255)
                                                                        cTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGt                                                                             >  1:3549427/1‑149 (MQ=255)
                                                                         tGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGGGGTGGCAGCGGCTGTTGGTGTGCACGCAGTc                                                                            >  2:1760487/1‑149 (MQ=255)
                                                                           ggCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGc                                                                          <  1:2237273/149‑1 (MQ=255)
                                                                            gCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGcc                                                                         <  1:3000600/149‑1 (MQ=255)
                                                                            gCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGcc                                                                         <  2:205131/149‑1 (MQ=255)
                                                                              tgtgCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAg                                                                       >  1:1966788/1‑149 (MQ=255)
                                                                               gtgCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAgc                                                                      >  1:610595/1‑149 (MQ=255)
                                                                                      aaaTGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTg                                                               <  2:1499151/149‑1 (MQ=255)
                                                                                        aTGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGAc                                                             >  1:4156458/1‑149 (MQ=255)
                                                                                         tGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGAcc                                                            >  2:1452821/1‑149 (MQ=255)
                                                                                           gTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGAcccg                                                          >  2:2504143/1‑147 (MQ=255)
                                                                                                 gggATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGa                                                    >  1:4593241/1‑149 (MQ=255)
                                                                                                 gggATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGa                                                    <  2:2522636/149‑1 (MQ=255)
                                                                                                 gggATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGa                                                    <  2:3756873/149‑1 (MQ=255)
                                                                                                 gggATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGa                                                    <  1:570445/149‑1 (MQ=255)
                                                                                                 gggATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAAACGCGCATGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCACTCGCCAGCGCCGTTGACCAGCGCAGa                                                    >  1:1361897/1‑149 (MQ=255)
                                                                                                  ggATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGAc                                                   <  2:3605982/149‑1 (MQ=255)
                                                                                                  ggATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAAACGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGAc                                                   >  1:450746/1‑149 (MQ=255)
                                                                                                   gATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACg                                                  >  2:1882199/1‑149 (MQ=255)
                                                                                                    aTCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGt                                                 >  2:4585610/1‑149 (MQ=255)
                                                                                                    aTCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGt                                                 <  1:1844008/149‑1 (MQ=255)
                                                                                                     tCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTc                                                <  2:3086245/149‑1 (MQ=255)
                                                                                                       gcagttcgttcTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAAACGCGCATGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGTCCAGCGCAGACGTCaa                                              >  1:1229774/4‑148 (MQ=255)
                                                                                                          gttcgttcTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAAc                                           >  2:277482/1‑149 (MQ=255)
                                                                                                           ttcgttcTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAAcc                                          <  2:61648/149‑1 (MQ=255)
                                                                                                            tcgttcTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCa                                         <  2:3965993/149‑1 (MQ=255)
                                                                                                             cgttcTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAg                                        <  1:832076/149‑1 (MQ=255)
                                                                                                               ttcTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCCGCGCCGTTGACCCGCGCAGACGTCATAACCAGCa                                      >  1:3807359/1‑149 (MQ=255)
                                                                                                                tcTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCaa                                     <  2:937113/149‑1 (MQ=255)
                                                                                                                tcTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGGTTAACCGCGCGTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCCGCGCCGTTGACCAGCGCCGACGTCATAAACAGCaa                                     >  1:1274854/1‑149 (MQ=255)
                                                                                                                 cTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAAc                                    >  2:2118953/1‑149 (MQ=255)
                                                                                                                 cTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAAc                                    <  2:2233847/149‑1 (MQ=255)
                                                                                                                    cAGCCGCGTGGTCGATAGTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTa                                 <  2:471443/149‑1 (MQ=255)
                                                                                                                     aGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGcaaccaac                                >  2:2833376/1‑149 (MQ=255)
                                                                                                                      gCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACa                               >  2:780034/1‑149 (MQ=255)
                                                                                                                      gCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACa                               >  1:2603855/1‑149 (MQ=255)
                                                                                                                           gTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAAc                          <  2:3140450/149‑1 (MQ=255)
                                                                                                                           gTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAAc                          >  2:3526754/1‑149 (MQ=255)
                                                                                                                               tCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGaa                      >  2:3844240/1‑149 (MQ=255)
                                                                                                                               tCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGaa                      <  2:1597438/149‑1 (MQ=255)
                                                                                                                               tCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGaa                      <  2:4332601/149‑1 (MQ=255)
                                                                                                                                cGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAAc                     <  1:649174/149‑1 (MQ=255)
                                                                                                                                 gATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACa                    <  2:2308865/149‑1 (MQ=255)
                                                                                                                                   tatTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATc                  <  2:4182837/149‑1 (MQ=255)
                                                                                                                                      tCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGc               <  2:4282516/149‑1 (MQ=255)
                                                                                                                                        cGCAAATGGGTTCTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCa             <  2:2185327/149‑1 (MQ=255)
                                                                                                                                        cGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCa             <  2:2469415/149‑1 (MQ=255)
                                                                                                                                         gCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTTACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAg            >  2:1757594/1‑149 (MQ=255)
                                                                                                                                         gCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAg            <  1:4345200/149‑1 (MQ=255)
                                                                                                                                             aTGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCaa        <  2:3200979/149‑1 (MQ=255)
                                                                                                                                              tGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAAt       <  1:3132227/149‑1 (MQ=255)
                                                                                                                                               gggTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACCGAAACCGAACATCAGCCAGGCAATg      >  2:3482745/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACCTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAGGA  >  CP009273/1028405‑1028700

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: