Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 224,915 A→C noncoding (1145/1542 nt) rrsH → 16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913224,9150AC90.1% 54.3 / 74.7 302noncoding (1145/1542 nt)rrsH16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (17/13);  new base C (107/165);  total (124/178)
Fisher's exact test for biased strand distribution p-value = 7.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

Only 100 of 308 total aligned reads displayed.

TTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATA  >  NC_000913/224778‑225018
                                                                                                                                          |                                                                                                       
           gATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGAc                                                                                            >  1:1588311/1‑140 (MQ=40)
            aTTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACt                                                                                           >  1:1752143/1‑140 (MQ=39)
                      tCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATa                                                                                 <  2:1737158/140‑1 (MQ=34)
                              ttcaagaCAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCggag                                                                         >  1:1142337/5‑140 (MQ=2)
                                tgAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTggagga                                                                       <  1:1704560/140‑1 (MQ=32)
                                tgAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTggagga                                                                       >  1:1410256/1‑140 (MQ=32)
                                   gaCAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAgg                                                                    <  2:616344/140‑1 (MQ=2)
                                        gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAAc                                                                              >  2:66936/5‑125 (MQ=2)
                                        gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  2:951339/5‑140 (MQ=2)
                                        gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  1:1211099/5‑140 (MQ=2)
                                        gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCggagg                                                                        >  2:87623/5‑131 (MQ=2)
                                        gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  1:4573/5‑140 (MQ=2)
                                        gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  1:1184577/5‑140 (MQ=2)
                                        gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  1:199144/5‑140 (MQ=2)
                                        gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTCGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  1:607671/5‑140 (MQ=2)
                                        gtggtgCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCggagga                                                                       >  1:864143/1‑133 (MQ=2)
                                        gtggtgCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  2:765329/1‑140 (MQ=2)
                                        gtggtgCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  1:129289/1‑140 (MQ=2)
                                        gtggtgCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  1:178076/1‑140 (MQ=2)
                                        gtggtgCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               <  1:283631/140‑1 (MQ=2)
                                        gTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTgggg                                                               >  2:1497198/1‑140 (MQ=32)
                                         tggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  1:832369/4‑140 (MQ=2)
                                         tggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGa                                                              <  2:934162/137‑1 (MQ=2)
                                          ggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGAt                                                             <  1:806140/138‑1 (MQ=2)
                                          ggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGAt                                                             <  1:110111/138‑1 (MQ=2)
                                           gtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGAt                                                             <  1:516475/139‑1 (MQ=2)
                                           gtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATg                                                            >  1:960082/2‑140 (MQ=2)
                                             gcATGGTTGTCGTCAGCTCGTGTCGTAAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGAc                                                          <  2:219871/140‑1 (MQ=2)
                                              cATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACg                                                         <  2:570152/140‑1 (MQ=2)
                                                  gttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaa                                                     <  2:129289/138‑1 (MQ=2)
                                                      tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGa                                                              >  1:850552/1‑128 (MQ=2)
                                                      tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtc                                                  >  1:902450/1‑140 (MQ=2)
                                                      tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtca                                                 >  1:644917/1‑140 (MQ=2)
                                                      tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtca                                                 >  2:806782/1‑140 (MQ=2)
                                                      tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGa                                                              >  2:477411/1‑127 (MQ=2)
                                                      tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtca                                                 >  1:1166517/1‑140 (MQ=2)
                                                      tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtca                                                 >  1:1231169/1‑140 (MQ=2)
                                                         tcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtcatca                                              <  1:1012359/140‑1 (MQ=2)
                                                           aGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATcatcat                                             >  1:583785/1‑140 (MQ=2)
                                                           aGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATg                                            <  1:1219998/140‑1 (MQ=2)
                                                           aGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATg                                            >  1:759649/1‑140 (MQ=11)
                                                            gCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                           <  2:318631/140‑2 (MQ=2)
                                                            gCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                           <  1:695428/140‑2 (MQ=2)
                                                            gCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                           >  2:1131151/1‑139 (MQ=2)
                                                            gCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                           <  2:178076/140‑2 (MQ=2)
                                                             cTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                           <  2:1114696/140‑2 (MQ=2)
                                                             cTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                           >  1:78243/1‑138 (MQ=2)
                                                                gtgtCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCt                                       >  1:386029/1‑140 (MQ=2)
                                                                gtgtCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCt                                       >  2:404948/1‑140 (MQ=2)
                                                                 tgtcgtGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                      <  1:476361/140‑1 (MQ=2)
                                                                 tgtcgtGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                      <  1:34586/140‑1 (MQ=2)
                                                                 tgtcgtGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                      <  2:657762/140‑1 (MQ=2)
                                                                 tgtcgtGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                      <  2:1149083/140‑1 (MQ=2)
                                                                 tgtcgtGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                      <  1:537028/140‑1 (MQ=2)
                                                                 tgtcgtGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                      <  2:1677397/140‑1 (MQ=2)
                                                                  gtcgtgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTa                                     >  1:665025/4‑140 (MQ=2)
                                                                     gtgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaaat                                                   <  1:257526/123‑3 (MQ=2)
                                                                      tgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTGTTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAc                                 <  2:1173602/140‑1 (MQ=2)
                                                                      tgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAc                                 <  1:1036868/140‑1 (MQ=2)
                                                                      tgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAc                                 <  1:1201321/140‑1 (MQ=2)
                                                                      tgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAc                                 <  2:737499/140‑1 (MQ=2)
                                                                       gagaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAc                                 <  2:710837/137‑1 (MQ=2)
                                                                       gagaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                           >  2:1024654/4‑128 (MQ=2)
                                                                       gagaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcc                                <  1:15439/137‑1 (MQ=2)
                                                                        agaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcct                               <  2:458/138‑2 (MQ=2)
                                                                        agaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcct                               <  2:686980/138‑2 (MQ=2)
                                                                        agaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAGTCATCATGCCCCTTATGAcct                               <  1:410015/138‑2 (MQ=2)
                                                                        agaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcct                               <  2:24334/138‑2 (MQ=2)
                                                                         gaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcctg                              <  2:553322/139‑3 (MQ=2)
                                                                         gaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcctg                              <  2:533615/139‑3 (MQ=2)
                                                                             ttGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAAACGGCGGGAGGTGGGGAAGACGTCAAATCATCATGCCCCTTATGACCTGGGCt                          >  2:592137/1‑140 (MQ=2)
                                                                              tGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                      <  1:970694/128‑1 (MQ=2)
                                                                              tGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTa                         <  2:1193238/140‑1 (MQ=2)
                                                                               gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  2:1231169/140‑1 (MQ=2)
                                                                                ggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  2:832369/140‑1 (MQ=2)
                                                                                 gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCTGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  1:502462/140‑1 (MQ=2)
                                                                                 gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  2:14476/140‑1 (MQ=2)
                                                                                  ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                     <  2:644917/140‑1 (MQ=2)
                                                                                  ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                     <  2:778466/140‑1 (MQ=2)
                                                                                  ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCGAATCATCATGCCCCTTATGACCTGGGCTacaca                     >  1:430534/1‑140 (MQ=2)
                                                                                  ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                     <  2:515424/140‑1 (MQ=2)
                                                                                   tAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                     <  2:1004608/140‑1 (MQ=2)
                                                                                    aaGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcc                                >  2:139174/1‑127 (MQ=2)
                                                                                     aGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGt                  >  2:1771625/1‑140 (MQ=32)
                                                                                     aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  2:356619/135‑1 (MQ=2)
                                                                                      gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 <  2:910073/140‑2 (MQ=2)
                                                                                      gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 <  1:729707/140‑2 (MQ=2)
                                                                                      gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 >  2:277854/1‑139 (MQ=2)
                                                                                       tCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAc                                 <  1:1076668/123‑1 (MQ=2)
                                                                                          cGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  1:1001282/131‑1 (MQ=2)
                                                                                          cGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGtactac             >  2:372530/1‑140 (MQ=2)
                                                                                          cGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGtactac             >  2:106839/1‑140 (MQ=2)
                                                                                           gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  2:1211099/140‑1 (MQ=2)
                                                                                           gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  2:1078955/140‑1 (MQ=2)
                                                                                            cAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACaa           >  2:394804/1‑140 (MQ=2)
                                                                                               cGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATgg        <  1:1497198/140‑1 (MQ=32)
                                                                                                gAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:869305/140‑1 (MQ=2)
                                                                                                gAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:977036/140‑1 (MQ=2)
                                                                                                gAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:466811/140‑1 (MQ=2)
                                                                                                gAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  1:1131151/140‑1 (MQ=2)
                                                                                                                                          |                                                                                                       
TTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATA  >  NC_000913/224778‑225018

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