Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,467,674 T→G L17R (CTG→CGG)  mgtA → Mg(2(+)) importing P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,467,6740TG97.8% 130.3 / ‑5.5 45L17R (CTG→CGG) mgtAMg(2(+)) importing P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (24/20);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

ATATGGAACCTGATCCCACGCCTCTCCCTCGACGGAGATTAAAACTTTTCCGGTAAGCCCGTCTTTTCACGGCGTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATC  >  NC_000913/4467429‑4467916
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
atatGGACCCTGATCCCACGCCTCGCCCCCGACGGAGATTAAAACTTTTCCGGGAAGCCCGGCTTTTCACGGCGTTACCGGATGCGTAAGTCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTGCCTTCAGGTAAGGCTGCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGGGGGACTCCTTATTTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGtg                                                                                                                                                                                                                                                  <  2:32546/250‑2 (MQ=255)
           gATCCCACGCCTCTCCCTCGACGGAGATTAAAACTTTTCCGGTAAGCCCGTCTTTTCACGGCGTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATcc                                                                                                                                                                                                                                       <  2:161620/250‑1 (MQ=255)
                  cGCCTCTCCCTCGACGGAGATTAAAACTTTTCCGGTAAGCCCGTCGTGTCACGGCGTTACCGGATGCGTAAGGTCGTGACGTTTTAACGTCCCTGCTCAGCTTTAGTACCTTCAGGTAATGCTTCGCCACGCCGGAAGACATTTCTGTACTGTTTCAGAGAGTGCGGAGGGACTCCTTATGTGGAAAGATATTTTTACCCGTCTCATTCTCCATTTACCTTCCCGGCTGTTTCAGCGGGATCCTTTGCCTg                                                                                                                                                                                                                               <  2:284892/251‑1 (MQ=255)
                   gCCTCTCCCTCGACGGAGATTAAAACTTTTCCGGTAAGCCCGTCTTTTCACGGCGTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTgg                                                                                                                                                                                                                              >  2:160061/1‑251 (MQ=255)
                   gCCTCTCCCTCGACGGAGATTAAAACTTTTCCGGTAAGCCCGTCTTTTCACGGCGTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATGTACCTTCCCGGCGGGTTCATCGTGATCCATTGCCTgg                                                                                                                                                                                                                              <  2:100193/251‑1 (MQ=255)
                    ccTCTCCCTCGACGGAGATTAAAACTTTTCCGGTAAGCCCGTCTTTTCACGGCGTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGgc                                                                                                                                                                                                                             >  2:250099/1‑251 (MQ=255)
                                         aatcTTTTCCGGTAAGCCCGTCTTTTCACGGCGTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCGTCCCGTCGGGTTCATCGTGATCCAGTGCCTGGCGCGCAGTAGACAGTGAATACg                                                                                                                                                                                                        <  2:106195/248‑1 (MQ=255)
                                                       aGCCCGTCTTTTCACGGCGTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTc                                                                                                                                                                                          >  1:386680/1‑251 (MQ=255)
                                                                   cACGGCGTTACCGGATGCGTAAGGCCTTGACTGTTTAACGTCCCTTCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCCCGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTAGGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTgcgc                                                                                                                                                                               <  1:568086/250‑1 (MQ=255)
                                                                      ggCGTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTg                                                                                                                                                                           <  1:521848/251‑1 (MQ=255)
                                                                       gCGTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTGCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCAGTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCTTAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGc                                                                                                                                                                          <  1:337104/251‑1 (MQ=255)
                                                                              cGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACCTTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTAGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGaaa                                                                                                                                                                    >  1:540290/1‑250 (MQ=255)
                                                                                                aCGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCgggggagcccga                                                                                                                                                   >  1:1242/1‑239 (MQ=255)
                                                                                                    tttAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCCCGACTGAAGCCATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGCATACGGTGGTCCCGCCGTCCTTAAGGGCGCATTGCCTGAAAATGGCGGTGGTGCCCggcaa                                                                                                                                                >  2:533805/1‑244 (MQ=255)
                                                                                                     ttAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGATTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGCGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGGGATGCCCgaggaaga                                                                                                                                             >  2:210700/1‑250 (MQ=255)
                                                                                                            cccTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTgg                                                                                                                                     <  1:288243/250‑1 (MQ=255)
                                                                                                             ccTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGaa                                                                                                                                   <  1:409864/251‑1 (MQ=255)
                                                                                                               tGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGAAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGGAGAAGAATTGTGGaaa                                                                                                                                  >  2:565833/1‑250 (MQ=255)
                                                                                                                        ttattaCCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAAACAGTGCGGCGGGACTCCTTATGTTTCAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGGTTCCCgaaactgacttggggaaaaagctcgac                                                                                                                          >  2:89232/1‑226 (MQ=255)
                                                                                                                          attaCCTTCAGGTAAGGCTTCGCCCCGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGacac                                                                                                                       <  1:458064/250‑1 (MQ=255)
                                                                                                                                  cAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCgg                                                                                                               <  1:233742/250‑1 (MQ=255)
                                                                                                                                       aaGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAggg                                                                                                          <  2:491739/250‑1 (MQ=255)
                                                                                                                                       aaGGCTTCGCCACGACTGAAGACATGTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAggg                                                                                                          <  1:200306/250‑1 (MQ=255)
                                                                                                                                            ttcgcaccgcattaggACATGTCTGTACTGTTTCAGACAGTGTGGAAGCACTCCTTATGTTTAAAGAACTTTTTACCCGGTTCATTAGCCATTTACCTTCCCGTCGGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGCGGTCCCGCCGTCGTTATGAGCGCATTGCCTGAAAATGGCGGAGTAGGCCGGAGAAGAATTGTTGGAAAAGTTCGACACCCCTCCTGAAGGGTTTAAt                                                                                                    >  2:257140/16‑251 (MQ=255)
                                                                                                                                                ccACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGGCCGAAGAAGAATTGTGGAAAAAGTTCGACACCCCTCCGGGAGGGTTAAATCAgg                                                                                                >  2:459371/1‑251 (MQ=255)
                                                                                                                                                ccACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAgg                                                                                                >  1:472650/1‑251 (MQ=255)
                                                                                                                                                     ccTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAACTTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTCAGTGCGCATTGCCTGAAAATGGCGGTGGTGCCCGAAGAAGAATTGTGGGAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCgg                                                                                             >  2:277483/1‑249 (MQ=255)
                                                                                                                                                     ccTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGGAGGGTTAAATCAGGCGGa                                                                                            >  1:113972/1‑250 (MQ=255)
                                                                                                                                                                    gTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGGTAAATCAGGCGGAAGTGGAATCTGCCcgg                                                                            >  1:388886/1‑250 (MQ=255)
                                                                                                                                                                             cagccAGTGCGGAGCGACTTCTTATGTCTAAAGAAGTTTTTCCCCGGCTCTTGCGCCATTGACCTTCCAGTCGGGCTCATCGTGATCCACTGTCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTACTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCACGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGAGTTAAATCAGGCGGAAGGGGAATCTGCCCGCTAACAACAt                                                                   <  2:132685/247‑1 (MQ=255)
                                                                                                                                                                             cAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGaa                                                                                                                                                                     <  1:423409/154‑1 (MQ=255)
                                                                                                                                                                             cAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGaa                                                                                                                                                                     >  2:423409/1‑154 (MQ=255)
                                                                                                                                                                                                            aGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCACGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGCAGAAAACTTGTGGAAAACGTTCGACACCCCTCCGGACGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGCACAAAATGGTGAAAATAAATTACCCGCACAACCACCgc                                   >  2:324387/1‑250 (MQ=255)
                                                                                                                                                                                                                aTTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCCTCCGGAAGGGGTAAATCAGGGGGAAGGGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCacaaacaccctagcc                               >  2:422320/1‑240 (MQ=255)
                                                                                                                                                                                                                aTTTTTACCAGGCTCATTCGCCATTTACCTTCTCGTCGGGTTCCTCGTGATCCCTTGCCTGGCGCGCAGCAGACAGTCAATACGGCGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGTTGCCCGCGGAAGAATTGTGGAAAACGTTCGACAACCATCCGGCAGGGGTCAATAAGGCGGAAGGGGAATCTGCCCGCGAACACCAAGGTGAAAAAAAATTTccccccacaaaaccatcgcc                               >  2:331528/1‑234 (MQ=255)
                                                                                                                                                                                                                      ccccGGCTCATTCGCCATTTACCTTCCCGGCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGGTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCgtggtg                         <  2:443753/250‑1 (MQ=255)
                                                                                                                                                                                                                        ccGGCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGGCACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTggg                       <  1:256395/251‑1 (MQ=255)
                                                                                                                                                                                                                         cGGCTCATTCGCCATTTACCTTCCCGTCGGGGTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGt                      <  1:333446/251‑1 (MQ=255)
                                                                                                                                                                                                                          ggCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGGCCGGCGGAGCATTGTGGAAAACCTTCGACCCCCATCCGGGCGGGATAAAACAGGCGGACGTGGGATCTGCCCGCCAACAACAAGGTGAGAATCAATTCCCCGCACAACAAccccagccccggtgggt                      >  1:274565/1‑235 (MQ=255)
                                                                                                                                                                                                                          ggCTCATTCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGGAGAAGAATTGTGGAAAACGTTCGAAACCCATCCGGAAGGGTTAAATCAGGCGGGAGTGGGATATGCCCGCGAACAACATGGTGAAAATAAATTACCCGCCCAACCACCGTCGCCGCGGTgggg                      >  1:298950/1‑249 (MQ=255)
                                                                                                                                                                                                                                 tCGCCATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGGGGTGCCCGAAGAAGAATTGTGGAAAACGTTCGAAACCCATCCGGAAGGGGTAAAACAGGCGGAAGGGGAATCTGACCGCGCACAACATGGGGAAAATAAATTACCCGCACAACAACCGTCGCCGGGGGGGGGACATTTAt              >  2:110015/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cGCCATTTTCCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATg             <  1:238451/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ccATTTACCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGCTTCCCCCAGAAGAAATGTGGAAACCGTTCACCACCCCTCCGGAAGGGTTAAATCAGGCGGGAGTGGCATCTGCCCGCGCACAACAAGGTGACAATAAATTACCCGCACAAAAACCACCGCCGGGGCGGGCACATTTATggg           >  1:278951/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tgaCCTTCCCGTCGGGTTCATCGTGATCCATTGCCTGGCGCTCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTTAAAATGTCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGTGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCt        <  2:361672/248‑1 (MQ=255)
                                                                                                                                                                                                                                             ttCCCGTCGAGTTCATCGTGATCAATTGCCTGGCGCGCAGCAGACAGAGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCACGAAGAAGGATTGTGGAAAACGTTCGACAACCATCCGGAAGGGGTAAAACAGGCGGAAGGGGAATCCGAACGCGGAAAAAATGGGGAAAATAACTTAACCGGCAAACAAACGGCGCCGGGGTGGGTGAATTTTTGGGGCTGTTATc  >  2:464336/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
ATATGGAACCTGATCCCACGCCTCTCCCTCGACGGAGATTAAAACTTTTCCGGTAAGCCCGTCTTTTCACGGCGTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATC  >  NC_000913/4467429‑4467916

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