Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,037,957 G→T 5.6% P50P (CCC→CCA recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,037,9570GT 5.6% 512.6 / 15.7 143P50P (CCC→CCArecJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base G (70/65);  new base T (4/4);  total (74/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

Only 100 of 147 total aligned reads displayed.

AACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCA  >  NC_000913/3037809‑3038105
                                                                                                                                                    |                                                                                                                                                    
        gagcgCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCCTGgcagcaa                                                                                                                                              <  4:517272/149‑2 (MQ=255)
          gcgcATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATAc                                                                                                                                            <  4:1145833/149‑1 (MQ=255)
          gcgcATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATAc                                                                                                                                            <  4:1460715/149‑1 (MQ=255)
                       aCGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCg                                                                                                                               >  1:1247261/1‑149 (MQ=255)
                           ttAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTaa                                                                                                                                      <  3:227526/138‑1 (MQ=255)
                           ttAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTcc                                                                                                                           >  3:1501273/1‑149 (MQ=255)
                           ttAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTcc                                                                                                                           >  2:38508/1‑149 (MQ=255)
                            tAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCa                                                                                                                          <  4:1501273/149‑1 (MQ=255)
                                    tGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTgc                                                                                                                  <  1:1182771/149‑1 (MQ=255)
                                    tGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTgc                                                                                                                  <  1:1225216/149‑1 (MQ=255)
                                       tGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCa                                                                                                               >  3:86154/1‑149 (MQ=255)
                                          tCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTg                                                                                                            >  2:1434149/1‑149 (MQ=255)
                                               ccGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTAc                                                                                                       >  1:106664/1‑149 (MQ=255)
                                                     gcgtcgaaaaggccgttgagatcctttacaacgcttttcgcgaaggaacgcggattattgtggtcggtgatttcgacgccgCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCg                                                                                                 >  2:260004‑M2/82‑149 (MQ=255)
                                                      cgtcgaaaaggccgttgagatcctttacaacgcttttcgcgaaggaacgcggattattgtggtcggtgatttcgacgccgCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTc                                                                                                                            >  1:1662787‑M2/81‑121 (MQ=255)
                                                      cgtcgaaaaggccgttgagatcctttacaacgcttttcgcgaaggaacgcggattattgtggtcggtgatttcgacgccgCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTc                                                                                                                            <  2:1662787‑M2/41‑1 (MQ=255)
                                                       gtcgaaaaggccgttgagatcctttacaacgcttttcgcgaaggaacgcggattattgtggtcggtgatttcgacgccgCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCtt                                                                                                                    <  4:321806‑M2/49‑1 (MQ=255)
                                                       gtcgaaaaggccgttgagatcctttacaacgcttttcgcgaaggaacgcggattattgtggtcggtgatttcgacgccgCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCtt                                                                                                                    >  3:321806‑M2/80‑128 (MQ=255)
                                                       GTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCg                                                                                                                               >  2:1141797‑M1/1‑117 (MQ=255)
                                                       GTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCg                                                                                                                               <  1:1141797‑M1/117‑1 (MQ=255)
                                                          GAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGgcagca                                                                                                                                               >  1:634779‑M1/1‑98 (MQ=255)
                                                           AAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCCtggcagcaac                                                                                                                                             <  2:708508‑M1/99‑3 (MQ=255)
                                                           AAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCCtggcagcaac                                                                                                                                             >  1:708508‑M1/1‑97 (MQ=255)
                                                              tCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTATTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGta                                                                                        >  4:1802465/1‑149 (MQ=255)
                                                                aCCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGtata                                                                                      >  1:127769/1‑149 (MQ=255)
                                                                 ccGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGtat                                                                                       >  3:1465294/1‑147 (MQ=255)
                                                                    aCCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACg                                                                                  >  3:1695402/1‑149 (MQ=255)
                                                                     ccACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACgg                                                                                 <  1:1075965/149‑1 (MQ=255)
                                                                       acaATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGtata                                                                                      >  4:1460482/1‑142 (MQ=255)
                                                                         aataatCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTaa                                                                                                                                      >  3:1812787/1‑92 (MQ=255)
                                                                         aataatCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTaa                                                                                                                                      <  4:1812787/92‑1 (MQ=255)
                                                                         aataatCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTc                                                                                                                      >  2:1199813/1‑108 (MQ=255)
                                                                         aataatCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTc                                                                                                                      <  1:1199813/108‑1 (MQ=255)
                                                                                cgcgTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGta                                                                                        >  4:538807/1‑131 (MQ=255)
                                                                                 gcgTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGa                                                                     <  3:1802465/149‑1 (MQ=255)
                                                                                 gcgTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGa                                                                     <  1:38508/149‑1 (MQ=255)
                                                                                   gTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAgg                                                                   <  2:1247261/149‑1 (MQ=255)
                                                                                        ttCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAAtt                                                              >  2:1856965/1‑149 (MQ=255)
                                                                                            cgAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGta                                                                                        <  3:977228/119‑1 (MQ=255)
                                                                                            cgAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGta                                                                                        >  4:977228/1‑119 (MQ=255)
                                                                                              gaaggaacgcggattattgtggtcggtgatttcgacgccgCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCgg                                                        <  1:260004‑M2/109‑1 (MQ=255)
                                                                                              aaaaGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCg                                                         <  3:884744/148‑1 (MQ=255)
                                                                                                 aGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTa                                                                                                        <  4:34711/98‑1 (MQ=255)
                                                                                                  gCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCaa                                                    <  1:1434149/149‑1 (MQ=255)
                                                                                                        tAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAAc                                                                                   >  2:1272813/1‑112 (MQ=255)
                                                                                                         aaaGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGcc                                             >  3:393187/1‑149 (MQ=255)
                                                                                                          aaGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCaa                                                                         >  2:1699577/1‑120 (MQ=255)
                                                                                                          aaGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCaa                                                                         <  1:1699577/120‑1 (MQ=255)
                                                                                                          aaGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAg                                                                    <  2:1562188/125‑1 (MQ=255)
                                                                                                          aaGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAg                                                                    >  1:1810925/1‑125 (MQ=255)
                                                                                                            ggATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGtt                                          >  2:790473/1‑149 (MQ=255)
                                                                                                            ggATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGtt                                          >  2:949474/1‑149 (MQ=255)
                                                                                                             gATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGGACTGCGTACTCTCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGttt                                         >  2:1344973/1‑149 (MQ=255)
                                                                                                             gATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGttt                                         >  3:1273093/1‑149 (MQ=255)
                                                                                                              aTCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTc                                        <  1:949474/149‑1 (MQ=255)
                                                                                                               tctcAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAg                                                                    <  3:1393101/120‑1 (MQ=255)
                                                                                                                ctcAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGgc                                                                                >  2:1178017/1‑107 (MQ=255)
                                                                                                                ctcAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGgc                                                                                <  3:986550/107‑1 (MQ=255)
                                                                                                                ctcAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAg                                                                    <  4:152173/119‑1 (MQ=255)
                                                                                                                 tcAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGa                                                                     >  4:1680400/1‑117 (MQ=255)
                                                                                                                 tcAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATc                                     <  4:1695402/149‑1 (MQ=255)
                                                                                                                 tcAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATc                                     <  4:1273093/149‑1 (MQ=255)
                                                                                                                   aaCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCaa                                                    <  1:1453125/132‑1 (MQ=255)
                                                                                                                   aaCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGa                                   >  3:452913/1‑149 (MQ=255)
                                                                                                                   aaCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGa                                   >  1:1298431/1‑149 (MQ=255)
                                                                                                                   aaCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGa                                   <  2:566095/149‑1 (MQ=255)
                                                                                                                    aCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGcc                                             <  2:1819238/138‑1 (MQ=255)
                                                                                                                     cGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACt                                 <  2:1693229/149‑1 (MQ=255)
                                                                                                                       gCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTgc                                                                                                                  >  4:356413/1‑66 (MQ=255)
                                                                                                                        ccTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCa                              >  1:540871/1‑149 (MQ=255)
                                                                                                                              cgacgccgCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGc                                                          <  3:686810‑M2/107‑1 (MQ=255)
                                                                                                                              cgacgccgCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGcga                        <  2:106664‑M2/141‑1 (MQ=255)
                                                                                                                                   ccgCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTgc                                                                                                                  >  4:871158‑M2/4‑54 (MQ=255)
                                                                                                                                     gCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAg                   >  3:72735‑M2/2‑147 (MQ=255)
                                                                                                                                     gCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAg                   <  4:72735‑M2/146‑1 (MQ=255)
                                                                                                                                     gCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGtt                 <  2:1208741‑M2/148‑1 (MQ=255)
                                                                                                                                     gCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGtt                 <  4:393187‑M2/148‑1 (MQ=255)
                                                                                                                                      cTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTg                >  2:227849/1‑149 (MQ=255)
                                                                                                                                       tCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGGACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGc                                                          >  4:1107048/1‑106 (MQ=255)
                                                                                                                                       tCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGGACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGc                                                          <  3:1107048/106‑1 (MQ=255)
                                                                                                                                         aGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTAt             <  1:790473/149‑1 (MQ=255)
                                                                                                                                          tatgctgcCCTGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATc                                     <  1:1686368/122‑1 (MQ=255)
                                                                                                                                          tatgctgcCCTGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGAc                                  <  4:716122/125‑1 (MQ=255)
                                                                                                                                          tatgctgcCCTGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGAc                                  >  3:716122/3‑127 (MQ=255)
                                                                                                                                          tatgctgcCCTGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATc            >  3:1771534/3‑149 (MQ=255)
                                                                                                                                          gTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCgca                                                                             <  4:1663867/84‑1 (MQ=255)
                                                                                                                                          gTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCgca                                                                             >  3:1663867/1‑84 (MQ=255)
                                                                                                                                          gTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATc            <  4:937507/149‑1 (MQ=255)
                                                                                                                                           atgctgcCCTGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGcgacga                     >  3:35563/2‑139 (MQ=255)
                                                                                                                                             gctgcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGcga                        >  4:1212891/1‑134 (MQ=255)
                                                                                                                                             gctgcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGcga                        <  3:1212891/134‑1 (MQ=255)
                                                                                                                                             gctgcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCtgt         >  3:351695/1‑149 (MQ=255)
                                                                                                                                             gctgcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCtgt         >  3:1541572/1‑149 (MQ=255)
                                                                                                                                              ctgcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCtgtt        <  1:227849/149‑1 (MQ=255)
                                                                                                                                                 ccAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTgttgtt     >  1:786563/1‑149 (MQ=255)
                                                                                                                                                  cAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGcc                                             >  3:868532/1‑108 (MQ=255)
                                                                                                                                                  cAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGcc                                             <  4:868532/108‑1 (MQ=255)
                                                                                                                                                   aGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCtggtgt      >  1:193037/1‑146 (MQ=255)
                                                                                                                                                   aGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTc   >  3:891629/1‑149 (MQ=255)
                                                                                                                                                    gggCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCa  <  2:137908/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCA  >  NC_000913/3037809‑3038105

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: