Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 224,910 2 bp→GT noncoding (1140‑1141/1542 nt) rrsH → 16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913224,9100CG91.1% 99.4 / 101.6 438noncoding (1140/1542 nt)rrsH16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (19/20);  new base G (174/225);  total (193/245)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01
*NC_000913224,9110CT91.1% 105.5 / 86.4 438noncoding (1141/1542 nt)rrsH16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (19/20);  new base T (174/225);  total (193/245)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-02

Only 100 of 439 total aligned reads displayed.

AACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAG  >  NC_000913/224774‑225028
                                                                                                                                         ||                                                                                                                     
                         ttCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGAt                                                                                            >  1:1897735/1‑140 (MQ=34)
                                cattcaagaCAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCgg                                                                                     <  2:1319051/134‑1 (MQ=2)
                                 cTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGGGACTGCCAGTGATAAACTgga                                                                                    >  1:1771882/1‑140 (MQ=18)
                                 attcaagaCAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCgga                                                                                    <  2:1518746/135‑1 (MQ=2)
                                  ttcaagaCAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCggag                                                                                   >  1:243978/5‑140 (MQ=2)
                                   tcaagaCAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCtgagg                                                                                  <  1:1447468/137‑1 (MQ=2)
                                    tgAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTggagga                                                                                 >  2:1620517/1‑140 (MQ=32)
                                      agaCAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAg                                                                               >  2:1042822/1‑140 (MQ=2)
                                       gaCAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAgg                                                                              <  2:831425/140‑1 (MQ=2)
                                        aCAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTTAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGt                                                                             >  2:1220722/1‑140 (MQ=2)
                                          aGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCggagg                                                                                  <  1:747085/133‑1 (MQ=2)
                                           ggttgtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTggg                                                                          >  2:652170/6‑140 (MQ=2)
                                           ggtggtGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTggg                                                                          >  1:1253238/1‑140 (MQ=2)
                                            gtggtgCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTggg                                                                          >  2:1539596/1‑140 (MQ=2)
                                            gtggtgCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                                         >  1:2293851/1‑140 (MQ=2)
                                            gTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTgggg                                                                         >  2:1706359/1‑140 (MQ=32)
                                             tggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGa                                                                        >  1:451761/4‑140 (MQ=2)
                                             tggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCAGCAACGAGCGCAACCCTTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGa                                                                        >  1:1421928/4‑140 (MQ=2)
                                              ggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGAt                                                                       <  2:451761/138‑1 (MQ=2)
                                                 gcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGAc                                                                    >  2:673291/1‑140 (MQ=11)
                                                  cATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACg                                                                   <  2:1026543/140‑1 (MQ=2)
                                                  cATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACg                                                                   >  2:2162061/1‑140 (MQ=32)
                                                      gttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaa                                                               <  1:215079/138‑1 (MQ=2)
                                                      gttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGAt                                                                       <  1:378154/130‑1 (MQ=2)
                                                      gttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaa                                                               <  2:616548/138‑1 (MQ=2)
                                                      gttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaa                                                               >  2:232811/3‑140 (MQ=2)
                                                      gttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCAGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaa                                                               >  1:1312796/3‑140 (MQ=2)
                                                      gttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaa                                                               >  1:1800566/3‑140 (MQ=11)
                                                      gttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaa                                                               >  1:2023120/3‑140 (MQ=11)
                                                      gttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCATATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACg                                                                   <  2:425964/134‑1 (MQ=2)
                                                      gttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCATATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACg                                                                   >  1:425964/3‑136 (MQ=2)
                                                       ttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaa                                                               >  2:929554/2‑140 (MQ=2)
                                                       ttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGAt                                                                       <  1:1066580/130‑1 (MQ=2)
                                                       ttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAATTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaaa                                                              <  2:1485163/139‑2 (MQ=2)
                                                        tGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaaa                                                              >  2:1192789/1‑139 (MQ=2)
                                                        tGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaaa                                                              <  1:2079/140‑2 (MQ=2)
                                                          tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtc                                                            >  2:1064/1‑140 (MQ=2)
                                                          tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTTCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtca                                                           >  1:572126/1‑140 (MQ=2)
                                                          tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtca                                                           >  2:78836/1‑140 (MQ=11)
                                                          tcgtcAGCTCGTGTAGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtca                                                           >  1:99624/1‑140 (MQ=2)
                                                           cgtcAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcat                                                          <  2:1897735/140‑1 (MQ=32)
                                                           cgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtcat                                                          <  1:673291/140‑1 (MQ=2)
                                                           cgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtcat                                                          <  2:1443907/140‑1 (MQ=2)
                                                            gtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtcat                                                          >  1:1210183/1‑140 (MQ=2)
                                                               aGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGTGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATg                                                      <  1:1001251/140‑1 (MQ=2)
                                                               aGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATg                                                      <  2:573227/140‑1 (MQ=2)
                                                               aGCTCGTGTAGTGAGATGTTGGGTTAAGTACCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATg                                                      <  1:760335/140‑1 (MQ=2)
                                                                gCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATg                                                      <  2:326425/140‑1 (MQ=2)
                                                                gCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                                     <  1:1530831/140‑2 (MQ=2)
                                                                gCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                                     <  2:1820196/140‑2 (MQ=2)
                                                                gCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                                     <  2:411908/140‑2 (MQ=2)
                                                                gCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                                     <  1:397489/140‑2 (MQ=2)
                                                                 cTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgcc                                                    >  1:104563/1‑138 (MQ=2)
                                                                 cTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgcc                                                    >  1:824503/1‑138 (MQ=2)
                                                                 cTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgcc                                                    >  1:899819/1‑138 (MQ=2)
                                                                     tgtcgtGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                                <  1:913327/140‑1 (MQ=2)
                                                                     tgtcgtGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                                <  2:549385/140‑1 (MQ=2)
                                                                      gtcgtgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                                     >  1:53369/4‑133 (MQ=2)
                                                                      gtcgtgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                                     <  2:53369/131‑2 (MQ=2)
                                                                      gtcgtgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTa                                               <  2:966293/137‑1 (MQ=2)
                                                                        cgtgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTatg                                             <  1:984957/139‑3 (MQ=2)
                                                                         gtgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTatg                                             >  1:117973/1‑138 (MQ=2)
                                                                          tgAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGAc                                           <  1:1620517/140‑1 (MQ=32)
                                                                           gagaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcc                                          <  2:1411019/137‑1 (MQ=2)
                                                                           gagaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcc                                          >  2:392565/4‑140 (MQ=2)
                                                                           gagaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcc                                          <  2:1134330/137‑1 (MQ=2)
                                                                                gTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTAAAATCATCATGCCCCTTATGAcct                                         >  1:1301372/1‑135 (MQ=2)
                                                                                gTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCt                                    >  2:585829/1‑140 (MQ=2)
                                                                                   gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATAATGCCCCTTATGACCTGGGCTac                                  <  1:497815/140‑1 (MQ=2)
                                                                                    ggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                                  <  2:117973/140‑1 (MQ=2)
                                                                                    ggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCt                                    >  1:1347507/1‑137 (MQ=2)
                                                                                    ggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCt                                    <  2:1347507/137‑1 (MQ=2)
                                                                                    ggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATAATGCCCCTTATGACCTGGGCTaca                                 >  2:1265682/1‑140 (MQ=2)
                                                                                     gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                                <  1:1232537/140‑1 (MQ=2)
                                                                                     gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTTACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                                <  1:642125/140‑1 (MQ=2)
                                                                                     gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                                <  1:585029/140‑1 (MQ=2)
                                                                                     gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                                <  1:115504/140‑1 (MQ=2)
                                                                                      ttAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTacaca                               >  1:2212793/1‑140 (MQ=33)
                                                                                      ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                                <  2:748102/140‑1 (MQ=2)
                                                                                      ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                               <  1:145435/140‑1 (MQ=2)
                                                                                      ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATAATCATGCCCCTTATGAc                                           >  2:130353/1‑128 (MQ=2)
                                                                                       tAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                               <  1:147754/140‑1 (MQ=2)
                                                                                        aaGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACg                             <  1:308122/140‑1 (MQ=2)
                                                                                          gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                           <  2:895913/140‑2 (MQ=2)
                                                                                           tCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGtac                          >  1:941206/1‑138 (MQ=2)
                                                                                              cGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTAc                       >  2:2282579/1‑140 (MQ=33)
                                                                                              cGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCATTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTa                                   >  2:1454463/1‑128 (MQ=2)
                                                                                              cGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATAATCATGCCCCTTATGACCTGGGc                                     <  1:874805/126‑1 (MQ=2)
                                                                                              cGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATAATCATGCCCCTTATGACCTGGGc                                     >  2:874805/1‑126 (MQ=2)
                                                                                               gCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACa                      <  1:1706359/140‑1 (MQ=32)
                                                                                                cAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACaa                     >  2:483219/1‑140 (MQ=2)
                                                                                                  aCGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATg                   <  2:348142/140‑1 (MQ=2)
                                                                                                   cGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGTGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATg                   >  1:346849/1‑139 (MQ=2)
                                                                                                   cGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATgg                  >  2:2191627/1‑140 (MQ=2)
                                                                                                    gAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc                 <  2:212025/140‑1 (MQ=2)
                                                                                                     aGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgcc                >  2:1102089/1‑139 (MQ=2)
                                                                                                     aGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgcc                >  2:298246/1‑139 (MQ=2)
                                                                                                     aGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgcc                >  2:2008989/1‑139 (MQ=2)
                                                                                                       cgcAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgccag              >  1:1550146/1‑137 (MQ=2)
                                                                                                       cgcAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                                >  2:702282/1‑122 (MQ=2)
                                                                                                                                         ||                                                                                                                     
AACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAG  >  NC_000913/224774‑225028

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