Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,627,027 T→C S372S (AGT→AGC radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,627,0270TC97.4% 110.5 / ‑5.9 40S372S (AGT→AGCradADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (16/21);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCG  >  NC_000913/4626788‑4627270
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
gAAGTGACCTACGGCAGCTCCGTGATGGTTGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAAGCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCGGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGc                                                                                                                                                                                                                                          <  2:168623/251‑1 (MQ=255)
   gTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGtt                                                                                                                                                                                                                                       >  1:129337/1‑251 (MQ=255)
      aCCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTAct                                                                                                                                                                                                                                    >  1:24155/1‑251 (MQ=255)
           cggcaGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGGCTTGCAACGTCCCGGAGGGCTGCAAATGGCCGATCAGGAGGTGTTTGGTAAAGTGGTCGGCGGCGTGTAGGTAACCACAACACAAGCCGGCTTAGCGATACTGCTgg                                                                                                                                                                                                                               >  2:10125/5‑251 (MQ=255)
              caGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGACGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCgg                                                                                                                                                                                                                            >  1:232332/2‑250 (MQ=255)
                                       ggaaCGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGAcc                                                                                                                                                                                                   <  1:128710/251‑1 (MQ=255)
                                           cgcgTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTg                                                                                                                                                                                               <  1:258667/251‑1 (MQ=255)
                                                    tgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGAt                                                                                                                                                                                      >  1:80088/1‑251 (MQ=255)
                                                     gctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATc                                                                                                                                                                                     <  1:259842/251‑1 (MQ=255)
                                                     gctgGTAGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGAGTGGCAGCGGGGCGGGACCACAAACGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTCGGCTGCAAATGACCGATCAGGATGTGTATGTGAACGTGGTCGACGGCGGCAAGGTAACCGAAACCAGCGCCGACTTAGCGTAACTACTAGCGATGGTTTCCAGACTGCGCGACAGACCgatggcacaggccc                                                                                                                                                                                     >  2:116860/1‑238 (MQ=255)
                                                         gtggAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATACGCGCCGCGGGGCAGTGGGGCTGGAACAAAACCGACTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAAAGGCCGATCAGGATGTGTTTGAGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGCTGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGAACtgg                                                                                                                                                                                  >  1:374905/1‑250 (MQ=255)
                                                            gagaTTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGAAGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCtggtggg                                                                                                                                                                              >  1:180558/1‑250 (MQ=255)
                                                                              gATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCt                                                                                                                                                            <  1:278325/251‑1 (MQ=255)
                                                                                   cTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCATGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGAGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGc                                                                                                                                                         >  1:251894/1‑249 (MQ=255)
                                                                                            ggCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATcc                                                                                                                                              <  1:10125/251‑1 (MQ=255)
                                                                                                         cgcgTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAg                                                                                                                                 <  1:314535/251‑1 (MQ=255)
                                                                                                            gTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCgg                                                                                                                              >  1:149417/1‑251 (MQ=255)
                                                                                                             tGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGACACCAGCGCCGACTTAGCGTTGCTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGGGGGGGTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCACAGCGGt                                                                                                                             >  1:63024/1‑251 (MQ=255)
                                                                                                                    ggggCTGGAACAAAACCGGCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGaa                                                                                                                       <  1:337782/250‑1 (MQ=255)
                                                                                                                        cTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATAAGGCTATGTTTGTGAACTAGGTCGGCGGCGTGAAGGTAACCGCAACCAGCGACGACTTAGCGATACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCGGCCACAGGATCTGGAGGTGTTTGGTGAATTCGGGCTGGCCGGGGAGATCCGCCAGGTGCCCAGCGGTCAGGCACgca                                                                                                                   >  2:371700/1‑248 (MQ=255)
                                                                                                                          ggAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtt                                                                                                                 <  1:77150/250‑2 (MQ=255)
                                                                                                                            aaCAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttca                                                                                                              <  1:362152/251‑5 (MQ=255)
                                                                                                                            aaCAAAACCGGCTGGCAATCCTGCTGGCGGTCTTGCACCGTCACGGTGGTCTGCAACTGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGGGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttca                                                                                                              <  1:382317/251‑5 (MQ=255)
                                                                                                                                aaaCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAg                                                                                                          <  1:265240/251‑1 (MQ=255)
                                                                                                                                      tCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGTTGCGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACACGCGCCGACTTCGCGTTACTGCTGGCGCTGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGTTCTGGCGGTGTTTGGTGAAGTCGGGCTGGCAGGGCAGATCCGCCCGGTGCCCCGTGGCCAGGCACGCATTTCAGAAgcggcg                                                                                                     >  1:54870/1‑250 (MQ=255)
                                                                                                                                               ccTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGtt                                                                                           >  1:191799/1‑251 (MQ=255)
                                                                                                                                                         gtcttGCACCGTCACGGTGGTCGGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGGGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGTCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCgg                                                                                   <  1:350029/246‑1 (MQ=255)
                                                                                                                                                                                gccatgggacggtcggattggttttGTGAACGTGGTTGGCGGATGAATGGTAACCAAATCGAGGCCCGGCTTGGCGTTAGTGCGGGCAATGGTTTCCGGCCTGCGGTCCAGTCCGCTTTCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCCGGGAGGGGAGATCCGCCCGGTGCCCAGCGGTGAGGAACGAATTTCAGAAGCGGCGAAACACCGTTTTCGCCGGGTGATTGTTCCGGCGGCTAACGtac                                                          <  2:3608/229‑3 (MQ=255)
                                                                                                                                                                                      ggctgATCATGATGTGTTGGTGACCGTGTTCGGCGGCTGGAGGGTACACGAAACCAGCGCCGTCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGCCAGACCGCTGCCCCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGACCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCTATTGTTCCGGCGGCTAACGTACCGaaa                                                     <  2:301318/246‑1 (MQ=255)
                                                                                                                                                                                      ggCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGaaaa                                                    <  1:301528/251‑1 (MQ=255)
                                                                                                                                                                                          gATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAgc                                                >  1:379366/1‑251 (MQ=255)
                                                                                                                                                                                            tCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAgcgc                                              >  1:188408/1‑251 (MQ=255)
                                                                                                                                                                                                    tgtTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCCGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCAACAGGATCTGGGGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAAATCCGCCCGGAGCCCAGCGGTCAGGCACGACTTTCAGAAGCGGCGAAACACGGTTTTCGCCGGACGATTGTACCGGCGGCTaaagaacagaaaaatgcgcaggacggg                                      >  1:307450/1‑226 (MQ=255)
                                                                                                                                                                                                      tttGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAACCCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAAc                                                                                                                      >  1:80629/1‑169 (MQ=255)
                                                                                                                                                                                                      tttGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCGGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCTGTGCCCAGCGGTCAGGAAc                                                                                                                      <  2:80629/169‑1 (MQ=255)
                                                                                                                                                                                                        tgtgAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATg                                   <  1:151472/250‑1 (MQ=255)
                                                                                                                                                                                                                         gcgTGAAGGTAACTGAACTCAGCGCCGCCTTTGCGTTACTGCTGGCGATGTTGTCCAGCCTGCGTGACAGCCCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTGCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGTATGCAGATTTTTGGCGTTaaa                 <  2:129337/251‑1 (MQ=255)
                                                                                                                                                                                                                          ggTGAAGGCAACCGAACACAGCGCCGCCTTTGCGTGTCTGCTTGCGATGGTTTCCAGCCTTCGTGACATACCGCTGACACAGGATCTGGGGGTGTTTGGGGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTaaaa                <  2:300843/250‑1 (MQ=255)
                                                                                                                                                                                                                                  gTAACCGAAACCAGCGCCGATTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACtctc         <  2:191799/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aaaCCAGCGCCGCCTTAGTGTTACTGTTGGCGATGGTGTCCAGCCTGCGTGACAGACCGCTGCCACAGGATTTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGCATTTCAGAAGCGGCGAAACTCGGTTTTCGCCGGGCGATTGTGCCGGCGGCTAACTTGCCTAAAAAAGTGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAcgcg  <  1:205088/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
GAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCG  >  NC_000913/4626788‑4627270

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