Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,453,175 A→G 5.2% T32T (ACT→ACC rpsJ ← 30S ribosomal subunit protein S10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,453,1750AG 5.2% 507.9 / 12.5 154T32T (ACT→ACCrpsJ30S ribosomal subunit protein S10
Reads supporting (aligned to +/- strand):  ref base A (70/76);  new base G (4/4);  total (74/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

Only 100 of 154 total aligned reads displayed.

CTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACC  >  NC_000913/3453028‑3453321
                                                                                                                                                   |                                                                                                                                                  
ctcGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAg                                                                                                                                                   >  1:150943/1‑149 (MQ=255)
 tcGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGt                                                                                                                                                  <  1:635022/149‑1 (MQ=255)
    gTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTgcg                                                                                                                                               <  1:1020945/149‑1 (MQ=255)
    gTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTgcg                                                                                                                                               >  1:546573/1‑149 (MQ=255)
        gCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTg                                                                                                                                           >  1:722731/1‑149 (MQ=255)
          tCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACTGGTgcgcag                                                                                                                                            >  1:5640/1‑144 (MQ=255)
          tCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACTGGTgcgcag                                                                                                                                            <  2:5640/146‑3 (MQ=255)
           cAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGc                                                                                                                                         >  4:564257/1‑148 (MQ=255)
            aaCGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAg                                                                                                                                       >  1:236380/1‑149 (MQ=255)
              cGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGtc                                                                                                                                     <  2:1016817/149‑1 (MQ=255)
                 tGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGtctc                                                                                                                                   <  4:887350/148‑1 (MQ=255)
                  gTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTcga                                                                                                                                 <  1:187615/149‑1 (MQ=255)
                     aaCCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCtt                                                                                                                                            <  2:957717/135‑1 (MQ=255)
                     aaCCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTcgacga                                                                                                                              <  2:546573/149‑1 (MQ=255)
                       ccAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGAtt                                                                                                                            >  2:805300/1‑149 (MQ=255)
                        cAGACGCAAGTGAGTACGGATTTCGTACCGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTgg                                                                                                                                          >  1:282058/1‑134 (MQ=255)
                        cAGACGCAAGTGAGTACGGATTTCGTACCGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTgg                                                                                                                                          <  2:282058/134‑1 (MQ=255)
                          gACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTc                                                                                                                          <  4:546321/148‑1 (MQ=255)
                                aGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGtt                                                                                                                   <  1:421691/149‑1 (MQ=255)
                                  tGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTg                                                                                                                  >  3:580207/1‑148 (MQ=255)
                                  tGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGc                                                                                                                 >  1:345090/1‑149 (MQ=255)
                                   gAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGc                                                                                                                 <  3:732099/148‑1 (MQ=255)
                                   gAGTACGGATTTCGTACTGATCGCGCGCGTCTTCGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGc                                                                                                                 <  3:210827/148‑1 (MQ=255)
                                   caGTACGAAATCCGTACTGATCGCGCGCGTCTTGGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGGTCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGc                                                                                                                 >  4:51097/2‑148 (MQ=255)
                                    aGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCtt                                                                                                               <  1:837095/149‑1 (MQ=255)
                                      tACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTga                                                                                                             >  1:241242/1‑149 (MQ=255)
                                           aTTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTcgac                                                                                                                                <  1:300798/125‑1 (MQ=255)
                                           aTTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTcgac                                                                                                                                >  2:300798/1‑125 (MQ=255)
                                            tttCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTcgacga                                                                                                                              <  1:995386/126‑1 (MQ=255)
                                            tttCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTgatcgat                                                                                                        <  4:5999/148‑1 (MQ=255)
                                            tttCGAACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTcgacga                                                                                                                              >  2:995386/1‑126 (MQ=255)
                                              tCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGTAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCa                                                                                                      >  4:690397/1‑148 (MQ=255)
                                                  aCTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACga                                                                                                 >  2:86825/1‑149 (MQ=255)
                                                  aCTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACga                                                                                                 <  1:304755/149‑1 (MQ=255)
                                                  aCTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACg                                                                                                  <  3:266598/148‑1 (MQ=255)
                                                     gATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGAGCTGCGCACCAGTCCGCTTGGCAGTCTCGACGATTTCAGCGGTTGCTTGATCGATCAGACgatg                                                                                               >  4:936008/1‑148 (MQ=255)
                                                     gATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTga                                                                                                             <  4:912575/134‑1 (MQ=255)
                                                     gATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTga                                                                                                             >  3:912575/1‑134 (MQ=255)
                                                      aTCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACgatgat                                                                                             <  2:100953/149‑1 (MQ=255)
                                                      aTCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACgatga                                                                                              >  4:371945/1‑148 (MQ=255)
                                                       tCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACgatgat                                                                                             >  4:766104/1‑148 (MQ=255)
                                                          cgcgcgTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGGTCGGACAACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGATTGCTTGATCGATCAGACGATGATCaa                                                                                          >  4:74815/1‑148 (MQ=255)
                                                          cgcgcgTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCaaa                                                                                         >  1:118177/1‑149 (MQ=255)
                                                          cgcgcgTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCaa                                                                                          >  4:214065/1‑148 (MQ=255)
                                                            cgcgTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAAc                                                                                        >  4:875812/1‑148 (MQ=255)
                                                              cgTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGc                                                                                      <  3:208264/148‑1 (MQ=255)
                                                                tCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCttt                                                                                   <  2:345090/149‑1 (MQ=255)
                                                                 cTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTc                                                                                  <  2:233792/149‑1 (MQ=255)
                                                                   ttgttgACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCa                                                                                 >  4:614667/1‑148 (MQ=255)
                                                                     gttgACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGc                                                                              >  1:1033119/1‑149 (MQ=255)
                                                                        gACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGa                                                                           <  1:486079/149‑1 (MQ=255)
                                                                        gACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGa                                                                           >  2:473890/1‑149 (MQ=255)
                                                                            tGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATAc                                                                        >  4:416685/1‑148 (MQ=255)
                                                                            tGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATAc                                                                        >  3:230761/1‑148 (MQ=255)
                                                                                gggAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGa                                                                           >  1:524581/1‑141 (MQ=255)
                                                                                gggAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGa                                                                           <  2:524581/141‑1 (MQ=255)
                                                                                     aTCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTg                                                              >  1:559124/1‑149 (MQ=255)
                                                                                      tCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTg                                                              >  3:801165/1‑148 (MQ=255)
                                                                                      tCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTg                                                              >  4:103524/1‑148 (MQ=255)
                                                                                      tCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTg                                                              >  4:811941/1‑148 (MQ=255)
                                                                                       cAGAACAGTGAGGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACTATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTgg                                                             >  4:665114/1‑148 (MQ=255)
                                                                                         gAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGtt                                                           <  3:800841/148‑1 (MQ=255)
                                                                                         gAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTc                                                          >  2:109132/1‑149 (MQ=255)
                                                                                          aaCAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTg                                                                                                                                           >  3:1015505/1‑67 (MQ=255)
                                                                                          aaCAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTg                                                                                                                                           <  4:1015505/67‑1 (MQ=255)
                                                                                          aaCAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCt                                                         >  2:141028/1‑149 (MQ=255)
                                                                                           aCAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCt                                                         >  4:1022288/1‑148 (MQ=255)
                                                                                             aGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACTGGTgcgcag                                                                                                                                            <  2:429304/63‑3 (MQ=255)
                                                                                             aGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACTGGTgcgcag                                                                                                                                            >  1:429304/1‑61 (MQ=255)
                                                                                             aGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGc                                                       <  4:912401/148‑1 (MQ=255)
                                                                                             aGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCa                                                      <  2:39308/149‑1 (MQ=255)
                                                                                                gcaGCGCTCTTTGCGTGTGGGCAGCGGGATCGTTCCCCGGGCCTGCGCACCAGTGCGCTTGGCAGTCTCGACTATTTCCGCGGTTGCTTGATTGATCATCCGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATgt                                                   <  2:840588/147‑2 (MQ=255)
                                                                                                gAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATga                                                   <  2:271050/149‑1 (MQ=255)
                                                                                                  aGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATgaga                                                 <  1:293309/149‑1 (MQ=255)
                                                                                                     gcTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCGGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCa                                              <  2:746079/149‑1 (MQ=255)
                                                                                                     gcTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGAcc                                               <  3:51097/148‑1 (MQ=255)
                                                                                                       tctTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCaga                                            <  1:805300/149‑1 (MQ=255)
                                                                                                              gtgtCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCaa                                     >  2:419633/1‑149 (MQ=255)
                                                                                                              gtgtCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCaa                                     <  2:559124/149‑1 (MQ=255)
                                                                                                               tgtCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAgg                                                                               >  3:708484/1‑106 (MQ=255)
                                                                                                               tgtCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAgg                                                                               <  4:708484/106‑1 (MQ=255)
                                                                                                               tgtCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCaa                                     >  4:565688/1‑148 (MQ=255)
                                                                                                                 tCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAatt                                   >  3:608403/1‑148 (MQ=255)
                                                                                                                     cAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTAttt                               <  4:801165/148‑1 (MQ=255)
                                                                                                                     cAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTAttt                               <  3:103524/148‑1 (MQ=255)
                                                                                                                      aGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTAtttt                              <  4:230761/148‑1 (MQ=255)
                                                                                                                        cGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTtat                            <  3:586022/148‑1 (MQ=255)
                                                                                                                          ggATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATaaa                         <  1:91409/149‑1 (MQ=255)
                                                                                                                          ggATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATaa                          <  3:936008/148‑1 (MQ=255)
                                                                                                                                gACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAattatt                                >  3:152536/1‑136 (MQ=255)
                                                                                                                                gACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAattatt                                <  4:152536/136‑1 (MQ=255)
                                                                                                                                gACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAatta                                  >  2:181765/1‑134 (MQ=255)
                                                                                                                                gACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAatta                                  <  1:181765/134‑1 (MQ=255)
                                                                                                                                   cACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAattatt                                <  2:496411/133‑1 (MQ=255)
                                                                                                                                   cACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAattatt                                >  1:496411/1‑133 (MQ=255)
                                                                                                                                   cACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGa                >  2:853876/1‑149 (MQ=255)
                                                                                                                                   cACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGCTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATg                 <  3:665114/148‑1 (MQ=255)
                                                                                                                                      ggACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGa                                                                     <  1:963809/93‑1 (MQ=255)
                                                                                                                                      ggACCTGCGCACCAGTGCGCTTGGCAGTCTAGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGa                                                                     >  2:963809/1‑93 (MQ=255)
                                                                                                                                          aagcgcACTGGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTActc          >  4:158379/3‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACC  >  NC_000913/3453028‑3453321

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