Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,206,176 C→A 11.9% A27A (GCG→GCT ttdR ← transcriptional activator of ttdABT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,1760CA11.9% 277.3 / 24.0 84A27A (GCG→GCTttdRtranscriptional activator of ttdABT
Reads supporting (aligned to +/- strand):  ref base C (35/39);  new base A (3/7);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.21e-01

CGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGG  >  NC_000913/3206038‑3206322
                                                                                                                                          |                                                                                                                                                  
cGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGgcgc                                                                                                                                          <  3:507058/149‑1 (MQ=255)
   aTGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTaa                                                                                                                                       >  4:122536/1‑149 (MQ=255)
   aTGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCtt                                                                                                                                        <  2:665850/148‑2 (MQ=255)
    tGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTaa                                                                                                                                       >  2:554961/1‑148 (MQ=255)
    tGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttgg                                                                                                                                      <  4:393000/149‑4 (MQ=255)
    tGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttg                                                                                                                                       <  2:896254/148‑3 (MQ=255)
    tGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAACGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttgg                                                                                                                                      <  3:762048/149‑4 (MQ=255)
      ttCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAAc                                                                                                                                    >  3:319849/1‑149 (MQ=255)
        cGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTc                                                                                                                                  <  4:598227/149‑1 (MQ=255)
           aacaacGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAAc                                                                                                                                    <  3:670226/144‑1 (MQ=255)
           aacaacGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAAc                                                                                                                                    >  4:670226/1‑144 (MQ=255)
             caacGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGc                                                                                                                              <  2:802239/148‑1 (MQ=255)
             caacGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGc                                                                                                                              <  1:495235/148‑1 (MQ=255)
              aacGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCt                                                                                                                             >  2:482396/1‑148 (MQ=255)
                  ttGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGa                                                                                                                         >  1:850747/1‑148 (MQ=255)
                  ttGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGa                                                                                                                         <  2:105009/148‑1 (MQ=255)
                   tGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCtt                                                                                                                                        >  2:121235/1‑131 (MQ=255)
                   tGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCtt                                                                                                                                        <  1:121235/132‑2 (MQ=255)
                      ccACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAt                                                                                                                    >  3:454967/1‑149 (MQ=255)
                       cACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAtt                                                                                                                   >  3:785552/1‑149 (MQ=255)
                        aCTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAtt                                                                                                                   <  2:782373/148‑1 (MQ=255)
                        aCTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAtt                                                                                                                   >  1:22714/1‑148 (MQ=255)
                          tttCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCt                                                                                                                <  4:11280/149‑1 (MQ=255)
                           ttCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTa                                                                                                               <  4:806490/149‑1 (MQ=255)
                             cagtcagCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTAcc                                                                                                             <  3:363901/149‑1 (MQ=255)
                              agtcagCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATTCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTAcc                                                                                                             <  2:22714/148‑1 (MQ=255)
                               gtcagCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGaa                                                                                                                        >  3:534463/1‑136 (MQ=255)
                               gtcagCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGaa                                                                                                                        <  4:534463/136‑1 (MQ=255)
                                tcagCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGt                                                                                                          >  4:857952/1‑149 (MQ=255)
                                tcagCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGt                                                                                                          >  3:883665/1‑149 (MQ=255)
                                    cgcCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCt                                                                                                      <  4:785552/149‑1 (MQ=255)
                                     gcCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTg                                                                                                     >  4:569021/1‑149 (MQ=255)
                                      ccACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAGGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTAc                                                                                                              >  3:852876/1‑139 (MQ=255)
                                      ccACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTAc                                                                                                              <  4:852876/139‑1 (MQ=255)
                                        aCGCCGCGAGCGGAGCGGTTCAGCAGTGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCa                                                                                                   <  2:499446/148‑1 (MQ=255)
                                        aCGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCa                                                                                                   <  1:575894/148‑1 (MQ=255)
                                         cgccgcGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCaa                                                                                                  >  2:357144/1‑148 (MQ=255)
                                            cgcgAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTc                                                                                               <  1:570929/148‑1 (MQ=255)
                                              cgAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCttt                                                                                            <  3:188259/149‑1 (MQ=255)
                                              cgAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCttt                                                                                            <  3:569021/149‑1 (MQ=255)
                                              cgAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCttt                                                                                            <  4:553265/149‑1 (MQ=255)
                                              cgAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCtt                                                                                             <  1:482396/148‑1 (MQ=255)
                                                  cggagcggTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttgggc                                                                                                                                    <  3:689062/105‑6 (MQ=255)
                                                  cggagcggTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttgggc                                                                                                                                    >  4:689062/1‑100 (MQ=255)
                                                          ttCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCa                                                                                <  4:504521/149‑1 (MQ=255)
                                                               cagcGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTa                                                                            >  2:413696/1‑148 (MQ=255)
                                                                 gcGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTAtt                                                                          >  1:22965/1‑148 (MQ=255)
                                                                 gcGTTGTTGCCAGCGTATTGTCAAGGGTTTGGAGCCGCTTAGTGAAGAGTGCAGGCGTCTGGCCACGATTCGGCGCTGCGGCGCTAGGACTCCCGCTGTGAACGATTTTTACCAGTACCTGCAATTCTTGGGCTAAAGGCCAGCTATTc                                                                         <  3:337530/149‑1 (MQ=255)
                                                                  cGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCa                                                                                >  2:362817/1‑141 (MQ=255)
                                                                  cGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCa                                                                                <  1:362817/141‑1 (MQ=255)
                                                                       ttgCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCAt                                                                    <  2:466661/148‑1 (MQ=255)
                                                                           cAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttg                                                                                                                                       >  3:81439/1‑75 (MQ=37)
                                                                           cAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttg                                                                                                                                       <  4:81439/77‑3 (MQ=37)
                                                                              cGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAAtt                                                            <  4:883665/149‑1 (MQ=255)
                                                                               gTATTTTCAAGGATTGGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTa                                                           >  4:667796/1‑149 (MQ=255)
                                                                                      cAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAt                                                                                                                    <  3:688191/85‑1 (MQ=255)
                                                                                      cAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAt                                                                                                                    >  4:688191/1‑85 (MQ=255)
                                                                                         ggATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGtt                                                  >  2:664430/1‑148 (MQ=255)
                                                                                           aTTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAg                                                                                    <  3:167633/112‑1 (MQ=255)
                                                                                           aTTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAg                                                                                    >  4:167633/1‑112 (MQ=255)
                                                                                               ggATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAAt                                                             >  4:537394/1‑131 (MQ=255)
                                                                                               ggATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAAt                                                             <  3:537394/131‑1 (MQ=255)
                                                                                               ggATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATa                                            <  1:413696/148‑1 (MQ=255)
                                                                                                 aTCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCAAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTATGGCTAAAGGCCAGCTATTCAGCATCGTTAATCATCCGCAGTTGTGATAAgc                                         <  4:569034/149‑1 (MQ=255)
                                                                                                  tCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTaa                                                              <  1:295270/127‑1 (MQ=255)
                                                                                                  tCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTaa                                                              >  2:295270/1‑127 (MQ=255)
                                                                                                  tCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAgc                                         <  1:357144/148‑1 (MQ=255)
                                                                                                   ccGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAgcgc                                       <  4:319849/149‑1 (MQ=255)
                                                                                                    cGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAgcgc                                       <  2:145308/148‑1 (MQ=255)
                                                                                                      cTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAg                                     >  1:391714/1‑148 (MQ=255)
                                                                                                            gACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTAtt                               >  2:512917/1‑148 (MQ=255)
                                                                                                              cAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTATCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTCAATTATCCGCAGTTGTGATAAGCCCAGTGTATTTCg                            <  4:566205/149‑1 (MQ=255)
                                                                                                               aaaTGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGc                           >  4:839733/1‑149 (MQ=255)
                                                                                                                    cAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCaaaa                       <  1:554961/148‑1 (MQ=255)
                                                                                                                     aGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGa                                              >  4:334169/1‑124 (MQ=255)
                                                                                                                     aGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGa                                              <  3:334169/124‑1 (MQ=255)
                                                                                                                             gCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCAt              >  2:483194/1‑148 (MQ=255)
                                                                                                                               ccAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTg                                                                                                     >  3:186478/1‑59 (MQ=255)
                                                                                                                               ccAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTg                                                                                                     <  4:186478/59‑1 (MQ=255)
                                                                                                                               ccAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTcc           >  3:108742/1‑149 (MQ=255)
                                                                                                                               ccAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTcc           >  3:166290/1‑149 (MQ=255)
                                                                                                                                 aaGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTccc          >  1:763537/1‑148 (MQ=255)
                                                                                                                                 aaGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCt         >  4:765364/1‑149 (MQ=255)
                                                                                                                                      tCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCAAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGctct    >  4:560787/1‑149 (MQ=255)
                                                                                                                                       cgccgcTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCg                            >  2:106761/1‑124 (MQ=255)
                                                                                                                                       cgccgcTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCg                            <  1:106761/124‑1 (MQ=255)
                                                                                                                                        gccgcTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTgg  <  4:108742/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
CGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGG  >  NC_000913/3206038‑3206322

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