Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 191,739 T→C 8.0% intergenic (+31/‑116) 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,7390TC 8.0% 215.7 / 6.7 75intergenic (+31/‑116)tsf/pyrHprotein chain elongation factor EF‑Ts/UMP kinase
Reads supporting (aligned to +/- strand):  ref base T (31/38);  new base C (2/4);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGC  >  NC_000913/191593‑191887
                                                                                                                                                  |                                                                                                                                                    
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:1055873/150‑1 (MQ=255)
                                    gCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTtgtg                                                                                                               <  2:550388/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                                                                                         <  1:1403232/144‑1 (MQ=255)
                                                                aGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGt                                                                                         >  2:1403232/1‑144 (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:1244457/150‑1 (MQ=255)
                                                                                             ctgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTg                                                      >  2:568917/1‑150 (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                                           <  2:33482/121‑1 (MQ=255)
                                                                                                                                     ccTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAgg                                           >  1:33482/1‑121 (MQ=255)
                                                                                                                                        gAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTc           <  2:294641/150‑1 (MQ=255)
                                                                                                                                        gAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTc           <  1:568917/150‑1 (MQ=255)
                                                                                                                                         caggcggcTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCt          <  1:1193585/148‑1 (MQ=255)
                                                                                                                                         aGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCt          <  1:980234/150‑1 (MQ=255)
                                                                                                                                          gggCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCta         >  1:455831/1‑150 (MQ=255)
                                                                                                                                           ggcggcTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTATCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCtat        >  1:911198/1‑150 (MQ=255)
                                                                                                                                            gcggcTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCtata       >  1:341926/1‑150 (MQ=255)
                                                                                                                                             cggcTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATaa      >  2:670313/1‑150 (MQ=255)
                                                                                                                                                cttcttTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACg   <  2:1086819/150‑1 (MQ=255)
                                                                                                                                                cttcttTTTGTGACCATCTTGTAAATTCATCTAACCCTTCTGGGGCTGCGCTGAAAAGAGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACg   >  2:894727/1‑150 (MQ=255)
                                                                                                                                                 ttcttTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGc  >  2:1423691/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGC  >  NC_000913/191593‑191887

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: