Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,639,764 G→A S181S (TCC→TCT arcA ← DNA‑binding transcriptional dual regulator ArcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,639,7640GA87.9% 136.1 / 2.3 59S181S (TCC→TCTarcADNA‑binding transcriptional dual regulator ArcA
Reads supporting (aligned to +/- strand):  ref base G (6/1);  new base A (27/24);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTG  >  NC_000913/4639533‑4640001
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
ttGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCCGTTCAGCACGAGACTGAATTTTGCCTgggg                                                                                                                                                                                                                             >  1:454547/1‑250 (MQ=255)
 tGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGGGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCCGTTCAGCACGAGACTGAATTTTGCCTggggt                                                                                                                                                                                                                            >  1:321264/1‑249 (MQ=255)
    aaaaTAAAAACGGCGCTAAAAAGAGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTCCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTGCTTCAGCAGTTCAGCTCGAGACTGAATTTTGCCTGGGTTTTc                                                                                                                                                                                                                         <  2:458224/251‑1 (MQ=255)
          aaaaCGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCAAGATCACCACAGAAGCGCTAACCTTCACCGTGGATGGTGGCGATGATTTCAGGCGTATCCGGCGTAGATTCGAGATGTTTACGAATACGGCGGATACTCAAGTCTAAAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGGCCTTTTTTTCCGCCGTTCAGCACGAGACTGAATTTTGCCTGGGGTTTCCCAGac                                                                                                                                                                                                                   >  2:214464/1‑250 (MQ=255)
                gcgcTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGTATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAg                                                                                                                                                                                                             <  1:428747/251‑1 (MQ=255)
                gcgcTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGa                                                                                                                                                                                                               >  1:437562/1‑249 (MQ=255)
                gcgcTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAg                                                                                                                                                                                                             <  2:25784/251‑1 (MQ=255)
                      aaaaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGGGTGGTTTCAGCTCCCGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATggg                                                                                                                                                                                                      >  1:45983/1‑250 (MQ=255)
                       aaaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGgcg                                                                                                                                                                                                     >  1:253364/1‑251 (MQ=255)
                            gcCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGaa                                                                                                                                                                                                <  1:350759/251‑1 (MQ=255)
                                   tttttGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGTAGCGATAACCTTCAACGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGGCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCATCAGTTCAACACGAGACTGTATTTTGCCTGGGTTTTTCCAGAAGTGAAGACTTGGGCGGAACTCGCTg                                                                                                                                                                                         >  2:373147/1‑251 (MQ=255)
                                   tttttGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTg                                                                                                                                                                                         >  1:344047/1‑251 (MQ=255)
                                    ttttGACGGTGGTAAAGCCGATTAAACTTCCAGATCACCACATAAGCGAGAACCTTCACCGTGGATTGTGGCGATGATTGCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTGCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGAGGAACTCGCTg                                                                                                                                                                                         <  2:107003/250‑1 (MQ=255)
                                                       gATTAATCTTCCAGATCAACACAGAAGCGTTAACCTTCACCGTGGATGGTGGCGATGGTTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATGGTCACGTCTACATTTATGTCGTGTGGGTTTATCTAACGGCAGGGCATTTTCTTTAGAAGTTCAGCACGGGAGTGGATTTTGCCTGGGTTTTCACAGATGTGAGTCATAGCGCGAAACtagcggacaggaagtttgcatataac                                                                                                                                                                     >  2:55671/1‑226 (MQ=255)
                                                            aTCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAt                                                                                                                                                                >  1:365049/1‑251 (MQ=255)
                                                                tCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCCGTTCAGCACGAGACTGAATTTTGCCTGGGGTTTCACAGAAGTGAAGAATGGCGCGGGACTCGCTGcgcgaggttggaaaggcgtcgtgtagggg                                                                                                                                                            >  2:480796/1‑226 (MQ=255)
                                                                tCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGcgcg                                                                                                                                                                                     <  1:480796/226‑1 (MQ=255)
                                                                 ccAGATCCCCACAGTAGCGATATCCTTCACCGTGGATGGGGGCGATGATTTCCGGCGTATTCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCGTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAgt                                                                                                                                                            <  2:13888/250‑2 (MQ=255)
                                                                 ccAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGACGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGCCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTTATTTTCTTCAGCAGTTCAGCACGCGACTGCATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGCACTCGCTGCGCGGCAGCTTGTACTGCTCGCCACCAggg                                                                                                                                                           >  1:224873/1‑251 (MQ=255)
                                                                           acAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCaacg                                                                                                                                                 >  1:255342/2‑251 (MQ=255)
                                                                                  cGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACt                                                                                                                                                                                                                                            >  1:77666/1‑154 (MQ=255)
                                                                                  cGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACt                                                                                                                                                                                                                                            <  2:77666/154‑1 (MQ=255)
                                                                                                gtggaTGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCCCGGCCGGGCATTTTCTTCAGCCGTTCCGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCCGCTTGTACTGCTCGCCATCAGGGCCGATCAACGGACGGCTggtgttggcccgtt                                                                                                                            >  2:41848/5‑238 (MQ=255)
                                                                                                  ggttggtggCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCTCGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTcc                                                                                                                          <  2:95570/248‑1 (MQ=255)
                                                                                                    agggtggCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATTTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCa                                                                                                                        <  2:396318/249‑1 (MQ=255)
                                                                                                       gtggCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAAcc                                                                                                                     <  2:255342/251‑1 (MQ=255)
                                                                                                       gtggCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAAcc                                                                                                                     >  1:265554/1‑251 (MQ=255)
                                                                                                        tggCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGGCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGAGGAACTCGCTGCGCGGCAGCTTGTAATGATCGCCATCAGGGCCGATCAAAGAAAGGCTGTTGATGTCCAGTTCCCAAcc                                                                                                                     >  1:413326/1‑250 (MQ=255)
                                                                                                             atgatTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGa                                                                                                                >  1:77227/1‑250 (MQ=255)
                                                                                                                    ccggcgtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAg                                                                                                        <  2:158772/251‑1 (MQ=255)
                                                                                                                      ggcgtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGGCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCCCAGAAGTGACACATGGCGCGGAACTCGCTGCGCGGCCGCTTGTACTGCTCGCCCTCAGGGCCGATCAACGGACGGCTGTTTATGTCCAGTTCCCAACCATTTAACATTTAGCt                                                                                                      >  2:141309/1‑251 (MQ=255)
                                                                                                                         gtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTc                                                                                                   >  1:213409/1‑251 (MQ=255)
                                                                                                                                 cgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGATCAGACCGGGACTGAATTTTGCCTGGGGTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCCTCAGGGCCGATCCACGAACGGCGGTGTATGTCCCGTTCCCAACCCTTGAACTTGTAGCTTTCCACGCtcc                                                                                           >  2:16055/1‑249 (MQ=255)
                                                                                                                                         cGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTACACGCTacgaaggtct                                                                                   >  1:17192/1‑245 (MQ=255)
                                                                                                                                             aTGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGGGGTTTGAGCTAACCGACGGTCATTTTCTTCAGCACTTCCGCAAGAGACTGAATTTTGCCTGGGTTTTAACAGAAGTGAAGCATGGCGCGGAAATCGCAGAGCGGCAGCTTGTACTGCTAGCCATCAGGGCCGATCAACGAACGGCTGTTGAGGTCCAGTTCCCAACAATTGGACTGGTAGATTTCAACACTACGCGGttctttgc                                                                               >  1:470216/1‑248 (MQ=255)
                                                                                                                                                  tACGAATACGGCGGATCGTCACGTCTTCAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGAc                                                                           <  1:307759/250‑1 (MQ=255)
                                                                                                                                                   aCGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGCTCaac                                                                                                                                                  >  1:166100/1‑178 (MQ=255)
                                                                                                                                                   aCGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCaac                                                                                                                                                  <  2:166100/178‑1 (MQ=255)
                                                                                                                                                     gAATACGGCGGATCGTCACGTCTCCAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTGCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTac                                                                                           <  2:260076/231‑1 (MQ=255)
                                                                                                                                                     gAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTac                                                                                           >  1:260076/1‑231 (MQ=255)
                                                                                                                                                         aCGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAAACCTTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGcctaccggtcacc                                                                   >  2:100964/1‑239 (MQ=255)
                                                                                                                                                         aCGGCGGATCGGCCTGTCTACAGGACGGGCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCCTCAGCAGTTCAGCACGAGACTGAATGTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTGCAACGCTACGACGTTCTTCGCTGACAGTACCCa                                                                   <  2:307787/251‑1 (MQ=255)
                                                                                                                                                          cggcggATCGTAACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAg                                                                  <  2:454547/251‑1 (MQ=255)
                                                                                                                                                              ggATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTc                                                              <  1:202259/251‑1 (MQ=255)
                                                                                                                                                                  cGTCATGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGTTCATTTTCTTCAGCATTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATg                                                           <  2:34198/250‑1 (MQ=255)
                                                                                                                                                                         tCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGaa                                                                                                               >  2:174989/1‑191 (MQ=255)
                                                                                                                                                                         tCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGaa                                                                                                               <  1:174989/191‑1 (MQ=255)
                                                                                                                                                                               gTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCACCGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAgc                                               >  1:397247/1‑248 (MQ=255)
                                                                                                                                                                                  cGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGtt                                          >  2:222959/1‑251 (MQ=255)
                                                                                                                                                                                                cAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATc                             <  2:45983/250‑1 (MQ=255)
                                                                                                                                                                                                   cTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAt                                                                                                                                                                                                         <  1:135118/75‑1 (MQ=255)
                                                                                                                                                                                                   cTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAt                                                                                                                                                                                                         >  2:135118/1‑75 (MQ=255)
                                                                                                                                                                                                        ggcaggTCATTTTATTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGAACTGCTCGCCATCAGGGCCGATCAACAAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTAAACGCTACGAAGTTCTTCGGTGACCGTACCCAGATTAATGGTCCGGGACAGCAGGTTGCGTGCACGAATCGTCCGTTc                     >  1:355555/5‑250 (MQ=255)
                                                                                                                                                                                                               cATTTTCTTCAACCGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCTCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCCACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGCTTCATGGTACGGGACAGACGGTTGCGAGCACGAATTGTCAATTCAcgcaggg             >  2:71947/1‑247 (MQ=255)
                                                                                                                                                                                                                aTTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGtt            <  1:222959/251‑1 (MQ=255)
                                                                                                                                                                                                                        cagcagTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCATGCGGGTTGAACGGtt    <  2:77227/251‑1 (MQ=255)
                                                                                                                                                                                                                        cagcagTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGtt    >  1:470228/1‑251 (MQ=255)
                                                                                                                                                                                                                             gTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTg  <  2:42445/248‑1 (MQ=255)
                                                                                                                                                                                                                             gTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTg  >  1:42445/1‑248 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
TTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTG  >  NC_000913/4639533‑4640001

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