Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,206,172 T→C 6.9% T29A (ACG→GCG)  ttdR ← transcriptional activator of ttdABT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,1720TC 6.9% 275.8 / 8.0 87T29A (ACG→GCG) ttdRtranscriptional activator of ttdABT
Reads supporting (aligned to +/- strand):  ref base T (42/39);  new base C (2/4);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

GTGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCT  >  NC_000913/3206024‑3206315
                                                                                                                                                    |                                                                                                                                               
gTGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGt                                                                                                                                                 <  2:1074859/149‑1 (MQ=255)
 tGAGGATTTCCAACGCATGTTCGTAAGAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTc                                                                                                                                                >  4:355408/1‑149 (MQ=255)
    ggATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTcgcc                                                                                                                                             >  1:1052056/1‑149 (MQ=255)
      aTTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTcgccgc                                                                                                                                           <  2:480378/149‑1 (MQ=255)
         tCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTgc                                                                                                                                        >  1:1130321/1‑149 (MQ=255)
           cAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTgcgg                                                                                                                                      >  3:460055/1‑149 (MQ=255)
           cAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTgcgg                                                                                                                                      >  3:295966/1‑149 (MQ=255)
            aaCGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCCgcggc                                                                                                                                     <  2:852912/149‑1 (MQ=255)
             aCGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTgcggcg                                                                                                                                    >  4:428058/1‑149 (MQ=255)
                cATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCtt                                                                                                                                 <  4:382992/149‑2 (MQ=255)
                  tGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttgg                                                                                                                               <  2:709797/149‑4 (MQ=255)
                   gTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTaaaa                                                                                                                              <  2:1130321/149‑1 (MQ=255)
                        tAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTccc                                                                                                                         >  4:1044073/1‑149 (MQ=255)
                         aacaacGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCg                                                                                                                        >  1:866290/1‑149 (MQ=255)
                         aacaacGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCg                                                                                                                        >  3:136891/1‑149 (MQ=255)
                           caacGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCt                                                                                                                      >  1:640483/1‑149 (MQ=255)
                            aacGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCtg                                                                                                                     >  2:375459/1‑149 (MQ=255)
                                ttGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGaa                                                                                                                 <  2:1052056/149‑1 (MQ=255)
                                 tGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAAc                                                                                                                <  3:407247/149‑1 (MQ=255)
                                 tGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAAc                                                                                                                <  1:1069374/149‑1 (MQ=255)
                                  ggCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCt                                                                                                                                  >  4:861719/1‑130 (MQ=255)
                                  ggCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCt                                                                                                                                  <  3:861719/130‑1 (MQ=255)
                                  ggCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCtt                                                                                                                                 >  2:896117/1‑130 (MQ=255)
                                  ggCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCtt                                                                                                                                 <  1:896117/131‑2 (MQ=255)
                                   gCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGa                                                                                                              <  2:931879/149‑1 (MQ=255)
                                   gCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGa                                                                                                              <  4:486289/149‑1 (MQ=255)
                                   gCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGa                                                                                                              <  4:264599/149‑1 (MQ=255)
                                       cTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTc                                                                                                          >  1:65843/1‑149 (MQ=255)
                                        tttCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAtt                                                                                                            >  4:287568/1‑146 (MQ=255)
                                        tttCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAtt                                                                                                            <  3:287568/146‑1 (MQ=255)
                                          tcagtcagCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTAc                                                                                                       <  2:556494/149‑1 (MQ=255)
                                                agCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTAc                                                                                                       >  4:1055073/1‑143 (MQ=255)
                                                agCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTAc                                                                                                       <  3:1055073/143‑1 (MQ=255)
                                                     cACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttg                                                                                                                                >  1:893189/1‑111 (MQ=255)
                                                     cACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttg                                                                                                                                <  2:893189/113‑3 (MQ=255)
                                                            cgAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCttt                                                                                     <  4:254530/149‑1 (MQ=255)
                                                             gagcggagcgGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTg                                                                                    <  2:318265/149‑1 (MQ=255)
                                                             gagcggagcgGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTg                                                                                    <  3:1044073/149‑1 (MQ=255)
                                                             gagcggagcgGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCCCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTTCCAGTGCCTGCAAGTCTTTg                                                                                    <  3:355408/149‑1 (MQ=255)
                                                               gcggagcggTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGc                                                                                  >  2:1101919/1‑149 (MQ=255)
                                                                    gcggTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAg                                                                             <  3:376088/149‑1 (MQ=255)
                                                                        ttCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCa                                                                         >  4:1118686/1‑149 (MQ=255)
                                                                        ttCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCa                                                                         >  4:1112763/1‑149 (MQ=255)
                                                                        ttCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCa                                                                         <  4:460055/149‑1 (MQ=255)
                                                                         tCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAg                                                                        <  4:131378/149‑1 (MQ=255)
                                                                          cagcagCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGc                                                                       >  1:778246/1‑149 (MQ=255)
                                                                            gcagcGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTa                                                                     >  3:179521/1‑149 (MQ=255)
                                                                             cagcGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCttt                                                                                     <  4:680283/132‑1 (MQ=255)
                                                                             cagcGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCttt                                                                                     <  3:22810/132‑1 (MQ=255)
                                                                             cagcGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCttt                                                                                     >  4:22810/1‑132 (MQ=255)
                                                                             cagcGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCttt                                                                                     >  3:680283/1‑132 (MQ=255)
                                                                             cagcGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTAt                                                                    <  2:866290/149‑1 (MQ=255)
                                                                                cGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCa                                                                 <  4:295966/149‑1 (MQ=255)
                                                                                        ccAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGGGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTATAGGCCAGCTATTCAGCATCGtt                                                         >  1:984911/1‑149 (MQ=255)
                                                                                         cAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCttt                                                                                     <  2:209725/120‑1 (MQ=255)
                                                                                         cAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCttt                                                                                     >  1:209725/1‑120 (MQ=255)
                                                                                         cAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTa                                                        >  1:290733/1‑149 (MQ=255)
                                                                                             gTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAAt                                                      <  1:221218/147‑1 (MQ=255)
                                                                                             gTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAAt                                                      >  2:221218/1‑147 (MQ=255)
                                                                                              tATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCt                                                                                               >  2:315331/1‑105 (MQ=255)
                                                                                              tATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCt                                                                                               <  1:315331/105‑1 (MQ=255)
                                                                                              tATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTAt                                                   <  4:917778/149‑1 (MQ=255)
                                                                                                    cAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAg                                             >  1:587917/1‑149 (MQ=255)
                                                                                                     aaGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGt                                            >  2:292264/1‑149 (MQ=255)
                                                                                                      aGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGtt                                           >  2:236303/1‑149 (MQ=255)
                                                                                                        gATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTtgt                                         <  3:150985/149‑1 (MQ=255)
                                                                                                          tttGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGa                                       <  2:640483/149‑1 (MQ=255)
                                                                                                           ttGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGAt                                      <  1:375459/149‑1 (MQ=255)
                                                                                                             ggATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATcc                                                 >  4:532/1‑136 (MQ=255)
                                                                                                             ggATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATcc                                                 <  3:532/136‑1 (MQ=255)
                                                                                                             ggATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATaa                                    <  3:1112763/149‑1 (MQ=255)
                                                                                                                   gCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATcc                                                 >  4:256622/1‑130 (MQ=255)
                                                                                                                   gCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATcc                                                 <  3:256622/130‑1 (MQ=255)
                                                                                                                          gACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTAttt                       <  4:136891/149‑1 (MQ=255)
                                                                                                                           aCAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTc                      <  2:65843/149‑1 (MQ=255)
                                                                                                                             aaaTGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGc                    >  1:7148/1‑149 (MQ=255)
                                                                                                                               aTGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCaa                  >  2:1121625/1‑149 (MQ=255)
                                                                                                                                   aGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAg                              >  4:863992/1‑133 (MQ=255)
                                                                                                                                   aGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAg                              <  3:863992/133‑1 (MQ=255)
                                                                                                                                          tGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAAtt                                                     >  4:334338/1‑103 (MQ=255)
                                                                                                                                          tGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAAtt                                                     <  3:334338/103‑1 (MQ=255)
                                                                                                                                          tGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTACCGCTGTGACCGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGGTGTGATAAGCGCAGTGTATTTCGCAACACATTGCCAt       >  2:588548/1‑149 (MQ=255)
                                                                                                                                             ccAAGCGTCGCCGCTGCGGCGCTAAAACTTCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGa                                       >  3:934045/1‑114 (MQ=255)
                                                                                                                                             ccAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGa                                       <  4:934045/114‑1 (MQ=255)
                                                                                                                                             ccAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTcc    >  4:79200/1‑149 (MQ=255)
                                                                                                                                             ccAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTcc    >  3:5923/1‑149 (MQ=255)
                                                                                                                                             ccAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTcc    >  3:991389/1‑149 (MQ=255)
                                                                                                                                               aaGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCt  >  3:287428/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GTGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCT  >  NC_000913/3206024‑3206315

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