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.9% 139.4 / ‑5.6 48S372S (AGT→AGCradADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (19/28);  minor base A (1/0);  total (20/28)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

GCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACG  >  NC_000913/4626783‑4627268
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
gagATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGCGATTCAGGCGCTGGTCGATCACTCGATGATGGCGACTCCGCGCCGCGTGGCAGTGGGGCTGGCACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGAc                                                                                                                                                                                                                                             >  2:128559/3‑251 (MQ=255)
 agATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACt                                                                                                                                                                                                                                            <  1:82887/250‑1 (MQ=255)
   aTGAAGTGACCGCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGGTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTa                                                                                                                                                                                                                                          <  2:17141/251‑1 (MQ=255)
          gACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGGCCGTTGCGGGTGGAGATGCAGGCGATGGGCTATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTATCGTTa                                                                                                                                                                                                                                    <  1:49111/250‑1 (MQ=255)
               ccggcagcgccgTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTActgctg                                                                                                                                                                                                                              <  2:211858/240‑1 (MQ=255)
                            gATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGTATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGTAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCa                                                                                                                                                                                                                 <  1:160873/251‑1 (MQ=255)
                             aTGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAg                                                                                                                                                                                                                <  1:27971/251‑1 (MQ=255)
                                gtggtgTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCt                                                                                                                                                                                                             <  1:98657/251‑1 (MQ=255)
                                     gtGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAGCAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTgcgtg                                                                                                                                                                                                        <  1:47408/250‑3 (MQ=255)
                                                cgcggccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGATCGCTg                                                                                                                                                                                             <  2:13818/242‑1 (MQ=255)
                                                   gtccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCCAGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGccc                                                                                                                                                                                          >  2:234857/7‑250 (MQ=255)
                                                    tccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGCTGGCGCATCAGCGAAGCGTGGCAGTGGCGCTGGCACAAAACCGTCTGGCACTCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCCAATGGCCGAACAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGAAACCGGAACCAGCGCCGACTTAGCGGTAATGCTGGCGATGGTTTCCAGCCTGCGTGACCGGCCGCTGccac                                                                                                                                                                                         >  2:242046/6‑249 (MQ=255)
                                                     cagttgctgGTGCAGATTCAGGCGAGCGTCGAGCACTCGATGATTGCGAATCCGCGCCACGTGGCAGTTGGGCTGGAACAACACCGTCGTGCAAGCCGGCTGGCTGTCGTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTGGGTGGCGGGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTCCTGGCGATGGTTTCCAGCCTGCTTGACAGACCGCTgactc                                                                                                                                                                                         <  2:16063/246‑5 (MQ=255)
                                                       gttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTATTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAgg                                                                                                                                                                                      >  1:152773/3‑251 (MQ=255)
                                                           ctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGAGGTCTGCAAATGGCCGATCAGGATGGGTTTGTGACCGTGTTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGAGTTACTGCTGGCGGTGGTTTCAGGCCTGCGGGACAGACCGCTGCCCCAGGGTCt                                                                                                                                                                                  >  2:42788/1‑251 (MQ=255)
                                                                 gagaTTCAGGCGCTGGTCGATCACTCGATGATGGCGAAGCCGCGCCGCGTGGCAGTGGGGCGGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTggtggt                                                                                                                                                                            <  2:200690/251‑1 (MQ=255)
                                                                               ggTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTc                                                                                                                                                               <  1:259074/250‑1 (MQ=255)
                                                                                                          gctgcgcgTGGCTGTGTGGCTGGTACAAACCCTTCTGGCAATCGTGCTGTCTGTCGTGCGCCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTATCCGAATCCAGCGCCGACTTAGCGTTACTGCTGGCGATGTTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGACCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGAAGATCCGCCCGGTGc                                                                                                                                   <  2:150488/247‑1 (MQ=255)
                                                                                                             ccgcGTGGCGGTGGGGCTGGAACAAAACCGGCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCa                                                                                                                                <  2:64895/251‑1 (MQ=255)
                                                                                                                 gTGGCAGTGGGGCTGGAACCAAACCGTCCGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGCTCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGCAACCAGCGCCGCCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGGGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCCGCgg                                                                                                                            >  2:172839/1‑251 (MQ=255)
                                                                                                                    gCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGACGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTc                                                                                                                          >  1:117287/1‑250 (MQ=255)
                                                                                                                           ggCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCCCGGTGGTCTGCAAATGGCCGCTCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGCCCGCTGCCACAGGATCTGGGGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCAGGGGCCCGGAGGTCAGGaccg                                                                                                                   >  2:262005/1‑247 (MQ=255)
                                                                                                                                  aCAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttc                                                                                                             >  1:166032/1‑246 (MQ=255)
                                                                                                                                    aaaaCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttcaga                                                                                                          >  2:221967/1‑244 (MQ=255)
                                                                                                                                       aCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcg                                                                                                      >  1:10762/1‑251 (MQ=255)
                                                                                                                                         cGTCTGGCAATCCTGCTGGCTGGCTTGCACCGTCACTGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcggc                                                                                                    <  2:82455/251‑1 (MQ=255)
                                                                                                                                            cTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGa                                                                                                  >  1:115028/1‑250 (MQ=255)
                                                                                                                                                   tCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGt                                                                                          >  1:131679/1‑251 (MQ=255)
                                                                                                                                                       gctggctgTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTc                                                                                      <  1:218039/251‑1 (MQ=255)
                                                                                                                                                                 ttGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGAt                                                                            >  1:216232/1‑251 (MQ=255)
                                                                                                                                                                   gCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGAtt                                                                           <  1:96571/250‑1 (MQ=255)
                                                                                                                                                                                      cAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTGCCGGCGGCTAACGtacc                                                       <  1:221967/251‑4 (MQ=255)
                                                                                                                                                                                       aaaTGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCg                                                      <  2:217413/251‑1 (MQ=255)
                                                                                                                                                                                                  cAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGAGCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGGTTGTTCCGGCGGCTAACGTTCCGAAAAAAGCGCCg                                          <  1:76894/251‑1 (MQ=255)
                                                                                                                                                                                                   aGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCg                                          <  1:244410/251‑1 (MQ=255)
                                                                                                                                                                                                           tttGTAAAAGTGGTCGGCGGCGTGCAGGTAACCGCAAACAGAGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGACACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTACAGAAGCGGCGAAACACGGCTTTAGCCGGGAGATAGTTCCGGCGACTAACGTAACGAAAAAAgacacggaaggggt                                  >  2:186850/1‑238 (MQ=255)
                                                                                                                                                                                                                  ccGTGGTCGGCGGCGTGAAGGTAACCGAACCCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGAtt                           <  2:166032/250‑1 (MQ=255)
                                                                                                                                                                                                                     tGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGAttttt                        >  1:79245/1‑251 (MQ=255)
                                                                                                                                                                                                                      gggcggcggCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTGACTGCTGGCGATGGTTTCCAGCCTGCGGGACAGACCGCTGCCACAGGATCTGGTGGTGGTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGTCCGGGGGCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTGCGGCGTCTAACGGACCTAAAAAAGCGCCGGAAGGGATGCAGAttgtt                        <  2:243160/247‑4 (MQ=255)
                                                                                                                                                                                                                        gcggcggCTTGAATGTAACCGAATCCTGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTTACAGACCTCTGCCACAGGATCTGGTGGTGTGTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCATATTTTTGGc                     <  1:96942/250‑1 (MQ=255)
                                                                                                                                                                                                                          ggcggcGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGt                   >  1:57114/1‑251 (MQ=255)
                                                                                                                                                                                                                          ggcggcGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGt                   >  2:165583/1‑251 (MQ=255)
                                                                                                                                                                                                                                gTGAGGGTAACCGAATCCAGCGCCGACTTAGCGTTACTGCTGGCGAGGGTTTCCAGGCTGCGTGACGGACCGCTGCCACAGGAGCTGGTGGTGTTTGGTGAAGTCGTGTTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAGTTCAGAAGCGGCGAAACACGGTTTTCGCGGGGGGATTGTTCCGGCGGCGAACGTACCGAAAAAATCGCCGGAAGGGATGCAGAGTTGTGGCGGTaaaaa             <  2:34466/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTAGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTAAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACAGAAAAAAGCGCCGGAAGGGATGCAGAATTTTGGAGTTAAAAAActct        >  2:72463/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aaCCGAAACCCGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGACACAGGATCGGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGCTCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCg     <  2:152773/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aaCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGTTCAGGAACGAATTTCAGAATCGGCGAAACACGGTTTGCGCCGGGGGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCg     <  1:242046/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cGAAACCAGCGCCGGCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAATCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAcg  <  1:228019/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cGAAACAAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAACGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAcg  >  2:235801/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
GCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACG  >  NC_000913/4626783‑4627268

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