Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 514,339 A→T 5.4% D116V (GAC→GTC)  cueR → DNA‑binding transcriptional dual regulator CueR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913514,3390AT 5.4% 285.3 / 5.1 93D116V (GAC→GTC) cueRDNA‑binding transcriptional dual regulator CueR
Reads supporting (aligned to +/- strand):  ref base A (43/45);  new base T (2/3);  total (45/48)
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.13e-01

GGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGG  >  NC_000913/514190‑514486
                                                                                                                                                     |                                                                                                                                                   
ggTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGa                                                                                                                                                     >  1:247113/1‑150 (MQ=255)
   gAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAg                                                                                                                                                  <  2:414574/150‑1 (MQ=255)
    aaTCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGCTGTCatc                                                                                                                                                 >  1:861966/1‑148 (MQ=255)
     aTCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGCTGTCAtcg                                                                                                                                                >  2:1597421/1‑150 (MQ=255)
     aTCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAgcg                                                                                                                                                <  1:1168352/150‑1 (MQ=255)
        tGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCg                                                                                                                                             <  1:73898/150‑1 (MQ=255)
         gTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGCTGTCATCGccag                                                                                                                                            <  1:889511/150‑3 (MQ=255)
          tttAACGACCCGCAGCGTCACAGCGCCGACGTCAAACGTCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTTCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGCTGTCATCGccagg                                                                                                                                           <  2:1403044/150‑4 (MQ=255)
           ttAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGCTGTCATCGccaggg                                                                                                                                          <  1:1597421/150‑5 (MQ=255)
           ttAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACt                                                                                                                                          >  1:1107350/1‑150 (MQ=255)
            tAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTg                                                                                                                                         >  2:1016025/1‑150 (MQ=255)
            tAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTg                                                                                                                                         <  2:1593281/150‑1 (MQ=255)
            tAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTg                                                                                                                                         >  2:215976/1‑150 (MQ=255)
              aCGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGcc                                                                                                                                       >  2:1151158/1‑150 (MQ=255)
               cGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGccc                                                                                                                                      >  1:750528/1‑150 (MQ=255)
               cGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGccc                                                                                                                                      >  1:1322351/1‑150 (MQ=255)
               cGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGccc                                                                                                                                      >  1:225214/1‑150 (MQ=255)
               cGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGccc                                                                                                                                      >  1:933074/1‑150 (MQ=255)
               cGACCCACAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGccc                                                                                                                                      >  1:525676/1‑150 (MQ=255)
                  cccGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGAt                                                                                                                                   <  2:247113/150‑1 (MQ=255)
                    cGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTa                                                                                                                                 <  1:1143044/150‑1 (MQ=255)
                         cGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGaa                                                                                                                            >  2:942678/1‑150 (MQ=255)
                         cGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGaa                                                                                                                            <  1:290664/150‑1 (MQ=255)
                           gCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGaaaa                                                                                                                          >  2:861654/1‑150 (MQ=255)
                            cacaGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAt                                                                                                                         >  2:35326/1‑150 (MQ=255)
                             acaGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAtc                                                                                                                        <  2:925237/150‑1 (MQ=255)
                                  gcCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTcc                                                                                                                   <  2:440549/150‑1 (MQ=255)
                                   ccGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCg                                                                                                                  >  2:936806/1‑150 (MQ=255)
                                    cGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCgg                                                                                                                 <  2:525676/150‑1 (MQ=255)
                                     gACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGgc                                                                                                                <  1:268732/150‑1 (MQ=255)
                                     gACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGgc                                                                                                                >  1:443300/1‑150 (MQ=255)
                                         tCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGgctgct                                                                                                            >  1:686714/1‑150 (MQ=255)
                                          cAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGctgctg                                                                                                           >  1:749910/1‑150 (MQ=255)
                                           aaaCGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGt                                                                                                          >  1:771868/1‑150 (MQ=255)
                                            aaCGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGtc                                                                                                         >  1:511010/1‑150 (MQ=255)
                                                     cGCTGGAGAAGGTGGCGGAGATCGAACGATACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCggg                                                                                                <  2:1606177/150‑1 (MQ=255)
                                                      gCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGc                                                                                               <  1:400792/150‑1 (MQ=255)
                                                              aGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGAt                                                                                       >  2:1117359/1‑150 (MQ=255)
                                                              aGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGAt                                                                                       >  1:1090336/1‑150 (MQ=255)
                                                               ggTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGAtt                                                                                      <  1:215976/150‑1 (MQ=255)
                                                               ggTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGAtt                                                                                      <  1:1016025/150‑1 (MQ=255)
                                                                  ggCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAag                                                                                   >  2:331057/1‑150 (MQ=255)
                                                                  ggCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAag                                                                                   <  2:750528/150‑1 (MQ=255)
                                                                  ggCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAag                                                                                   >  1:1073715/1‑150 (MQ=255)
                                                                    cGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagac                                                                                 >  2:193820/1‑150 (MQ=255)
                                                                      gagaTCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacga                                                                               <  2:497750/150‑1 (MQ=255)
                                                                       agaTCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacgag                                                                              <  1:861654/150‑1 (MQ=255)
                                                                         aTCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacgagac                                                                            <  1:169018/150‑1 (MQ=255)
                                                                               cGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCg                                                                      <  1:779845/150‑1 (MQ=255)
                                                                                 acacATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGa                                                                    >  1:698194/1‑150 (MQ=255)
                                                                                   acaTTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGata                                                                  <  1:1126803/150‑1 (MQ=255)
                                                                                   acaTTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGata                                                                  <  1:1132459/150‑1 (MQ=255)
                                                                                    caTTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGatat                                                                 <  2:749910/150‑1 (MQ=255)
                                                                                    caTTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCCCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTGGGatat                                                                 >  2:1430787/1‑150 (MQ=255)
                                                                                      ttGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGt                                                               <  2:283764/150‑1 (MQ=255)
                                                                                       tGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTa                                                              >  2:1261918/1‑150 (MQ=255)
                                                                                        gaggagCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCAGCGGGCAGGGTGATTAAGACGAGCCGGCTCGGATATGTAg                                                             >  2:1088780/1‑150 (MQ=255)
                                                                                         aggagCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAgg                                                            <  2:1489613/150‑1 (MQ=255)
                                                                                           gagCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGt                                                          <  1:1261918/150‑1 (MQ=255)
                                                                                           gagCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGt                                                          <  2:58941/150‑1 (MQ=255)
                                                                                             gCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTa                                                        >  2:530520/1‑150 (MQ=255)
                                                                                                gCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATcc                                                     >  1:796722/1‑150 (MQ=255)
                                                                                                   aTCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCtt                                                  >  2:1279380/1‑150 (MQ=255)
                                                                                                   aTCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCtt                                                  >  1:286201/1‑150 (MQ=255)
                                                                                                   aTCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCtt                                                  <  1:217546/150‑1 (MQ=255)
                                                                                                     ccATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCt                                                <  1:484845/150‑1 (MQ=255)
                                                                                                     ccATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCt                                                <  2:511010/150‑1 (MQ=255)
                                                                                                      cATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTa                                               >  2:527091/1‑150 (MQ=255)
                                                                                                             gACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAAt                                        <  1:1520002/150‑1 (MQ=255)
                                                                                                                 aGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACt                                    <  2:276809/150‑1 (MQ=255)
                                                                                                                  gctgctGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACtt                                   >  1:1597391/1‑150 (MQ=255)
                                                                                                                    tgctgGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCa                                 >  1:264693/1‑150 (MQ=255)
                                                                                                                    tgctgGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCa                                 <  1:189512/150‑1 (MQ=255)
                                                                                                                      ctgGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAAc                               <  2:1322351/150‑1 (MQ=255)
                                                                                                                       tgGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACa                              <  2:1039287/150‑1 (MQ=255)
                                                                                                                        ggCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACAt                             <  1:1279380/150‑1 (MQ=255)
                                                                                                                         gCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATg                            >  1:1199476/1‑150 (MQ=255)
                                                                                                                          cACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATg                                       <  1:678795/138‑1 (MQ=255)
                                                                                                                          cACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATg                                       >  2:504347/1‑138 (MQ=255)
                                                                                                                          cACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATg                                       >  2:678795/1‑138 (MQ=255)
                                                                                                                          cACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATg                                       <  1:504347/138‑1 (MQ=255)
                                                                                                                          cACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGc                           <  1:331057/150‑1 (MQ=255)
                                                                                                                           aCTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCg                          >  1:1550330/1‑150 (MQ=255)
                                                                                                                           aCTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCg                          <  1:528839/150‑1 (MQ=255)
                                                                                                                           aCTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCg                          >  2:704702/1‑150 (MQ=255)
                                                                                                                                     gcctgccCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCAcc                <  1:402388/150‑1 (MQ=255)
                                                                                                                                      cctgccCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCg               >  1:75998/1‑150 (MQ=255)
                                                                                                                                      cctgccCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCg               >  2:1097915/1‑150 (MQ=255)
                                                                                                                                      cctgccCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCg               <  1:980402/150‑1 (MQ=255)
                                                                                                                                       ctgccCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCga              <  2:483258/150‑1 (MQ=255)
                                                                                                                                         gccCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCgaga            >  2:1186256/1‑150 (MQ=255)
                                                                                                                                             tGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCt        >  1:498158/1‑150 (MQ=255)
                                                                                                                                               gCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCg      <  2:933074/150‑1 (MQ=255)
                                                                                                                                                   tGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTgg  <  1:510451/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGG  >  NC_000913/514190‑514486

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