Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,433,459 C→T G189D (GGC→GAC)  qorB ← NAD(P)H:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,433,4590CT97.3% 119.1 / NA 37G189D (GGC→GAC) qorBNAD(P)H:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (26/10);  total (26/11)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGCCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTG  >  NC_000913/4433216‑4433705
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
gggtgccaaatcaatTTGCTAAGTGTTTTGCTTTCATCAACCAGACCGCCTTTCGATGCGCCAACGTCAGAACCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGt                                                                                                                                                                                                                                                  <  1:238388/242‑1 (MQ=255)
   tGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGAcc                                                                                                                                                                                                                                              >  1:275804/1‑251 (MQ=255)
        aatcagtTTGCTAAGAGTTTTGCTGTAATCAAACAGACCGCATTTCGATGCGCAAACGACAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTACGACGCTTTTCCGTGCCGCGGCGAACTCGGCTTCGCTCAGATTTTGATCGGCAACCTGTTTGCCGCTCTGTTTGGTCCGCTCTGCCGCTAACTGTGTCAACGTCCAGGCAATATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACAGGCtt                                                                                                                                                                                                                                         >  2:361296/1‑251 (MQ=255)
        aatcaatTTGCTCAGTGTTTTGCTGTCATCAAACAGACAGCCTTTCGATGAGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTACGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGCTTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGAAAACCTTGTCTTCGGGACCGGCtt                                                                                                                                                                                                                                         >  2:79577/1‑251 (MQ=255)
                  cTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGCTAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCATGCACTCTCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCAc                                                                                                                                                                                                                               >  1:5244/1‑251 (MQ=255)
                                 tcatcaAACAGACCGCCTTTCGATGCGCAATCGTCAGACTACCCAAACATATCAGCCACTCAGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGCACTCGGCTTAGCTCAGATTTTGATAGGTCACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGACAACGTCCAGGCCCTATCGCCCGCCAGTTGGAAAACCTTGTCTTCGCGACGGGCTTCGCTAATCACGCGTGCCGCAgtcgc                                                                                                                                                                                                                >  2:184780/1‑247 (MQ=255)
                                        aCAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATa                                                                                                                                                                                                          >  2:294000/1‑250 (MQ=255)
                                           gcccGCCTTTCGATGCCCCACCGTCAGAACCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCt                                                                                                                                                                                                      <  1:161397/249‑1 (MQ=255)
                                                            ccagcgcgaaaaTCCGCCAGCTTACCCGCCATTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCTCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAGCCGTGTCTTCGTGACCGGCGTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGc                                                                                                                                                                                     <  1:64359/241‑1 (MQ=255)
                                                                      aaTCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGcc                                                                                                                                                                            >  1:94506/1‑250 (MQ=255)
                                                                           gcatagattcacgcaggTCCGTCGGGTAGGCCGACTCGTTTCAGGGCCGCGGCGAAATCGGCTTCCCACAGATTTTGATAGTTACCCTGTTTGCCGCGCTGTTTGGACAGCTCTTCCGCTAACTGTGTCAGCGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCCATTTTGCTATcgcc                                                                                                                                                                      <  2:271424/235‑1 (MQ=255)
                                                                               gCATATCAGCCAGCCCGCCGGGCAGTCCGCCGCTTTTCAGTGCCGCGGCGAAATCGGCCTCGCTCAGATTTTGCTAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATAGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATcgccagcc                                                                                                                                                                  >  2:161762/1‑251 (MQ=255)
                                                                                       cccAGTCCGTCGGGCAGTCCGCCGCTTTTCGGTGCCGGGGCGAAATCGGCTCCGCTCAGATGTGGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACGGTGCCAACGCCCAGGCACTATCGCCCGCCAGTGCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGGGATTTTGCCATCGCCCGCCGCACCGAt                                                                                                                                                          <  2:199917/250‑1 (MQ=255)
                                                                                                     cggTCCGACGCTTTCCAGGCCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATGGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGtt                                                                                                                                            <  1:266982/249‑1 (MQ=255)
                                                                                                         ccGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTcc                                                                                                                                          >  1:234850/1‑249 (MQ=255)
                                                                                                                     attGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGccc                                                                                                                                                                                                  <  2:106675/179‑1 (MQ=255)
                                                                                                                     aGTGCCGCGGCGAACTCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGccc                                                                                                                                                                                                  >  1:106675/1‑181 (MQ=255)
                                                                                                                                gAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCCACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGGTAAATACGCCGTGCTCCAAGGCTGCCGGGGCGCCGGCGAggg                                                                                                                 >  2:68004/1‑250 (MQ=255)
                                                                                                                                 aaaTCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGt                                                                                                                 >  1:364713/1‑250 (MQ=255)
                                                                                                                                         ttCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGc                                                                                                        >  2:262509/1‑251 (MQ=255)
                                                                                                                                             cTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCAAGGCACTATCGCCCGCCAGTTCGTAAAACTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGACATCGCCCGCCGCACCGATAAATACGCCGTGTTCCCGTGCTGCCGGGGCGCTGGCGAGGAAGTTTTCGCTGt                                                                                                     >  1:371364/1‑250 (MQ=255)
                                                                                                                                                          aGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTgcgc                                                                                       >  1:162127/1‑251 (MQ=255)
                                                                                                                                                             tAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCgcagc                                                                                    >  2:120668/1‑251 (MQ=255)
                                                                                                                                                               aCCTGTTTGCCGCTCTGTTTGGTCCGCTCTGCCGCTAACTTTGTCACCGTCCAGGCACTATCCCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGGGCCGCAGCGGCCGCATAATCGTCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCTCTGTACCAGCCGTTGCGcagcag                                                                                  <  1:93631/251‑1 (MQ=255)
                                                                                                                                                                         cGCTCTGATTGGTCCGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACg                                                                        >  1:109098/1‑251 (MQ=255)
                                                                                                                                                                                 ttGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGcc                                                                   >  1:10544/1‑248 (MQ=255)
                                                                                                                                                                                          tctGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATcagcc                                                         >  1:306142/1‑249 (MQ=255)
                                                                                                                                                                                           ctGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCACCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGTCGATGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCCGCAGGGTGTAAACGATGCCAGACTcagccagc                                                      >  1:286256/1‑251 (MQ=255)
                                                                                                                                                                                                gCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGCGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAAc                                                                         >  2:5299/1‑227 (MQ=255)
                                                                                                                                                                                                gCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAAc                                                                         <  1:5299/227‑1 (MQ=255)
                                                                                                                                                                                                          tCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCg                                       >  1:33294/1‑251 (MQ=255)
                                                                                                                                                                                                          caaaCGTCCAGCCACTATCCCCCGCAATTTCGTAAACCATGCCTTCGTGACCGGCTCCGCTAATCCCGCGTGCCGCAGCTGCCGAAAAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAAACAGCCAGCATTTTCTCCGTCTCg                                       <  1:91609/249‑1 (MQ=255)
                                                                                                                                                                                                           cAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCg                                       >  1:17606/1‑250 (MQ=255)
                                                                                                                                                                                                                                cgccAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCgg                 >  1:45939/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gccAGTTCGTAAACCTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCg                  >  1:292439/1‑249 (MQ=255)
                                                                                                                                                                                                                                    aGTTCTTAACCCTTGTCTTCGTGACCGGCTTCGCTATCCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGGGGAGGGGTTTTTGCCATGGCCCGCCGCACCGATAAATCCGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCgg             <  2:364713/251‑1 (MQ=255)
                                                                                                                                                                                                                                               cTTGTCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTg  >  2:264539/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
GGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGATTTTGATAGGTAACCTGTTTGCCGCTCTGTTTGGTCAGCTCTGCCGCTAACTGTGTCAACGTCCAGGCACTATCGCCCGCCAGTTCGTAAACCTTGCCTTCGTGACCGGCTTCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTG  >  NC_000913/4433216‑4433705

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