Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,453,174 C→T 5.2% G33S (GGT→AGT)  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,1740CT 5.2% 506.3 / 12.1 153G33S (GGT→AGT) rpsJ30S ribosomal subunit protein S10
Reads supporting (aligned to +/- strand):  ref base C (70/75);  new base T (4/4);  total (74/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

Only 100 of 156 total aligned reads displayed.

TCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACC  >  NC_000913/3453027‑3453321
                                                                                                                                                   |                                                                                                                                                   
tctcGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCAcc                                                                                                                                                     <  3:667416/148‑1 (MQ=255)
tctcGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCAcc                                                                                                                                                     <  4:730513/148‑1 (MQ=255)
 ctcGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAg                                                                                                                                                   >  1:150943/1‑149 (MQ=255)
  tcGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGt                                                                                                                                                  <  1:635022/149‑1 (MQ=255)
     gTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTgcg                                                                                                                                               <  1:1020945/149‑1 (MQ=255)
         gCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTCTTGCCGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCtt                                                                                                                                            <  4:1024183/148‑1 (MQ=255)
          cTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACTGGTgcgcagg                                                                                                                                           <  3:597486/148‑4 (MQ=255)
          cTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTgg                                                                                                                                          <  2:40150/149‑1 (MQ=255)
           tCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACTGGTgcgcag                                                                                                                                            <  2:5640/146‑3 (MQ=255)
           tCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGc                                                                                                                                         >  2:837095/1‑149 (MQ=255)
             aaCGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAg                                                                                                                                       >  1:236380/1‑149 (MQ=255)
              aCGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAg                                                                                                                                       >  4:468473/1‑148 (MQ=255)
              aCGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGGGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGGGCGCTTGGCAg                                                                                                                                       >  4:556739/1‑148 (MQ=255)
                  tGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGtctc                                                                                                                                   <  4:887350/148‑1 (MQ=255)
                   gTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTcga                                                                                                                                 <  1:187615/149‑1 (MQ=255)
                      aaCCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCtt                                                                                                                                            >  1:957717/1‑135 (MQ=255)
                      aaCCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTcgacga                                                                                                                              <  2:546573/149‑1 (MQ=255)
                      aaCCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTcgacga                                                                                                                              >  1:7747/1‑149 (MQ=255)
                        ccAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGAtt                                                                                                                            >  2:805300/1‑149 (MQ=255)
                         cAGACGCAAGTGAGTACGGATTTCGTACCGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTgg                                                                                                                                          <  2:282058/134‑1 (MQ=255)
                         cAGACGCAAGTGAGTACGGATTTCGTACCGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTgg                                                                                                                                          >  1:282058/1‑134 (MQ=255)
                              gCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCcgcg                                                                                                                      >  2:992739/1‑149 (MQ=255)
                                 aGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGtt                                                                                                                   <  1:421691/149‑1 (MQ=255)
                                   tGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGc                                                                                                                 >  1:345090/1‑149 (MQ=255)
                                    gAGTACGGATTTCGTACTGATCGCGCGCGTCTTCGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGc                                                                                                                 <  3:210827/148‑1 (MQ=255)
                                    caGTACGAAATCCGTACTGATCGCGCGCGTCTTGGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGGTCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGc                                                                                                                 >  4:51097/2‑148 (MQ=255)
                                       tACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACTGGTgcgcagg                                                                                                                                           >  2:845076/1‑116 (MQ=255)
                                       tACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACTGGTgcgcagg                                                                                                                                           <  1:845076/119‑4 (MQ=255)
                                       tACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTga                                                                                                             >  1:581801/1‑149 (MQ=255)
                                       tACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTga                                                                                                             >  1:241242/1‑149 (MQ=255)
                                        aCGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTgat                                                                                                            <  1:293472/149‑1 (MQ=255)
                                        aCGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTga                                                                                                             >  3:87918/1‑148 (MQ=255)
                                            aTTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTcgac                                                                                                                                <  1:300798/125‑1 (MQ=255)
                                               tCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGTAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCa                                                                                                      >  4:690397/1‑148 (MQ=255)
                                                   aCTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACga                                                                                                 >  2:86825/1‑149 (MQ=255)
                                                   aCTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACg                                                                                                  <  3:266598/148‑1 (MQ=255)
                                                    cTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACgat                                                                                                >  2:377506/1‑149 (MQ=255)
                                                      gATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGAGCTGCGCACCAGTCCGCTTGGCAGTCTCGACGATTTCAGCGGTTGCTTGATCGATCAGACgatg                                                                                               >  4:936008/1‑148 (MQ=255)
                                                      gATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTga                                                                                                             >  3:912575/1‑134 (MQ=255)
                                                      gATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTga                                                                                                             <  4:912575/134‑1 (MQ=255)
                                                      gATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGa                                                                                                    >  2:497831/1‑143 (MQ=255)
                                                       aTCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACgatgat                                                                                             <  2:1003977/149‑1 (MQ=255)
                                                       aTCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACgatgat                                                                                             <  2:100953/149‑1 (MQ=255)
                                                       aTCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACgatga                                                                                              >  4:371945/1‑148 (MQ=255)
                                                        tCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACgatgat                                                                                             >  4:766104/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)
                                                              gcgTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGc                                                                                      >  2:293309/1‑149 (MQ=255)
                                                              gcgTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGc                                                                                      >  2:632025/1‑149 (MQ=255)
                                                                 tCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCttt                                                                                   <  2:345090/149‑1 (MQ=255)
                                                                  cTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCttt                                                                                   >  4:679796/1‑148 (MQ=255)
                                                                    ttgttgACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCa                                                                                 >  4:614667/1‑148 (MQ=255)
                                                                      gttgACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGc                                                                              >  1:1033119/1‑149 (MQ=255)
                                                                      gttgACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGc                                                                              >  1:319618/1‑149 (MQ=255)
                                                                         gACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGa                                                                           >  2:473890/1‑149 (MQ=255)
                                                                         gACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGa                                                                           <  1:486079/149‑1 (MQ=255)
                                                                            gTGCGGGGAGATCAGAACAGTGAAGCGCACTTTGCGTGTCGGCAGCGGGATCGGACCAGGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATCATCAAACGCTTTCAGGCGCATa                                                                         <  4:719247/148‑1 (MQ=255)
                                                                             tGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATAc                                                                        >  3:230761/1‑148 (MQ=255)
                                                                              gCGGGGAGATCAGAACAGTGAAGCGGTCTTTGCGTGTCGGCAGCGGGATCGGAGCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTGGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACg                                                                       >  3:1021173/1‑148 (MQ=255)
                                                                               cGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGa                                                                     >  2:91409/1‑149 (MQ=255)
                                                                                 gggAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGa                                                                           >  1:524581/1‑141 (MQ=255)
                                                                                  ggAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGAtt                                                                   <  3:564257/148‑1 (MQ=255)
                                                                                   gagaTCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTc                                                                  >  4:120594/1‑148 (MQ=255)
                                                                                      aTCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTg                                                              >  1:559124/1‑149 (MQ=255)
                                                                                      aTCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTg                                                              >  2:448746/1‑149 (MQ=255)
                                                                                       tCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTg                                                              >  4:811941/1‑148 (MQ=255)
                                                                                        cAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTgg                                                             >  3:774446/1‑148 (MQ=255)
                                                                                         aGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGt                                                            <  3:945420/148‑1 (MQ=255)
                                                                                          gAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGtt                                                           <  3:800841/148‑1 (MQ=255)
                                                                                          gAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTc                                                          >  2:109132/1‑149 (MQ=255)
                                                                                           aaCAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTg                                                                                                                                           <  4:1015505/67‑1 (MQ=255)
                                                                                           aaCAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTg                                                                                                                                           >  3:1015505/1‑67 (MQ=255)
                                                                                              aGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACTGGTgcgcag                                                                                                                                            >  1:429304/1‑61 (MQ=255)
                                                                                              aGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGc                                                       <  4:912401/148‑1 (MQ=255)
                                                                                               gTGAAGCTCTCTTTGCGTTTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCAt                                                     <  2:269162/149‑1 (MQ=255)
                                                                                               gTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCAt                                                     <  2:1070268/149‑1 (MQ=255)
                                                                                                      gcTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGAcc                                               <  3:51097/148‑1 (MQ=255)
                                                                                                      gcTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCa                                              <  2:172697/149‑1 (MQ=255)
                                                                                                        tctTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCaga                                            <  1:805300/149‑1 (MQ=255)
                                                                                                          tttGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCagag                                           <  4:17786/148‑1 (MQ=255)
                                                                                                               gtgtCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCaa                                     >  2:419633/1‑149 (MQ=255)
                                                                                                                tgtCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAgg                                                                               >  3:708484/1‑106 (MQ=255)
                                                                                                                tgtCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCaa                                     >  4:565688/1‑148 (MQ=255)
                                                                                                                    ggCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAattat                                 <  4:580207/148‑1 (MQ=255)
                                                                                                                      cAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTAttt                               <  4:801165/148‑1 (MQ=255)
                                                                                                                         cGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTtat                            <  3:586022/148‑1 (MQ=255)
                                                                                                                            gATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAAc                        >  1:558141/1‑149 (MQ=255)
                                                                                                                                 gACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAattatt                                <  4:152536/136‑1 (MQ=255)
                                                                                                                                 gACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAattatt                                >  3:152536/1‑136 (MQ=255)
                                                                                                                                 gACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAatta                                  <  1:181765/134‑1 (MQ=255)
                                                                                                                                 gACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAatta                                  >  2:181765/1‑134 (MQ=255)
                                                                                                                                    cACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCGGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGTGCTCCAATTATTTTATAAACGAAAATg                 <  4:650288/148‑1 (MQ=255)
                                                                                                                                    cACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAattatt                                >  1:496411/1‑133 (MQ=255)
                                                                                                                                    cACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGCTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATg                 <  3:665114/148‑1 (MQ=255)
                                                                                                                                       ggACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGa                                                                     <  1:963809/93‑1 (MQ=255)
                                                                                                                                       ggACCTGCGCACCAGTGCGCTTGGCAGTCTAGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGa                                                                     >  2:963809/1‑93 (MQ=255)
                                                                                                                                             gcgcACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTActcctc       <  2:558141/149‑1 (MQ=255)
                                                                                                                                                cACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCac     <  4:254368/148‑2 (MQ=255)
                                                                                                                                                   cAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGAcc  <  3:190724/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCTCGGTTGGCTCAACGATGTCAACCAGACGCAAGTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACC  >  NC_000913/3453027‑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: