Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,849,854 A→T 7.7% H27L (CAC→CTC)  yidI → inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,849,8540AT 7.7% 830.3 / 44.6 234H27L (CAC→CTC) yidIinner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (118/98);  new base T (7/11);  total (125/109)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 235 total aligned reads displayed.

AGCGGTTTTTTTAAGAATTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGAT  >  CP009273/3849705‑3850002
                                                                                                                                                     |                                                                                                                                                    
 gCGGTTTTTTTAAGAATTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCAc                                                                                                                                                     >  2:4290997/1‑150 (MQ=255)
 gCGGTTTTTTTAAGAATTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCAc                                                                                                                                                     <  1:390589/150‑1 (MQ=255)
   ggTTTTTTTAAGAATTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCTCtt                                                                                                                                                   >  2:742646/1‑150 (MQ=255)
    gTTTTTTTAAGAATTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCTCttt                                                                                                                                                  <  1:4641217/150‑1 (MQ=255)
    gTTTTTTTAAGAATTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACttt                                                                                                                                                  <  2:2295124/150‑1 (MQ=255)
     tttttttAAGAATTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACtttt                                                                                                                                                 <  2:2869243/150‑1 (MQ=255)
         tttAAGAATTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCAtt                                                                                                                                             <  2:4368718/150‑1 (MQ=255)
         tttAAGAATTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCAtt                                                                                                                                             <  1:2798808/150‑1 (MQ=255)
          ttAAGAATTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTc                                                                                                                                            >  2:2172379/1‑150 (MQ=255)
          ttAAGAATTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTc                                                                                                                                            >  2:4554758/1‑150 (MQ=255)
              gAATTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGcc                                                                                                                                        <  2:3049770/150‑1 (MQ=255)
               aaTTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGccc                                                                                                                                       >  2:347824/1‑150 (MQ=255)
               aaTTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGccc                                                                                                                                       >  1:470247/1‑150 (MQ=255)
               aaTTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGccc                                                                                                                                       >  2:4675743/1‑150 (MQ=255)
               aaTTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGccc                                                                                                                                       >  2:2719384/1‑150 (MQ=255)
                 tttCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCtt                                                                                                                                     >  2:3208030/1‑150 (MQ=255)
                 tttCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCtt                                                                                                                                     >  1:3988697/1‑150 (MQ=255)
                 tttCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTGTTCATTCGGCCCtt                                                                                                                                     >  1:3988690/1‑150 (MQ=255)
                   tCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTtc                                                                                                                                   <  2:3524978/150‑1 (MQ=255)
                     gTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTtctc                                                                                                                                 >  2:3646400/1‑150 (MQ=255)
                     gTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTtctc                                                                                                                                 >  1:4503802/1‑150 (MQ=255)
                      tCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTcc                                                                                                                                >  1:4173141/1‑150 (MQ=255)
                      tCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTcc                                                                                                                                >  1:39511/1‑150 (MQ=255)
                      tCCGCCCGTTTCCTGAGATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTcc                                                                                                                                <  1:3730159/150‑1 (MQ=255)
                       ccGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCg                                                                                                                               >  1:3512075/1‑150 (MQ=255)
                       ccGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCg                                                                                                                               >  1:1119887/1‑150 (MQ=255)
                         gCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGca                                                                                                                             >  1:3057797/1‑150 (MQ=255)
                          cccGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGcac                                                                                                                            <  2:3721710/150‑1 (MQ=255)
                           ccGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGcacc                                                                                                                           >  1:4671967/1‑150 (MQ=255)
                             gTTTCCTGATATGCTTTTGCCATCCGCTACAATGAGGTAAGATGGGTAACATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCTGTGCCATCGCTTTAATGATGTGGATTAACAACTTTTCATTCGGCCCTTTCTCCGcaccac                                                                                                                         >  1:3500421/1‑150 (MQ=255)
                               ttCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCg                                                                                                                       <  1:1723035/150‑1 (MQ=255)
                                  cTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCa                                                                                                                    >  2:4394870/1‑150 (MQ=255)
                                          cTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGa                                                                                                            <  1:41005/150‑1 (MQ=255)
                                               gCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCa                                                                                                       >  1:4689475/1‑150 (MQ=255)
                                                ccATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCAt                                                                                                      <  2:2495641/150‑1 (MQ=255)
                                                 cATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATc                                                                                                     >  1:3183686/1‑150 (MQ=255)
                                                        tACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCgg                                                                                              >  2:4062289/1‑150 (MQ=255)
                                                             ggaggTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACaaa                                                                                         <  2:24042/150‑1 (MQ=255)
                                                              gaggTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACaaaa                                                                                        <  2:2687135/150‑1 (MQ=255)
                                                                ggTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACg                                                                                      <  2:2211411/150‑1 (MQ=255)
                                                                   cAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCTCTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCt                                                                                   <  2:875000/150‑1 (MQ=255)
                                                                       tGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGaaa                                                                               >  1:4334464/1‑150 (MQ=255)
                                                                         ggTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAAtt                                                                             <  1:1352094/150‑1 (MQ=255)
                                                                           tATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTaa                                                                           >  1:228629/1‑150 (MQ=255)
                                                                            aTCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCTCTTTTCATTAGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTaaa                                                                          >  1:4020195/1‑150 (MQ=255)
                                                                            aTCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTaaa                                                                          <  2:3297678/150‑1 (MQ=255)
                                                                            aTCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTaaa                                                                          <  2:3296778/150‑1 (MQ=255)
                                                                             tCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTaaaa                                                                         <  2:4592686/150‑1 (MQ=255)
                                                                              cATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAag                                                                        <  2:3670402/150‑1 (MQ=255)
                                                                              cATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAag                                                                        >  1:1787614/1‑150 (MQ=255)
                                                                                  gCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGt                                                                    <  1:4268283/150‑1 (MQ=255)
                                                                                   cTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGtc                                                                   >  1:198444/1‑150 (MQ=255)
                                                                                   cTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGtc                                                                   >  1:4368041/1‑150 (MQ=255)
                                                                                    tCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGtct                                                                  >  2:1819376/1‑150 (MQ=255)
                                                                                    tCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGtct                                                                  <  2:433910/150‑1 (MQ=255)
                                                                                    tCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGtct                                                                  <  2:4570503/150‑1 (MQ=255)
                                                                                     cAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGtctt                                                                 <  1:1416143/150‑1 (MQ=255)
                                                                                     cAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGtctt                                                                 <  1:3808912/150‑1 (MQ=255)
                                                                                      aaaataaaatCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGtcttc                                                                <  2:1392858/150‑1 (MQ=255)
                                                                                      aaaataaaatCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGtcttc                                                                >  1:2019322/1‑150 (MQ=255)
                                                                                      aaaataaaatCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGtcttc                                                                <  2:3784150/150‑1 (MQ=255)
                                                                                       aaataaaatCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCTCTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGtcttct                                                               <  2:2903154/150‑1 (MQ=255)
                                                                                        aataaaatCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTtctc                                                              <  1:3340678/150‑1 (MQ=255)
                                                                                              atCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCTCTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGa                                                        >  1:1088326/1‑150 (MQ=255)
                                                                                              atCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGa                                                        <  2:341768/150‑1 (MQ=255)
                                                                                                 aGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATa                                                     >  2:2015429/1‑150 (MQ=255)
                                                                                                  gTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATaa                                                    <  2:3232481/150‑1 (MQ=255)
                                                                                                   ttCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATaaa                                                   <  2:3025709/150‑1 (MQ=255)
                                                                                                   ttCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATaaa                                                   <  1:2422294/150‑1 (MQ=255)
                                                                                                   ttCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGAATATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAAAGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATGAAAAGAGGTCTTCTCGCCGGAATaaa                                                   >  2:4125214/1‑150 (MQ=255)
                                                                                                     cATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAg                                                 <  2:276843/150‑1 (MQ=255)
                                                                                                     cATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAg                                                 >  2:1802711/1‑150 (MQ=255)
                                                                                                     cATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAg                                                 >  1:478333/1‑150 (MQ=255)
                                                                                                      aTTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCTCTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAgg                                                <  1:288773/150‑1 (MQ=255)
                                                                                                       ttGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGc                                               >  1:4532089/1‑150 (MQ=255)
                                                                                                          ggCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGaa                                            <  1:1643688/150‑1 (MQ=255)
                                                                                                            cATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGTCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGaaaa                                          <  2:888669/150‑1 (MQ=255)
                                                                                                            cATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGaaaa                                          <  2:1127540/150‑1 (MQ=255)
                                                                                                               gCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAAt                                       <  1:2678996/150‑1 (MQ=255)
                                                                                                                cTGTTCGGTGCCATCGCTTTACTGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAAtt                                      <  2:3820414/150‑1 (MQ=255)
                                                                                                                    tCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTacca                                  <  2:3863366/150‑1 (MQ=255)
                                                                                                                      ggTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTaccacc                                <  1:81813/150‑1 (MQ=255)
                                                                                                                      ggTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTaccacc                                >  1:1283687/1‑150 (MQ=255)
                                                                                                                          ccATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACACAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTgg                            <  1:3646400/150‑1 (MQ=255)
                                                                                                                          ccATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTgg                            <  1:4572999/150‑1 (MQ=255)
                                                                                                                             tCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGgaga                         <  1:910835/150‑1 (MQ=255)
                                                                                                                             tCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGgaga                         <  2:470247/150‑1 (MQ=255)
                                                                                                                             tCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGATGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGgaga                         <  2:4272293/150‑1 (MQ=255)
                                                                                                                               gCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGaaa                       >  1:4299671/1‑150 (MQ=255)
                                                                                                                                   tAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGaacaatg                   >  2:3843864/1‑148 (MQ=255)
                                                                                                                                       gatgGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGACATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGa               >  1:3459059/1‑150 (MQ=255)
                                                                                                                                       gatgGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAACTTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAAGTACCACCGTGGAGACAAAAGAGGa               >  2:3087448/1‑150 (MQ=255)
                                                                                                                                          ggggATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGt            >  1:2921923/1‑150 (MQ=255)
                                                                                                                                           gggATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAGGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTa           <  1:1077561/150‑1 (MQ=255)
                                                                                                                                             gATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGa         >  1:625351/1‑150 (MQ=255)
                                                                                                                                             gATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGa         >  1:2994588/1‑150 (MQ=255)
                                                                                                                                             gATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGa         <  2:156392/150‑1 (MQ=255)
                                                                                                                                               ttATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATg       >  2:178444/1‑150 (MQ=255)
                                                                                                                                               ttATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATg       >  2:1967337/1‑150 (MQ=255)
                                                                                                                                                  tCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCg    <  2:4268786/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCGGTTTTTTTAAGAATTTCGTCCGCCCGTTTCCTGATATGCTTTTGCCATCCGCTACAAGGAGGTCAGATGGGTATCATAGCTCAAAATAAAATCAGTTCATTGGGCATGCTGTTCGGTGCCATCGCTTTAATGATGGGGATTATCCACTTTTCATTCGGCCCTTTCTCCGCACCACCGCCAACGTTTGAAAGCATCGTGGCGGACAAAACGGCTGAAATTAAAAGAGGTCTTCTCGCCGGAATAAAAGGCGAAAAAATTACCACCGTGGAGAAAAAAGAGGACGTAGATGTCGAT  >  CP009273/3849705‑3850002

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: