Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,093,770 C→T G168D (GGC→GAC)  rhaD ← rhamnulose‑1‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,093,7700CT98.2% 198.9 / NA 55G168D (GGC→GAC) rhaDrhamnulose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (29/25);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACTTTTACGATGCCTAAGTTAGCCGC  >  NC_000913/4093532‑4094014
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
cATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                                                                                                                                                                                                           <  1:263146/250‑1 (MQ=255)
                tAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgc                                                                                                                                                                                                                          >  2:303967/1‑251 (MQ=255)
                 aaaCCTTCACTACTACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGACGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgcg                                                                                                                                                                                                                         >  2:208402/1‑251 (MQ=255)
                                  tgtgCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTAc                                                                                                                                                                                                        >  1:176804/1‑251 (MQ=255)
                                              tCTGCGGTGTCGATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGc                                                                                                                                                                                            >  2:141121/1‑251 (MQ=255)
                                                 gCGGTGTCGATTAAACCGAAGGTTTCATCCAGCGTCGTTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCg                                                                                                                                                                                           <  2:369480/249‑1 (MQ=255)
                                                   ggTGTCGATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCa                                                                                                                                                                                       <  2:405425/251‑1 (MQ=255)
                                                      gtCGATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGACGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt                                                                                                                                                                                    >  1:372067/1‑251 (MQ=255)
                                                        cGATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGCTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGCTCCGGTtg                                                                                                                                                                                  >  2:118038/1‑251 (MQ=255)
                                                        cGATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGGCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtg                                                                                                                                                                                  >  2:42264/1‑251 (MQ=255)
                                                          aTTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCAGTGGAAAGGCCACAAAACCAGCGAATGTTTTTGCATATCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTACCGTGCACCATACGCGGCAAAATGCAAACGCCATCCGGGAATACCACCAGACACTCGGGGCTGCCTTCCCACAGTTGGAGCGTGAAGACCGCGGGGTCGTTTTAAAGTAACTCGGTGCGGGCGATCAGGAtggg                                                                                                                                                                                >  2:319297/1‑250 (MQ=255)
                                                            tAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGGGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggtg                                                                                                                                                                               <  1:193985/250‑1 (MQ=255)
                                                             aaaCCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggc                                                                                                                                                                             >  1:132283/1‑251 (MQ=255)
                                                                 cGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtgg                                                                                                                                                                         >  1:273123/1‑251 (MQ=255)
                                                                  gAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtggc                                                                                                                                                                        >  2:291362/1‑251 (MQ=255)
                                                                  gAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtggc                                                                                                                                                                        >  1:68229/1‑251 (MQ=255)
                                                                  gAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcg                                                                                                                                                                            <  1:139484/247‑1 (MQ=255)
                                                                  gAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcg                                                                                                                                                                            >  2:139484/1‑247 (MQ=255)
                                                                              cAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATc                                                                                                                                                             >  1:333446/1‑250 (MQ=255)
                                                                                 cGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGGTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCAccc                                                                                                                                                         <  1:362297/251‑1 (MQ=255)
                                                                                      gTCCGCTTCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTGGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTGTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCt                                                                                                                                                    <  2:142923/251‑1 (MQ=255)
                                                                                                     cgtcGTGGAAGGGCCACAACACCAGTGAATGTTTTGGCTTCTCTTGTGAGTTTGCCGGGCCGATGTCGTCTGTGACGGGCACCAGCCACGGCAAAATGCCAAGGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGGCGTTTTCGAGGACATAGGTGAGGGTGATCAGGTTGGTGGCGTGTCAGTGCATGATCGCCCGATCTTTGCCGTTGGTGGCt                                                                                                                                     <  2:260902/248‑1 (MQ=255)
                                                                                                        cgaggAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGGATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGGCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGAGGCGTGGCAGTGCCTGATCACCCGATCTTTGCCGGTGGTGGCttt                                                                                                                                   >  2:309838/4‑250 (MQ=255)
                                                                                                          tGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAAt                                                                                                                                >  2:400017/1‑251 (MQ=255)
                                                                                                               gggCCACAACTCCAGCGATTGTTTTTGCATCTCTTGTGCGGTAGACTGGCAGTTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCGTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGATGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCt                                                                                                                           <  2:226677/251‑1 (MQ=255)
                                                                                                                           cAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCAATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCCCAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGCTCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGaa                                                                                                                >  1:368184/1‑250 (MQ=255)
                                                                                                                               gAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGtt                                                                                                                                                                                   >  1:174164/1‑179 (MQ=255)
                                                                                                                               gAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGtt                                                                                                                                                                                   <  2:174164/179‑1 (MQ=255)
                                                                                                                                            tctcTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGCCCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGt                                                                                              >  1:248709/1‑251 (MQ=255)
                                                                                                                                                 tgtgCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAATTGAGCCGGAAGTTc                                                                                            <  1:327110/248‑1 (MQ=255)
                                                                                                                                                 tgtgCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCgg                                                                                          <  1:84191/250‑1 (MQ=255)
                                                                                                                                                 tgtgCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCg                                                                                           >  1:5494/1‑249 (MQ=255)
                                                                                                                                                 tgtgCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCg                                                                                           <  2:5494/249‑1 (MQ=255)
                                                                                                                                                  gtgCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGaa                                                                                        <  1:42264/251‑1 (MQ=255)
                                                                                                                                                     cGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAACTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTg                                                                                     >  2:173082/1‑251 (MQ=255)
                                                                                                                                                     cGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTg                                                                                     >  2:293008/1‑251 (MQ=255)
                                                                                                                                                          gCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGAc                                                                                <  1:303967/251‑1 (MQ=255)
                                                                                                                                                            cTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACgg                                                                              >  2:122807/1‑251 (MQ=255)
                                                                                                                                                                  gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAAt                                                                                                                                >  2:87706/1‑195 (MQ=255)
                                                                                                                                                                  gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAAt                                                                                                                                <  1:87706/195‑1 (MQ=255)
                                                                                                                                                                     ttCGGCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTg                                                                      <  2:370852/250‑1 (MQ=255)
                                                                                                                                                                      tcggccGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTgg                                                                     <  1:174659/246‑1 (MQ=255)
                                                                                                                                                                              gCCGGGCACCATCCACGGCACAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAAcccc                                                             <  1:372837/250‑1 (MQ=255)
                                                                                                                                                                                     aCCATCCGCGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGTGGACGGCTTCGTTGGTTAACCCCCGAAGAAt                                                     <  1:118038/251‑1 (MQ=255)
                                                                                                                                                                                                 aaaaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCcgc                                         <  1:141121/251‑1 (MQ=255)
                                                                                                                                                                                                         aaCGCCATCCGAGAATAACACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGAGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGCGGGCGATCAGGATGGTGGCGTGGCAGAGCATGATCACCCGATCTTAGACGTTGGAGGCTTTAATGCGCTAGCAGTGGGAAAGGCAGTGAGCCGGAAGTTCGGAAGCGGGGCCGGCTTGGTTGGTTAACCCCCAAAGACTGTGGTAGCCCGAGCCGTCGc                                 >  2:389650/1‑251 (MQ=255)
                                                                                                                                                                                                           cGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTg                               <  1:208803/251‑1 (MQ=255)
                                                                                                                                                                                                               aTCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGCCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGCAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTCGCCCGCGCCGTCGCTGTCg                          >  2:68986/1‑250 (MQ=255)
                                                                                                                                                                                                                    ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGGTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCGTCGTTGGTTAACCCCCAAAGAATGTGGGAGCCCGCGCCGTCGCTGTCGACgcat                      <  2:342299‑M1/251‑5 (MQ=255)
                                                                                                                                                                                                                    ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTGTTCAAGTACATATGTGAGGGCGATCAGGTTGGTGGCGGGGCAGTGCATGATCTCCCGATCTGTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCGTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCTCCGTCGCTGTCGACgcat                      <  2:169563‑M1/251‑5 (MQ=255)
                                                                                                                                                                                                                      aaTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgc                    >  2:85875‑M1/1‑245 (MQ=255)
                                                                                                                                                                                                                              caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctc            >  1:97704‑M1/1‑237 (MQ=255)
                                                                                                                                                                                                                               aGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACCTAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctct           >  2:274475‑M1/1‑236 (MQ=255)
                                                                                                                                                                                                                                 acacTCTGTGCGGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTTCATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAATTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctag         <  1:349523‑M1/251‑18 (MQ=255)
                                                                                                                                                                                                                                        gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTGTCAGGCACATAGGCGAGGGCGATCAGGTTGGCGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCGCGCAGTGGGAAAGGAAGAGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgac  >  2:257973‑M1/1‑227 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
CATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCATCCAGCGTCGGTCCGCTGCCGAAGACGCCGTGGAAGGGCCACAACACCAGCGAATGTTTTTGCATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACTTTTACGATGCCTAAGTTAGCCGC  >  NC_000913/4093532‑4094014

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