Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,408,129 A→G 6.4% K27R (AAA→AGA)  ydaN → putative Zn(II) transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,408,1290AG 6.4% 200.8 / 2.6 62K27R (AAA→AGA) ydaNputative Zn(II) transporter
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base G (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

TCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGG  >  NC_000913/1407869‑1408388
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                   
tCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATa                                                                                                                                                                                                                                                       >  1:189956/1‑275 (MQ=255)
     tAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGCTCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCTGTAAACCGCTGGAAAATACAGAt                                                                                                                                                                                                                                                  >  1:331995/1‑275 (MQ=255)
         gaaagaagtcaTTTAACACTTCACACTGTTTGTATACGGAAAAGCATTTTTCTTTTTTTATTCAATTTAGGCAGAATTTGATTAATCATTTCAGGGTAATGGGGTGATGAGAGGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCTATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGtga                                                                                                                                                                                                                                              <  1:169858/265‑1 (MQ=255)
          aTAGACGCTATTTAACAATGAACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGtgat                                                                                                                                                                                                                                             <  2:213104/275‑1 (MQ=255)
              aCGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAACTACAGAtgtgct                                                                                                                                                                                                                                             >  1:243243/1‑269 (MQ=255)
               cGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATtgatg                                                                                                                                                                                                                                        <  2:20623/275‑1 (MQ=255)
                    tttAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAgcg                                                                                                                                                                                                                                   <  1:225725/275‑1 (MQ=255)
                            ttCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTg                                                                                                                                                                                                                           >  1:120838/1‑275 (MQ=255)
                                        gTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAAGCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCt                                                                                                                                                                                                               <  2:440889/275‑1 (MQ=255)
                                               gAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGCTGAAGGGCATCCCTGTTGCCTCCACCTTAATTa                                                                                                                                                                                                         >  2:15372/1‑274 (MQ=255)
                                               gAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTa                                                                                                                                                                                                         >  1:157212/1‑274 (MQ=255)
                                               gAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGGTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCACGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAACATACAGATGCGATTGATGAAGGGCATACCTGTTGCCTCCACCTTGATTAt                                                                                                                                                                                                        >  1:511212/1‑275 (MQ=255)
                                                   aGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTa                                                                                                                                                                                                     <  2:157212/274‑1 (MQ=255)
                                                          tGCTTCTTGTATTTAATTTTGACAGGTTGTTAGTAGTGAGTGGAGGGTAATAGGGGGATGAGAGGTGCGGAAACAGGTCAAGAGGGGTAAACGACAAAATAATGCGTTGATGAGGGGGGTTCTGTGGAGGCGGTTAAGGGGTAGGACGTTAATGTCCAGGAGCCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAGCCGCTGGAAAATATAGATGTGGTTGATGATGCGCATCTCTGTTGGCTCCACCTTAATTATGTACCCCATGa                                                                                                                                                                                             <  2:351691/275‑1 (MQ=255)
                                                                          tttAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCATAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCCCCTTAATTATGTACACCATGATAGCGCCCAAtgg                                                                                                                                                                                >  1:376003/1‑272 (MQ=255)
                                                                                        tattaATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGc                                                                                                                                                                                                                                                                >  2:278214/1‑179 (MQ=255)
                                                                                           taATCATTTCAGGGTGATGGGGTGATGAGATGTTGAGTAACAAGGCCAGAAGGCTAGACTACAAAAGAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTGGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCGTCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCgctgct                                                                                                                                                            <  1:256075/275‑1 (MQ=255)
                                                                                           taATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCgctgct                                                                                                                                                            <  2:58033/275‑1 (MQ=255)
                                                                                            aaTCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGa                                                                                                                                                                                                                                                            <  1:278214/179‑1 (MQ=255)
                                                                                                       ggTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAacgt                                                                                                                                                <  2:189956/275‑1 (MQ=255)
                                                                                                            tGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTg                                                                                                                                           >  1:507085/1‑275 (MQ=255)
                                                                                                            tGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTg                                                                                                                                           >  1:239774/1‑275 (MQ=255)
                                                                                                            tGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGACTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTg                                                                                                                                                                                                                        >  1:134548/1‑198 (MQ=255)
                                                                                                                gTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTc                                                                                                                                                                                                                    <  2:134548/198‑1 (MQ=255)
                                                                                                                 tgatgaGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAAt                                                                                                                                                                                                                                                        >  1:433987/1‑161 (MQ=255)
                                                                                                                     gagaTGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAg                                                                                                                                                                                                                                                    <  2:433986/161‑1 (MQ=255)
                                                                                                                        aTGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAg                                                                                                                                                                                                                                     >  1:460688/1‑173 (MQ=255)
                                                                                                                        aTGGTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTggcg                                                                                                                                <  2:392230/274‑1 (MQ=255)
                                                                                                                          gTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTcc                                                                                                                                                                                                                   >  2:23338/1‑189 (MQ=255)
                                                                                                                            tGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCa                                                                                                                                                                                                                                 <  2:460687/173‑1 (MQ=255)
                                                                                                                              cGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATTATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCt                                                                                                                                                                                                               <  1:23338/189‑1 (MQ=255)
                                                                                                                                      ggCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCAcgcg                                                                                                                  >  2:445089/1‑274 (MQ=255)
                                                                                                                                      ggCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGACGGTAGAGGAGGCGTTAAACCGCTGGAAAATACAGATGTGATTGAGGAAGCGCATCCCTGTTGGCTCCACCTTAATTATTTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGt                                                                                                                 >  1:376957/1‑275 (MQ=255)
                                                                                                                                        ccAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGTTGGCGGGCGCGAGCACCCGt                                                                                                                 >  1:503568/1‑273 (MQ=255)
                                                                                                                                        ccAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTc                                                                                                                >  1:72137/1‑274 (MQ=255)
                                                                                                                                                cTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCa                                                                                                       >  1:394652/1‑275 (MQ=255)
                                                                                                                                                    aCTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAGACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGcc                                                                                                    >  1:247170/1‑274 (MQ=255)
                                                                                                                                                                 tGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCCCTGGCGGGCGAGAGAACGCGTCCCCGAGTAAGCCGGCCTGGTGGAgg                                                                                       >  2:317584/1‑275 (MQ=255)
                                                                                                                                                                   cGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGATGGCGCCAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGATCGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGc                                                                                      >  1:319684/1‑274 (MQ=255)
                                                                                                                                                                    gTTGATGGTGGGGGGACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGc                                                                                      <  2:511211/273‑1 (MQ=255)
                                                                                                                                                                         tggtgGAGGCACGGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGTATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGa                                                                               <  2:376002/272‑1 (MQ=255)
                                                                                                                                                                             ggAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTAc                                                                           >  1:55861/1‑275 (MQ=255)
                                                                                                                                                                                   aCTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAGACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTgcg                                                                     <  2:247170/275‑1 (MQ=255)
                                                                                                                                                                                   aCTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTgcg                                                                     <  2:486689/275‑1 (MQ=255)
                                                                                                                                                                                    cTGTGGAATCGATTAAGGGACAGGAGGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTgcg                                                                     <  1:317584/274‑1 (MQ=255)
                                                                                                                                                                                     tgtgGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTgcgc                                                                    >  1:428955/1‑274 (MQ=255)
                                                                                                                                                                                     tgtgGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCt                                                                   <  1:445090/275‑1 (MQ=255)
                                                                                                                                                                                            gCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGAGGGTAGAGGCGGCGTTAAACCGCTTGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGtata                                                              <  1:15372/273‑1 (MQ=255)
                                                                                                                                                                                               aTTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACgg                                                         >  2:56324/1‑275 (MQ=255)
                                                                                                                                                                                                          ggACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGGCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCCGTATAAACGGCAGCAcctaa                                               >  1:107557/1‑271 (MQ=255)
                                                                                                                                                                                                                      cccGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAGACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGAt                                   >  1:108864/1‑274 (MQ=255)
                                                                                                                                                                                                                                 tATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGcc                       <  1:56324/275‑1 (MQ=255)
                                                                                                                                                                                                                                  aTTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGcc                       <  2:72137/274‑1 (MQ=255)
                                                                                                                                                                                                                                   tttGCCTGGATGCTGGATGGTAGAGGCGGCGTTAGACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCAt                     <  2:108864/275‑1 (MQ=255)
                                                                                                                                                                                                                                   tgtGCCTGGATGCGTGGTGGTAGAGGCGGCGTTAACCCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCAt                     <  2:319684/273‑1 (MQ=255)
                                                                                                                                                                                                                                     tGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCAt                     >  2:135150/1‑273 (MQ=255)
                                                                                                                                                                                                                                     tGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCAt                     <  1:135150/273‑1 (MQ=255)
                                                                                                                                                                                                                                     tGCCGGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGAGGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGAGAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGc                   <  2:120838/275‑1 (MQ=255)
                                                                                                                                                                                                                                        cgggaggctggatgGTAGAGGCGGCGGTAAACCGCTGGAAAATACAGATGGGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCgt                 <  2:394651/272‑1 (MQ=255)
                                                                                                                                                                                                                                                tggtggggagagGCGGGGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGGGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGtatata          <  2:376956/265‑1 (MQ=255)
                                                                                                                                                                                                                                                      gTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGATGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGTTGAACGCCCCGATCAACTGGTCGCCATGCGCGTATATATGGACggg  >  1:86316/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                cgTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAacgta                                                                                                                                               >  2:51363/1‑124 (MQ=255)
                                                                                                                                                                                                                                                                    aaaCCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTg                                                                                                                                           <  1:51363/124‑1 (MQ=255)
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                   
TCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGTTAATGTCCCGGATGCAGTATTTGCCTGGATGCTGGATGGTAGAGGCGGCGTTAAACCGCTGGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGG  >  NC_000913/1407869‑1408388

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: