Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,578,610 C→A 6.2% K39N (AAG→AAT 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,6100CA 6.2% 210.6 / 5.9 82K39N (AAG→AATgntRDNA‑binding transcriptional repressor GntR
Reads supporting (aligned to +/- strand):  ref base C (41/35);  new base A (2/3);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-01

TTTCGATTCCGCGTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAAGAGATAGA  >  NC_000913/3578468‑3578752
                                                                                                                                              |                                                                                                                                              
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)
                                                                                                                                              cTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAAgagaatct  <  1:701248/143‑5 (MQ=255)
                                                                                                                                              |                                                                                                                                              
TTTCGATTCCGCGTAATACTTCCGCGAAAACCTGGTTGGTGAGAGAAGGTAACAGGACGCCAATCGCCCGGCTGGTGGCGTTAGAGAGGATATCGGGCGCACGATTGGGAATATAGCCCAGTTCATCAAGAGCCGCGGCAATCTTGCCGCGTAGAGCGACGGAAACCTGCTCCGGGTTGCGTAAAAAACGGCTGACCGTCATTTTGGTCACGCCTACACGGTCAGCCACATCCTGAAGTACGGGTCTTTTCTTTTTCATCGTCCTGAAGGTACAAAAGAGATAGA  >  NC_000913/3578468‑3578752

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: