Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,871,153 T→G N51T (AAC→ACC)  truD ← tRNA pseudouridine(13) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,871,1530TG98.2% 185.3 / NA 54N51T (AAC→ACC) truDtRNA pseudouridine(13) synthase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (29/24);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACG  >  NC_000913/2870918‑2871397
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
gTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGAAGCAGGTTTTGAGGATTCTa                                                                                                                                                                                                                                        >  2:483156/1‑250 (MQ=255)
 ttACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGGGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCCTGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGTGATTGAAGCCGGTTTTGAGGATTCTacc                                                                                                                                                                                                                                      >  2:459496/1‑249 (MQ=255)
    ccTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCCTACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGCGCGCATAACCACTATTCCGTAACAGCATGTTTGTCTTTTTTCACAGCGCAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCATTGCATCCTCCACAATCCGGGTATTGCAGCAGGTTTTGatgtttttaacc                                                                                                                                                                                                                                     >  1:348970/1‑238 (MQ=255)
    ccTTTCAGCGAGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGCAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATCACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCAAGCGACGCTGACTTCACGGGCATGAATTTTCAGGACTTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAAcc                                                                                                                                                                                                                                     >  2:222415/1‑249 (MQ=255)
       ttCAGCGCGCCTAAACGCAGCTCGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATGTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAat                                                                                                                                                                                                                                >  1:434226/1‑251 (MQ=255)
         cAGCGCGCCTAAACGCAGCTTGCGCTTGTGCAGCGCATACTCAAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCCGGCATTTCCTTGCACGGCACGCGAGCGCCTAACCACTGTTACGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGCATTTCGCCAGTGCATCCGACACAAACCGGGAATAGCAGCCGGTTTTGAGGATTATAACCAGActat                                                                                                                                                                                                                              >  2:373531/1‑251 (MQ=255)
                          gttggcgtttgtgccgcccaacacccAGCCCCTGGCCTCCTTCCAGTTGAAAGGCGCTAAGCTCGGGTTTTTCTTGCCCCGGCTCGCGAGCCTATACCCCCTTTCCCGTAACAGAATTTTTGTTTTTTTCCCCAGCGAATCTGACTTCACGGGCATGAATTTTAAGGATTTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCa                                                                                                                                                                                                             <  1:311076/227‑1 (MQ=255)
                                    gtgCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGTATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGc                                                                                                                                                                                                   >  1:105411/1‑251 (MQ=255)
                                            cATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCAAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCt                                                                                                                                                                                           >  2:509541/1‑251 (MQ=255)
                                                 cccAGCCCCTGGCAGCTTTCCAGGTGAAAGGCGCCAAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCAAAACCCCGGTCCCTTAACATCATGTTTGTTTTTTTGCCCAGCGAAGCTGACTCCACGGGCATGATTTTTCAGGAATTTGGCCAGGGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTccac                                                                                                                                                                                      <  1:560313/250‑1 (MQ=255)
                                                        ccgggcccccTTCCAGTTAAAAGCCGCTCAGATGGGGCATTTCCTTCCCCGGCACGCGCGCGCATAACCACTTTTCCGTACCAGCTTGTTTGTCTTTTTGCCCAGCGAATCTGACTCCCCGGGCATAAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACcaca                                                                                                                                                                               <  1:225091/243‑1 (MQ=255)
                                                           ggCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACaa                                                                                                                                                                              >  1:563521/1‑249 (MQ=255)
                                                               gCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCtt                                                                                                                                                                        >  1:199926/1‑251 (MQ=255)
                                                                ccTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTc                                                                                                                                                                       >  1:154989/1‑251 (MQ=255)
                                                                    ccAGTTGAAAGGCGCTCAGATGGGGCTTTTCCTTCCCCGGCACGCGAGCGCATACCCACTTTCCCGTAACAGCATGTTTTTCTTTTTCCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCgg                                                                                                                                                                    <  1:509541/250‑1 (MQ=255)
                                                                        ttGAAGGGCGTTCAGATGGGTCATTTCCTGGCCCCTCCCGCGAGCTCATACCCTGTGTTCCGTAACATCTTTTTTGTCTTTTCGCCCAGCAAAGCGGATTGCGCGGACATTAATTTTCGGGATTTTCGCCAGTGCTTCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATTCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTg                                                                                                                                                                <  1:326487/250‑1 (MQ=255)
                                                                        tggaaaggcgtataGATCGGGCTTTTCCTTGCCCGTCACGCGAGCGCATAACCTCTGTTCCGAAACAGCATGTTTGTTTTTTTGCCCAGCGAAGCTGACTTCACGGGCATAAATTTTCAGGATTTTCGCCAGGGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCGCCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTg                                                                                                                                                                <  1:373649/237‑1 (MQ=255)
                                                                           aaaGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCAACCCAAAGTCTTCCGGATTGGCtt                                                                                                                                                            >  2:183959/1‑251 (MQ=255)
                                                                            aaGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCttt                                                                                                                                                           >  2:567410/1‑251 (MQ=255)
                                                                                         ccggcctttCCCTGCCCGGGCACGCGGGCGCATAACCCTTTTCCCGTAACAGCTTGTTTGTCTTTTGGCCCAGCGAAGCTGACTTCCCGGGCATGAATTTTCAGAAATTTCGCCAGTCCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGt                                                                                                                                               <  1:150929/245‑1 (MQ=255)
                                                                                          gggaaTTTCCTTGCCGGGCACGAGAGCGTATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGAGAAGCTGACTGCACGGGCATGAATTTTGAGGAATTGCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGc                                                                                                                                             <  1:212186/247‑1 (MQ=255)
                                                                                                                   gAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTa                                                                                                                       >  1:248505/1‑250 (MQ=255)
                                                                                                                        aaaaacctcgttcccGTAACAGCTTGTTTGTTTTTTTGCCCAGCGAAGCTGCCTTAAGGGGTATGAATTTTAAGGAATTCCGCAGGGGAATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTgag                                                                                                                 <  1:68827/238‑1 (MQ=255)
                                                                                                                           aaCCACTGTTCCGTAACAGAATGTTTGGCTTTTCGCCCCGCGAAGCTGAATTCACGGGCATAAGTTTTCAGGACTTTCGCCCGTGCATCCGCCACAAAACGGGTATTGCAGACGGCTTTGAGGATTATAACCAGAATATGCTCACCTTCACCATCAGGCTCAACGCCCAAATCTTCCACCACCCCAATGTCTTCCGGACTGACTTTCAGACGCCCGGTGCCTGGCGTTTTACCGTGGGGGCAAGTGAGAtt                                                                                                              >  2:280207/1‑251 (MQ=255)
                                                                                                                                  tttCCTTAACAGCATTTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTGAAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCaaa                                                                                                        <  1:483156/249‑1 (MQ=255)
                                                                                                                                  ggtcccgaaacatcgtttttGTTTTTTTGCCAAGCGAAGCGTACTCCCCGGGCATAAATTTTCAGGATTTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCaa                                                                                                         <  1:36760/233‑1 (MQ=255)
                                                                                                                                  gTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAAc                                                                                                       >  1:393361/1‑251 (MQ=255)
                                                                                                                                  gTTCCGTAACAGCAGGTTTGTTTTTTGGCCCAGAGAAGCTGACTTCACGGGCATGAATTTACAGAAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCTGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCaa                                                                                                         <  1:277854/249‑1 (MQ=255)
                                                                                                                                   taccGTAACAGAATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGAAATGCAGACGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGACTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGCATGGCTTTCAGCAGCCCGGGGCCTTACGGTTTACCGCGGAGGGAAGTGAGATTATCACACt                                                                                                      >  2:529741/3‑251 (MQ=255)
                                                                                                                                    tCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGACGGATTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTACACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCAGTGCCTTGCGGATTACCGTGGAGGTAAGTGAGATTATCAAACTc                                                                                                     >  1:339626/1‑251 (MQ=255)
                                                                                                                                        gTAACGGCATTTTTGTTTTTTGGCCCAGCGAGGTTGATTTCCCGGCCAGGAATTTTCAGGAATTGCGCCGGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAAt                                                                                                  <  1:87397/251‑1 (MQ=255)
                                                                                                                                            cAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCAGCCACAAAACGGGTATTGCACCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGACTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGCGATTATCAAACTCAATCAtt                                                                                              >  1:534040/1‑251 (MQ=255)
                                                                                                                                               cATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCAtttt                                                                                            >  1:113532/1‑250 (MQ=255)
                                                                                                                                                    ttGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGcc                                                                                      >  1:158978/1‑251 (MQ=255)
                                                                                                                                                         ttttcGCCCTGCGAAGCTGACTTCCTGGGCATGAATTTTAAGGATTTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTaa                                                                                  <  1:288285/250‑1 (MQ=255)
                                                                                                                                                             tGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGt                                                                             >  1:172520/1‑251 (MQ=255)
                                                                                                                                                               cccAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCt                                                                                     >  2:237533/1‑241 (MQ=255)
                                                                                                                                                               cccAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTCCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCt                                                                                     <  1:237533/241‑1 (MQ=255)
                                                                                                                                                                ccAGCGAAGCGGACTTCACGGGAATGAATTTTCAGAAATTGCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAgc                                                                          <  1:349172/251‑1 (MQ=255)
                                                                                                                                                                  aGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAgcgc                                                                        >  1:28283/1‑251 (MQ=255)
                                                                                                                                                                            aCTTCCCGGGCATGAATTTTCAGAATTTTCCCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTca                                                              <  1:299470/251‑1 (MQ=255)
                                                                                                                                                                               aaagggtcATGAATTTAAAGAAATTGCGCCAGTCCACCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAg                                                           <  2:128205/244‑1 (MQ=255)
                                                                                                                                                                                 ccGGGCATGAATTTTCAGGATTTCGCCCAGTCCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGc                                                         <  1:529150/250‑1 (MQ=255)
                                                                                                                                                                                         gAATTTTCAGGAATTGCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCt                                                  <  1:183959/250‑1 (MQ=255)
                                                                                                                                                                                                taggtaactctccgcggtatccctccACAACACGGGTATTGCAGCCGTTTTTGAGGAGTCTAACCAGAATAGTCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCTCCACAAAGTCTTCCGGATTGGCTTGCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTAGCAAACGCAATCAGTTTGTTGCGTTAATGAGTAGCGTCGCCGGTTCACAGGCATTCCCTTCCCCa                                            <  2:55947/226‑1 (MQ=255)
                                                                                                                                                                                                cagtaaTTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCAc                                           <  1:352646/246‑1 (MQ=255)
                                                                                                                                                                                                     aTTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATAAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCgg                                     >  2:176490/1‑251 (MQ=255)
                                                                                                                                                                                                           ccAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATc                               >  1:480578/1‑251 (MQ=255)
                                                                                                                                                                                                              gctccATCCCCCACAAATCCGGTATTGCAGCCGGTCGTGAAGATTCTAACCAGAATATGCTCGTCGTCCCCAGCAGGCTCAAAGTCCAAACCTTCCACCGCCACAGAGTCTACCGGATTGGCTTTCAGCAGCACGGTGCCTTTCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCCATCATGTTGTTGCCTTAATACGTATCGCAATCGCTTCACAGGCAATCCCTTCCCCACTTcctggaaacccc                             <  2:434226/247‑11 (MQ=255)
                                                                                                                                                                                                                 gCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAg                           <  1:400268/250‑1 (MQ=255)
                                                                                                                                                                                                                         cacaAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCgt                  <  1:15508/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gggTATTGCAGCAGGTTTTGAGGATTCTAACCAGACTATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATATTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGCTTATCACACTCAATCATTTTGTTTCCTTAATGAGTAGCGCCACCGCTTCACAGGCACTCCCTTCCCCACGTCCGGTAAATCCCCATTTTTCCGTACTAGTgg          >  2:339543/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCAc   >  1:499214/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCCAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCACACCCTTCCCCACGTCCGGACAATCCCGGCTTTTCCGTAGTAGAGGCTTTCACg  >  2:553452/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
GTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACG  >  NC_000913/2870918‑2871397

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