Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 191,731 G→A 7.0% intergenic (+23/‑124) tsf → / → pyrH protein chain elongation factor EF‑Ts/UMP kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913191,7310GA 7.0% 217.7 / 4.7 71intergenic (+23/‑124)tsf/pyrHprotein chain elongation factor EF‑Ts/UMP kinase
Reads supporting (aligned to +/- strand):  ref base G (27/39);  new base A (2/3);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

TGAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTA  >  NC_000913/191593‑191880
                                                                                                                                          |                                                                                                                                                     
tgAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGcttctt                                                                                                                                            <  1:124001/150‑1 (MQ=255)
   aaGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGAAGCCGCCCTCAGGCGGCTCCttttt                                                                                                                                         <  1:329063/150‑1 (MQ=255)
    agagCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTtg                                                                                                                                        <  2:1117834/150‑1 (MQ=255)
      agCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTtgtg                                                                                                                                      <  1:76962/150‑1 (MQ=255)
      agCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGAAGCCGCCCTCAGGCGGCTCCTTTTtgat                                                                                                                                      <  2:477787/150‑3 (MQ=39)
       gCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGc                                                                                                                                     >  1:661542/1‑150 (MQ=255)
        cATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGcc                                                                                                                                    <  2:267732/150‑1 (MQ=255)
        cATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGcc                                                                                                                                    <  1:323201/150‑1 (MQ=255)
         aTAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGAAGCCGCCCTCAGGCGGCTCCTTTTtgataat                                                                                                                                   <  2:43281/150‑6 (MQ=38)
          tAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCa                                                                                                                                  <  2:387403/150‑1 (MQ=255)
             cGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCt                                                                                                                               >  2:546257/1‑150 (MQ=255)
                    gTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAAt                                                                                                                        >  2:437021/1‑150 (MQ=255)
                      gACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTc                                                                                                                      <  1:521542/150‑1 (MQ=255)
                         tGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGc                                                                                                                   <  1:437021/150‑1 (MQ=255)
                          ggCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCt                                                                                                                  <  1:811222/150‑1 (MQ=255)
                              tCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAAcc                                                                                                              >  1:1061904/1‑150 (MQ=255)
                                    gCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTtgtg                                                                                                        <  2:550388/150‑1 (MQ=255)
                                    gCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTtgtg                                                                                                        <  2:1055873/150‑1 (MQ=255)
                                        cGAAGTGGGTGAAGGCATCGAGAAAGTTGATACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGc                                                                                                    <  2:1213165/150‑1 (MQ=255)
                                        cGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGc                                                                                                    <  1:1274251/150‑1 (MQ=255)
                                         gAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCt                                                                                                   <  2:1318911/150‑1 (MQ=255)
                                         gAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCt                                                                                                   <  2:300733/150‑1 (MQ=255)
                                          aaGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTg                                                                                                  >  1:329777/1‑150 (MQ=255)
                                          aaGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGATAACCCTTGTGGGGGTg                                                                                                  >  1:324586/1‑150 (MQ=255)
                                           aGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTgc                                                                                                 >  1:1217450/1‑150 (MQ=255)
                                              gggTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCt                                                                                              <  2:1382286/150‑1 (MQ=255)
                                                gTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGa                                                                                            >  1:508501/1‑150 (MQ=255)
                                                    aGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAg                                                                                        >  2:907957/1‑150 (MQ=255)
                                                       cATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGa                                                                                     >  1:1177957/1‑150 (MQ=255)
                                                         tCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACg                                                                                   >  1:1258552/1‑150 (MQ=255)
                                                          cGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGt                                                                                  >  2:111993/1‑150 (MQ=255)
                                                            tgaAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTAc                                                                                <  2:209968/149‑1 (MQ=255)
                                                              aaaGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACaa                                                                              <  1:1389796/150‑1 (MQ=255)
                                                               aaGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAAt                                                                             >  1:10990/1‑150 (MQ=255)
                                                                aGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGt                                                                                  >  2:1403232/1‑144 (MQ=255)
                                                                aGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGt                                                                                  <  1:1403232/144‑1 (MQ=255)
                                                                     agaCTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCGTGAGGGCTGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTTCGCTGAAAAGCGACGTACAATGTCGCt                                                                       <  2:154927/150‑1 (MQ=255)
                                                                          gactTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTAt                                                                  >  1:346529/1‑150 (MQ=255)
                                                                              ttGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAAt                                                              <  1:1235658/150‑1 (MQ=255)
                                                                               tGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAAtt                                                             <  1:1249006/150‑1 (MQ=255)
                                                                                 cagcagAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCa                                                           >  1:1227623/1‑150 (MQ=255)
                                                                                 cagcagAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCa                                                           <  1:1082518/150‑1 (MQ=255)
                                                                                  agcagAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCAt                                                          >  2:255757/1‑150 (MQ=255)
                                                                                   gcagAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCAtt                                                         <  2:909560/150‑1 (MQ=255)
                                                                                   gcagAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCAtt                                                         >  1:1086819/1‑150 (MQ=255)
                                                                                    cagAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCAttt                                                        <  1:488214/150‑1 (MQ=255)
                                                                                          ttGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCg                                                  <  2:1217450/150‑1 (MQ=255)
                                                                                             ctgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTa                                                                   <  2:736072/130‑1 (MQ=255)
                                                                                             ctgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTa                                                                   >  1:736072/1‑130 (MQ=255)
                                                                                             ctgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTg                                               >  2:568917/1‑150 (MQ=255)
                                                                                             ctgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTg                                               <  2:1244457/150‑1 (MQ=255)
                                                                                               gcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGAc                                             >  2:83859/1‑150 (MQ=255)
                                                                                                cGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACa                                            >  2:459149/1‑150 (MQ=255)
                                                                                                 gATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAg                                           <  2:31227/150‑1 (MQ=255)
                                                                                                      ccAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCa                                      <  2:735764/150‑1 (MQ=255)
                                                                                                            aGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaag                                <  2:1061904/150‑1 (MQ=255)
                                                                                                             gTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaaga                               <  1:546257/150‑1 (MQ=255)
                                                                                                               cTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaagaaa                             >  1:815852/1‑150 (MQ=255)
                                                                                                                 tAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACa                           >  2:422692/1‑150 (MQ=255)
                                                                                                                 tAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACa                           >  2:577561/1‑150 (MQ=255)
                                                                                                                     gcacAAAAAGAAGCCGCCCTCAGGCGGCTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGc                       >  1:445970/4‑150 (MQ=39)
                                                                                                                        cAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAc                    >  2:437307/1‑150 (MQ=255)
                                                                                                                        cAAAAAGAAGCCGCCCTCAGGCGGCTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAc                    >  1:1015100/1‑150 (MQ=255)
                                                                                                                             aGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCTCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATg               <  1:418964/150‑1 (MQ=255)
                                                                                                                             aGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATg               <  1:907957/150‑1 (MQ=255)
                                                                                                                                 gccgccTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCaaa           <  2:858030/150‑1 (MQ=255)
                                                                                                                                     ccTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAgg                                    >  1:33482/1‑121 (MQ=255)
                                                                                                                                     ccTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAgg                                    <  2:33482/121‑1 (MQ=255)
                                                                                                                                        gAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTc    <  1:568917/150‑1 (MQ=255)
                                                                                                                                        gAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTc    <  2:294641/150‑1 (MQ=255)
                                                                                                                                         aGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCt   <  1:980234/150‑1 (MQ=255)
                                                                                                                                          gggCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCta  >  1:455831/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                     
TGAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTA  >  NC_000913/191593‑191880

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