Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,042 T→C M44T (ATG→ACG)  radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,0420TC90.0% 92.8 / 2.4 40M44T (ATG→ACG) radADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (19/17);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

CACCATCCGTCACGCTGACCTGATGGGGGTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGC  >  NC_000913/4625797‑4626275
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
cACCATCCGTCACGCGTACCTGATGGGGGTATTCTGCATCCTCTCAGGGAGCCTGAATCAGAAGTAAGTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAAcggtgg                                                                                                                                                                                                                                      <  2:61037/251‑1 (MQ=255)
     tCCGTCACGCTGACCTGATGGGGGTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGt                                                                                                                                                                                                                                 >  1:370163/1‑251 (MQ=255)
        gTCACGCTGACCTGATGGGGGTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAAc                                                                                                                                                                                                                              >  1:101988/1‑251 (MQ=255)
         caccGCTGACCGAATGGGGGTATTCTACATCCTCTCAGGCAGCCTAAATCATAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGGTTGTAATGAATGCGGTGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAAc                                                                                                                                                                                                                              <  2:257600/247‑1 (MQ=255)
          cACGCTGACCTGATGGGGGTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACg                                                                                                                                                                                                                             <  1:188084/250‑1 (MQ=255)
                                            cAGGCAGCCTGAACCAGAAGTGATGTCTCGCCCGTAAGCCTGCCGGCGGCACAGTATGAACGAGGTACACCGTGGAAAAAGCTCCCAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTTTCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCGTGCTGCGGCGCCAATGGTGGTGCGTAACGAGCGTCGCAGCGGCTATGCCGGTAGCGCCGGGGTg                                                                                                                                                                                          <  2:374119/251‑1 (MQ=255)
                                                  gCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCaaaa                                                                                                                                                                                    <  1:292565/251‑1 (MQ=255)
                                                      gAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTcc                                                                                                                                                                                <  1:223090/251‑1 (MQ=255)
                                                                         gcccgccATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATc                                                                                                                                                               >  1:218897/1‑249 (MQ=255)
                                                                           ccgccATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGAGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGGCCAGAAACTCTCCAATATCAGcc                                                                                                                                                           >  1:85979/1‑251 (MQ=255)
                                                                             gccATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCCGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAAATCTCCGATATCAGCCtt                                                                                                                                                         >  1:197336/1‑251 (MQ=255)
                                                                                  ccTGCGGGCAGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCCACGGTGGCGCGAAACCAACGTCTCCGCGGCTAAGCCGAGAGAGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGGTATCAGcactga                                                                                                                                                       >  2:156921/1‑243 (MQ=255)
                                                                                               caGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGtttt                                                                                                                                       <  1:323855/251‑1 (MQ=255)
                                                                                                  aaTTAAAGAGGTACACCCGTGCAAAAGCTCCTAAACGCGCCTTGTTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTTGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTTGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGTTGCCGCGTTTTTcc                                                                                                                                    <  2:167047/250‑1 (MQ=255)
                                                                                                    ttAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCAc                                                                                                                                  >  1:306084/1‑251 (MQ=255)
                                                                                                           gggactacGTGGCAAAAGCGCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGTTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCTTTTTTCCACCGGAttt                                                                                                                           <  2:218897/244‑1 (MQ=255)
                                                                                                           ggTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTATCGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGAttt                                                                                                                           >  1:171989/1‑251 (MQ=255)
                                                                                                                acCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTaa                                                                                                                         >  1:228873/1‑248 (MQ=255)
                                                                                                                           aaGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACcgc                                                                                                            >  1:103998/1‑251 (MQ=255)
                                                                                                                            aGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACcgcg                                                                                                           >  1:3352/1‑251 (MQ=255)
                                                                                                                                        cgcCTTTGTTTGTAATGAATGTGGGGCCGATTACCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCCCCGAGGTGCGTCTTGCTGCTTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcg                                                                                               <  2:300039/251‑1 (MQ=255)
                                                                                                                                             ttgtttgtAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGgcggcggcg                                                                                            <  1:156921/249‑1 (MQ=255)
                                                                                                                                                      atgaatgCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGa                                                                                 <  1:194751/251‑1 (MQ=255)
                                                                                                                                                                    gATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGAt                                                                   >  1:75724/1‑251 (MQ=255)
                                                                                                                                                                     aTTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGAtt                                                                  >  1:257419/1‑251 (MQ=255)
                                                                                                                                                                          ccGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCgg                                                             <  1:116630/251‑1 (MQ=255)
                                                                                                                                                                           cgcgCTAGCAGGGGCAGTGTAGTGCCTGTCATGCCAGGAACCCCAACTCCGAGGTGCGGCTTGCTGCGTCGCCAACGGAGGCGCGTAACGAGCGTCTCAGCGGCTATACCGGTAGCGCCGGGGTGGCAAAAGTCCAGAACCTCTCCGATCTCAGACTTGAAGAGCTGCCGCGTTTTTGCACCGGCTTTTACGAGCTCGACCGCGTACTAGGCGGCGGCGAGGTGCCAGGAAGGTCCCTTCTGATTTGCggg                                                            >  2:311356/1‑250 (MQ=255)
                                                                                                                                                                            gcgcTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGt                                                            >  1:220261/1‑250 (MQ=255)
                                                                                                                                                                                           gtgcagtggCTGTCATGCTCGGAACACCATCACCGAGGTGCGTCTTGCGGCGTCGCCAACGGTGGCGCGTAACGAGCGACTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTTCGGGGa                                            <  2:308967/251‑1 (MQ=255)
                                                                                                                                                                                           gtgcagtgcCTGTCAGTCCTGGAAAGCCATCACCGTGGTGCGTCGTGGTGCGTCGCCATCGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCGCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCGCCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGa                                            <  2:184658/251‑1 (MQ=255)
                                                                                                                                                                                             gcagtgcCTGTCATGCCTGGACCACCATCACCGAGGGGCGTCTTGCTGCGGCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTGTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGaa                                           <  2:101988/250‑1 (MQ=255)
                                                                                                                                                                                                      gTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAAAGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGGAACCCTGGGGCGGGGGAATCCAAgctg                                 >  2:78805/1‑251 (MQ=255)
                                                                                                                                                                                                      gTCATGACTGGAACACCATCACCGAGGGGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGCCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGCGCTGCCGCGGTTTTCCACCGGATTTAATGAGTTCGACCGCGTACTAGGCGGGGGCGTGGTGCCAGGAAGTGACATTCTGATTGGCGGCAAccactggacggggaatgtcaccgtg                                 >  2:6628/1‑228 (MQ=255)
                                                                                                                                                                                                        cATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctgtt                               <  1:69083/251‑3 (MQ=255)
                                                                                                                                                                                                        cATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctgt                                <  1:87702/250‑2 (MQ=255)
                                                                                                                                                                                                            ccTGGAACACCATCACCGAGGTGCGTCTTTCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctgttgctg                           <  1:25996/251‑7 (MQ=255)
                                                                                                                                                                                                            ccTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCt                                                      >  1:111305/1‑224 (MQ=255)
                                                                                                                                                                                                            ccTGGAACACCAGCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCATCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCt                                                      <  2:111305/224‑1 (MQ=255)
                                                                                                                                                                                                                        tCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGGGCGGGGGAATCCACGCTGTTGCTGCAAACGCTGTGc               >  2:66711/1‑251 (MQ=255)
                                                                                                                                                                                                                             gAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACt          >  1:61531/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCcagc  >  1:254462/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
CACCATCCGTCACGCTGACCTGATGGGGGTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGC  >  NC_000913/4625797‑4626275

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: