Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,578,609 T→G 6.4% I40L (ATT→CTT)  gntR ← DNA‑binding transcriptional repressor GntR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,578,6090TG 6.4% 205.8 / 6.5 84I40L (ATT→CTT) gntRDNA‑binding transcriptional repressor GntR
Reads supporting (aligned to +/- strand):  ref base T (42/36);  new base G (3/3);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-02

TTTCGATTCCGCGTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAAGAGA  >  NC_000913/3578468‑3578748
                                                                                                                                             |                                                                                                                                           
tttCGATTCCGCGTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATc                                                                                                                                            >  2:882809/1‑143 (MQ=255)
tttCGATTCCGCGTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCACAGTTCATCAGGAGCCGCGGCAATc                                                                                                                                            >  1:137353/1‑143 (MQ=255)
 ttCGATTCCGCGTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCt                                                                                                                                           >  2:351990/1‑143 (MQ=255)
    gATTCCGCGTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCACTCTCGc                                                                                                                                        >  1:140199/1‑143 (MQ=255)
    gATTCCGCGTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGc                                                                                                                                        >  1:319605/1‑143 (MQ=255)
     aTTCCGCGTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGcc                                                                                                                                       <  1:1221194/143‑1 (MQ=255)
      ctCCGCGTAATACTCCCCCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCcg                                                                                                                                      <  2:716311/142‑1 (MQ=255)
        ccGCGTAATACTTCCGCGAAAACCTGGGTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCcgcg                                                                                                                                    <  2:789037/143‑1 (MQ=255)
         cgcgTAGTACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAAGATTGCcgcgg                                                                                                                                   <  2:1136733/143‑2 (MQ=255)
         cgcgTAATTCTTCCGCCAAAACCTGCTTGGTGAGAGACGGTAACAGGACGCCAATCGCCCGGCTGGTTTCGTTAGCTATGTTATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAAGATTGCcgcgg                                                                                                                                   <  2:1037106/143‑2 (MQ=255)
         cgcgTAATGCGTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAAGATTGCcgcgg                                                                                                                                   <  1:316385/143‑2 (MQ=255)
         cgcgTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGCCGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCGAGAGCCGCGGCAATCTTGCCGCGt                                                                                                                                   >  2:532250/1‑143 (MQ=255)
         cgcgTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGt                                                                                                                                   >  1:592018/1‑143 (MQ=255)
          gcgTAATACTTCCGCGAAAAAAAAAATCTTTAGTTACGGTCACACGTCGCTATTCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTa                                                                                                                                  <  1:1139557/143‑1 (MQ=255)
                 agttCCGCGAAATCCTGGTTGGTGAGAGAAGGTAACAGGACGCCAGTCGCCCGGCTGGGGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGAc                                                                                                                           <  2:715882/141‑1 (MQ=255)
                 aCTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGAc                                                                                                                           <  1:1266540/143‑1 (MQ=255)
                  cgtCCGGGAAACCCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACg                                                                                                                          <  2:366837/141‑1 (MQ=255)
                   ttCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGTCGCCTAGAGCGACgg                                                                                                                         >  2:468362/1‑143 (MQ=255)
                       gcgAGAACGTGGTGTGTGAGCGAAGGTACCAGGACGCAGGTCGCCGAGATGGTGGCGTTAAAGAGGCTATCGGGCGCACGATTAGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGAGTAGAGCGACGGAAAc                                                                                                                     <  2:601051/143‑1 (MQ=255)
                       gcgAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGGAGAGCGACGGAAAc                                                                                                                     >  1:359024/1‑143 (MQ=255)
                          aaaaCCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTg                                                                                                                  >  1:376561/1‑143 (MQ=255)
                          aaaaCCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTg                                                                                                                  >  1:89773/1‑143 (MQ=255)
                           aagcgtggttgggGAGCGAAGGTAACAGGACGCCAATCGCCCGGCTCGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGc                                                                                                                 <  1:1155636/138‑1 (MQ=255)
                           aacccTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGc                                                                                                                 <  1:1143317/140‑1 (MQ=255)
                           aaaCCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCGTCAAGAGCCGCGGCAATCTTGCCGGGTAGATGGACGGAAACCTGc                                                                                                                 >  1:470365/1‑143 (MQ=255)
                             aCCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTc                                                                                                               <  2:19394/143‑1 (MQ=255)
                                tggttggtGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCgg                                                                                                            >  1:892040/1‑143 (MQ=255)
                                          gagaAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTaca                                                                                                  >  1:1258032/1‑141 (MQ=255)
                                          gagaAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGCCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAAGCTTGCCGAGTAGCGCAACGGAAACCGACTCCGGGTGGCGTaaa                                                                                                  >  2:160225/1‑143 (MQ=255)
                                             aaGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATCGGGAATATAGCCCAGTTCGTCAAGAGCCGCGGCAATCGTGGCGCGTAGAGCGTCGGAAACCTGCTCCGGGTCGCGTtaacca                                                                                               >  1:210142/1‑140 (MQ=255)
                                              gtgTAGCAGGCCGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCAGCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAAc                                                                                              <  2:1265330/141‑1 (MQ=255)
                                               ggTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCGGGGTTGCGTAAGAAACg                                                                                             >  1:85878/1‑143 (MQ=255)
                                                 gaaCGGGACGGCAACCGCCCGGCGGGTGGCGTGAGAGAGGAGATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATGTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGc                                                                                           <  1:401956/142‑1 (MQ=255)
                                                  aaCAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCt                                                                                          >  1:346640/1‑143 (MQ=255)
                                                  aaCAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCt                                                                                          <  1:882809/143‑1 (MQ=255)
                                                          gCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCa                                                                                  >  1:976205/1‑143 (MQ=255)
                                                            cAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCAtt                                                                                >  1:130467/1‑143 (MQ=255)
                                                              accGCCCGGCTGGTGGCGTTAGAAGGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCAtttt                                                                              <  2:231616/141‑1 (MQ=255)
                                                              aTCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCAtttt                                                                              >  1:282426/1‑143 (MQ=255)
                                                              aTCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCAttct                                                                              >  2:6147/1‑141 (MQ=255)
                                                                  cccGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGCCCTTTTGGTc                                                                          >  2:498830/1‑143 (MQ=255)
                                                                  cccGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGATGCGTAAAAAACGGCTGACGGTCATTTTGGTc                                                                          >  1:1168703/1‑143 (MQ=255)
                                                                   ccGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCa                                                                         >  1:125084/1‑143 (MQ=255)
                                                                     ggCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACg                                                                       >  1:95481/1‑143 (MQ=255)
                                                                         cttcgCGTTAGCGCCGATAGCGGGCGCATGATTGGCAACATAGGCCTGTTCATCAAGGGCCGCGGCAATCTTGCCGCGTGGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGTCTa                                                                   <  2:400709/139‑1 (MQ=255)
                                                                                     gagGATGTCGGGCGCACGAGTGAGAAGATAGCTCAGTTCACCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCGCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCca                                                       <  1:168333/143‑1 (MQ=255)
                                                                                     gagGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCca                                                       >  1:966156/1‑143 (MQ=255)
                                                                                          tatCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATcc                                                  <  2:346640/143‑1 (MQ=255)
                                                                                            tCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGGAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACAGCCTg                                                >  2:446967/1‑143 (MQ=255)
                                                                                              gggCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGaa                                              >  2:165986/1‑143 (MQ=255)
                                                                                               ggCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTggac                                             >  1:692422/1‑140 (MQ=255)
                                                                                                  gcACGGTTGGGAACAGACCCCAGTCCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTAc                                          <  1:503103/143‑1 (MQ=255)
                                                                                                  gcACGATTGGGAATATAGCCCAGGTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTAc                                          <  1:505216/143‑1 (MQ=255)
                                                                                                     cGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACggg                                       <  2:89773/143‑1 (MQ=255)
                                                                                                     cGATGGGGAATATAGCCCAGTTCAGCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACggg                                       <  1:738265/143‑1 (MQ=255)
                                                                                                      gATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGt                                      <  2:592018/143‑1 (MQ=255)
                                                                                                      gATTGGGAATATAGCCCAGTCCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGt                                      <  2:382779/143‑1 (MQ=255)
                                                                                                         tGGGAATATAGCCCAGTTCATCAAGGGGCGCGGCAATCTTGCCGCGGAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGtctt                                   <  1:171812/143‑1 (MQ=255)
                                                                                                           ggAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGtctttt                                 <  2:846705/143‑1 (MQ=255)
                                                                                                              atatAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGtcttttctt                              >  1:1171935/1‑143 (MQ=255)
                                                                                                                ataGCCCAGTTCCTCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTcttttctttt                            <  2:176113/143‑1 (MQ=255)
                                                                                                                ataGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTTCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTcttttctttt                            <  1:79081/143‑1 (MQ=255)
                                                                                                                 taGCCCAGTTCATCAAGAGCCGCGGCAATCGTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCATAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCttttt                           <  1:1009563/143‑1 (MQ=255)
                                                                                                                    cccAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCGGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCAt                        >  1:1155600/1‑143 (MQ=255)
                                                                                                                    cccAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCAt                        >  1:1111955/1‑143 (MQ=255)
                                                                                                                    cccAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTATTTTTCAt                        >  1:832910/1‑143 (MQ=255)
                                                                                                                     ccAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATc                       >  2:307255/1‑143 (MQ=255)
                                                                                                                         ttCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTcc                   >  1:900316/1‑143 (MQ=255)
                                                                                                                            atcaAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCGCCGGGTTGCGTAAAAAACGGATGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGACGTACGGGCCTTTTCTTTGTCCTCGTCCTGa                >  2:975318/1‑143 (MQ=255)
                                                                                                                             tcaAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGaa               <  2:19438/143‑1 (MQ=255)
                                                                                                                             tcaAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGaa               <  1:1021899/143‑1 (MQ=255)
                                                                                                                              caAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAg              >  1:822466/1‑143 (MQ=255)
                                                                                                                               aaGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAgg             >  1:1054369/1‑143 (MQ=255)
                                                                                                                                 tctacgcgGCAAGATTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTa           >  2:1010600/5‑143 (MQ=255)
                                                                                                                                 gagCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTa           <  2:892040/143‑1 (MQ=255)
                                                                                                                                   tacgcgGCAAGATTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGGCCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACa         >  2:735300/3‑143 (MQ=255)
                                                                                                                                    ccGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACaa        >  1:214572/1‑143 (MQ=255)
                                                                                                                                    ccGCGGCAATCTTGCCGCGCAGAGCGACGGAAACCTGCTCAGGGTTGCGTAAAAAGCGGCTGACCGTCATTGGGGTCACGGCTACACGGTCAGCCGCACCCTGCAGGACGGGGCTTTCCGTCTTCATCGCCCCGAAggcgaaa        >  2:945634/1‑138 (MQ=255)
                                                                                                                                     cgcgGCAATCTTGCCGCGTAGAGCGCCGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGCCCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTGTTTCAGCGGCCTGGAGGTGCaga       >  1:161041/1‑141 (MQ=255)
                                                                                                                                      gcgGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACaaaa      <  2:966156/143‑1 (MQ=255)
                                                                                                                                        ggCGATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAaga    <  2:47084/143‑1 (MQ=255)
                                                                                                                                        ggCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAaga    <  2:290429/143‑1 (MQ=255)
                                                                                                                                         gCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAagag   >  1:241791/1‑143 (MQ=255)
                                                                                                                                          cAATCTGGCCGCGTAGAGCGACGGAAACCGGCTCCGGGGTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAAgaga  <  2:633687/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                           
TTTCGATTCCGCGTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAAGAGA  >  NC_000913/3578468‑3578748

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