Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,990,933 G→T 7.4% intergenic (+41/+21) tyrP → / ← yecA tyrosine:H(+) symporter/UPF0149 family protein YecA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,990,9330GT 7.4% 264.7 / 7.6 81intergenic (+41/+21)tyrP/yecAtyrosine:H(+) symporter/UPF0149 family protein YecA
Reads supporting (aligned to +/- strand):  ref base G (34/41);  new base T (2/4);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01

GGTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAG  >  NC_000913/1990787‑1991077
                                                                                                                                                  |                                                                                                                                                
ggTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTg                                                                                                                                               <  2:1076882/150‑1 (MQ=255)
     aaGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACAcc                                                                                                                                          >  2:1242689/1‑150 (MQ=255)
      aGGTGGTCCACCGGAGCGGCTGGTGGTATTTATCTGTGGTATTGATGTGATTGGCGTGAAATGTTTGATTGCGGACGGGTTGATACAAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTCGGTGGCAGAATTAACGCTGACGccc                                                                                                                                         >  2:105070/1‑150 (MQ=255)
           gtcgtcCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGTCAGCGTTAATTCTGGCACCCAAtaa                                                                                                                                    <  2:302982/150‑3 (MQ=255)
            gcgtcCGGCGCTGGGGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGTCAGCGTTAATTCTGGCACCCAAtaag                                                                                                                                   <  2:1043949/149‑4 (MQ=255)
             cgtcCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGTCAGCGTTAATTCTGGCACCCAAtaagg                                                                                                                                  <  2:5059/150‑5 (MQ=39)
             cgtcCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTgg                                                                                                                                  >  2:832641/1‑150 (MQ=255)
                 cGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCctc                                                                                                                              >  1:679870/1‑150 (MQ=255)
                  ggcgggggtggtggtGTTTCTCTGTGGTATTGCTGTTATTGTCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAATTGGGTTTAGCAGTTTTCCTCATATTCCTTATTGTGTGCCTGAATTTAAGCTTACACCCAATTTGGActct                                                                                                                             <  2:163674/144‑1 (MQ=255)
                   gcgcTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCtt                                                                                                                            >  2:961152/1‑150 (MQ=255)
                     gcTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTaa                                                                                                                          <  1:476928/150‑1 (MQ=255)
                     gcTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTaa                                                                                                                          >  1:1472704/1‑150 (MQ=255)
                      cTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGTCAGCGTTAATTCTGGCACCCAAta                                                                                                                                     <  1:845520/138‑2 (MQ=39)
                      cTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGTCAGCGTTAATTCTGGCACCCAAta                                                                                                                                     >  2:845520/1‑137 (MQ=39)
                      cTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAAt                                                                                                                         <  1:433152/150‑1 (MQ=255)
                        ggtggtggtGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACACAATTTGGCCTCTTAATgc                                                                                                                       >  1:1063343/1‑150 (MQ=255)
                         gtggtggtgTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATgca                                                                                                                      <  2:1393405/150‑1 (MQ=255)
                           ggtggtgTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTaa                                                                                                                          >  1:835398/1‑144 (MQ=255)
                           ggtggtgTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTaa                                                                                                                          <  2:835398/144‑1 (MQ=255)
                            gtggtgTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGGTACCAGAAGGGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGACTCTTAATgcaggc                                                                                                                   <  2:31815/150‑1 (MQ=255)
                             gggtgTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATgcaggca                                                                                                                  <  2:1105091/149‑1 (MQ=255)
                              ggtgTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATgcaggcag                                                                                                                 >  1:1094761/1‑150 (MQ=255)
                                tgtTTCTCTGTGGAATTGCTCTGATTGCCCTGCAATATTTGATTGCGGCAGGGTACTTACCAGAAGTGGCGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCAAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGgcagca                                                                                                               >  2:1056563/1‑150 (MQ=255)
                                    tctcTGTGGTATTGCTGTGATTGGTGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGc                                                                                                           >  1:462522/1‑150 (MQ=255)
                                     ctctGTGGTATTGCTGTGATTGGCGTGCAATTTATGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCt                                                                                                          <  2:1118313/150‑1 (MQ=255)
                                        tgtgGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTaa                                                                                                       >  2:1087781/1‑150 (MQ=255)
                                        tgtgGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGACAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTaa                                                                                                       >  2:784532/1‑150 (MQ=255)
                                          tgGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAAt                                                                                                     >  2:305499/1‑150 (MQ=255)
                                                  cTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTAcc                                                                                             >  1:963870/1‑150 (MQ=255)
                                                       aTTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTAc                                                                                        <  1:1354697/150‑1 (MQ=255)
                                                            cGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCAc                                                                                   >  2:1045643/1‑150 (MQ=255)
                                                             gTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACg                                                                                  >  1:235186/1‑150 (MQ=255)
                                                             gTGCAATTTTTGATTGCGGCAGGGTTGTTACAAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACg                                                                                  <  1:1101669/150‑1 (MQ=255)
                                                                   ttgttGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAg                                                                            <  2:1488717/147‑1 (MQ=255)
                                                                      ttGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGAt                                                                         <  1:220716/150‑1 (MQ=255)
                                                                            gCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTAc                                                                                              >  1:855481/1‑123 (MQ=255)
                                                                            gCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTAc                                                                                              <  2:855481/123‑1 (MQ=255)
                                                                            gCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTAc                                                                   >  1:655645/1‑150 (MQ=255)
                                                                             cGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACg                                                                  <  2:462522/150‑1 (MQ=255)
                                                                                 aGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGcccc                                                              <  1:1242689/150‑1 (MQ=255)
                                                                                   ggTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCgg                                                            <  1:1440661/150‑1 (MQ=255)
                                                                                      tgttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGttt                                                         >  2:1586287/1‑150 (MQ=255)
                                                                                       gttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtc                                                        >  2:634316/1‑150 (MQ=255)
                                                                                         tACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtctc                                                      <  1:1087781/150‑1 (MQ=255)
                                                                                            cAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACACAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCtcttc                                                   <  2:1063343/150‑1 (MQ=255)
                                                                                                    gggTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGa                                           >  1:1355899/1‑150 (MQ=255)
                                                                                                     ggTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGAt                                          >  2:965213/1‑150 (MQ=255)
                                                                                                     ggTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGAt                                          <  1:454997/150‑1 (MQ=255)
                                                                                                      gTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATc                                         <  1:837635/150‑1 (MQ=255)
                                                                                                      gTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATc                                         <  1:350932/150‑1 (MQ=255)
                                                                                                        gATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCgg                                       <  1:1243134/150‑1 (MQ=255)
                                                                                                         aTCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGt                                      <  2:1453152/150‑1 (MQ=255)
                                                                                                          tCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGtt                                     <  2:1433668/150‑1 (MQ=255)
                                                                                                                 gCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTATCTTCTGCTTTCATCGGTTGGTGaaa                              >  2:1239305/1‑149 (MQ=255)
                                                                                                                 gCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACAACCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAAc                              >  2:1357251/1‑150 (MQ=255)
                                                                                                                  ccTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTa                                                                    <  1:1626204/111‑1 (MQ=255)
                                                                                                                  ccTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTa                                                                    >  2:1626204/1‑111 (MQ=255)
                                                                                                                   cTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAg                            <  2:1080222/150‑1 (MQ=255)
                                                                                                                      aaaTTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCa                                                                                    >  1:553421/1‑91 (MQ=255)
                                                                                                                      aaaTTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCa                                                                                    <  2:553421/91‑1 (MQ=255)
                                                                                                                           ccTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCTCGGGCAAGGATCGTTACGCCCCGTTTTCTCTTCTTCTTTGATCGGTTGCTGAACAGCTTTTTcc                    <  1:1056563/150‑1 (MQ=255)
                                                                                                                           ccTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc                    <  1:1045643/150‑1 (MQ=255)
                                                                                                                           ccTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc                    <  2:1472704/150‑1 (MQ=255)
                                                                                                                            cTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGcc                                                                <  1:526673/105‑1 (MQ=255)
                                                                                                                            cTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGcc                                                                >  2:526673/1‑105 (MQ=255)
                                                                                                                               aTTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCg                <  2:875339/150‑1 (MQ=255)
                                                                                                                               aTTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCg                <  2:655645/150‑1 (MQ=255)
                                                                                                                                ttGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCgg               >  1:1568133/1‑150 (MQ=255)
                                                                                                                                ttGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCgg               <  1:1586287/150‑1 (MQ=255)
                                                                                                                                 tGGGTGCCAGAATTAACGCTGACACCCAATTTTGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGa              <  2:235186/150‑1 (MQ=255)
                                                                                                                                  gggTGTCAGCGTTAATTCTGGCACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGAt             >  2:644271/1‑150 (MQ=255)
                                                                                                                                  gggTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCCCGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGAt             >  2:438405/1‑150 (MQ=255)
                                                                                                                                  gggTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGAt             >  1:330602/1‑150 (MQ=255)
                                                                                                                                     tGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATgcg          <  1:1103374/150‑1 (MQ=255)
                                                                                                                                     tGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATgcg          >  1:909441/1‑150 (MQ=255)
                                                                                                                                      gCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATgcgc         <  1:832641/150‑1 (MQ=255)
                                                                                                                                        cAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCa       >  1:1597171/1‑150 (MQ=255)
                                                                                                                                          gAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATc     >  2:555262/1‑150 (MQ=255)
                                                                                                                                           aaTTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATcc    <  1:305499/150‑1 (MQ=255)
                                                                                                                                             ttAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAg  <  2:679870/150‑1 (MQ=255)
                                                                                                                                             ttAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAg  <  2:1094761/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAG  >  NC_000913/1990787‑1991077

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