Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,378,286 C→A 10.6% L122L (CTG→CTT rplM ← 50S ribosomal subunit protein L13

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,378,2860CA10.6% 632.5 / 43.4 185L122L (CTG→CTTrplM50S ribosomal subunit protein L13
Reads supporting (aligned to +/- strand):  ref base C (87/78);  new base A (10/10);  total (97/88)
Fisher's exact test for biased strand distribution p-value = 8.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 189 total aligned reads displayed.

CCGGTTTGATGAAAACGCGAGCTGCGGAACTTTTGCGGCGACCAGTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATACCACCGATGTGGCCGGTGTGGTGATAGTACA  >  NC_000913/3378138‑3378434
                                                                                                                                                    |                                                                                                                                                    
ccGGTTTGATGAAAACGCGAGCTGCGGAACTTTTGCGGCGACCAGTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTc                                                                                                                                                      <  2:980366/149‑1 (MQ=255)
   gTTTGATGAAAACGCGAGCTGCGGAACTTTTGCGGCGACCAGTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGt                                                                                                                                                   >  3:727695/1‑149 (MQ=255)
          gaAAACGCGAGCTGCGGAACTTTTGCGGCGACCAGTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGa                                                                                                                                            >  2:895065/1‑149 (MQ=255)
          gaAAACGCGAGCTGCGGAACTTTTGCGGCGACCAGTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGa                                                                                                                                            >  2:953216/1‑149 (MQ=255)
          gaAAACGCGAGCTGCGGAACTTTTGCGGCGACCAGTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTACTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGa                                                                                                                                            >  2:1202706/1‑149 (MQ=255)
              aCGCGAGCTGCGGAACTTTTGCGGCGACCAGTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACAt                                                                                                                                        <  1:953216/149‑1 (MQ=255)
                    gCTGCGGAACTTTTGCGGCGACCAGTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACg                                                                                                                                  >  3:757731/1‑149 (MQ=255)
                       gCGGAACTTTTGCGGCGACCAGTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGAcc                                                                                                                               <  1:1202706/149‑1 (MQ=255)
                          gAACTTTTGCGGCGACCAGTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTGAAAGTTTACg                                                                                                                                              <  3:965691/131‑1 (MQ=255)
                            aCTTTTGCGGCGACCAGTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCg                                                                                                                          <  3:1172344/149‑1 (MQ=255)
                             cTTTTGCGGCGACCAGTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCgg                                                                                                                         <  1:895065/149‑1 (MQ=255)
                              ttttGCGGCGACCAGTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCa                                                                                                                                                     <  3:982549/120‑1 (MQ=255)
                                            gTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCtttt                                                                                                                  <  1:162178/141‑1 (MQ=255)
                                            gTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCtttt                                                                                                                  >  2:162178/1‑141 (MQ=255)
                                                  tagtaTTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATa                                                                                                                                       >  2:435699/1‑114 (MQ=255)
                                                      aTTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGaa                                                                                                                                           <  1:1430970/106‑1 (MQ=255)
                                                       ttGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCAc                                                                                                                                   <  2:150422/113‑1 (MQ=255)
                                                       ttGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCAc                                                                                                                                   >  1:150422/1‑113 (MQ=255)
                                                       ttGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACcg                                                                                               >  2:958313/1‑149 (MQ=255)
                                                         gATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCtt                                                                                                     >  1:279234/1‑141 (MQ=255)
                                                                    aTTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATcac                                                                                  >  4:27966/1‑149 (MQ=255)
                                                                          tataATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGa                                                                                                                                 >  1:505822/1‑96 (MQ=255)
                                                                          tataATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGa                                                                                                                                 <  2:505822/96‑1 (MQ=255)
                                                                          tataATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCtttt                                                                                                                  <  1:1010655/111‑1 (MQ=255)
                                                                          tataATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCttt                                                                                                                   <  4:1083310/110‑1 (MQ=255)
                                                                          tataATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCttt                                                                                                                   >  3:1083310/1‑110 (MQ=255)
                                                                          tataATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTg                                                                                                                 >  2:1338331/1‑112 (MQ=255)
                                                                          tataATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTg                                                                                                                 <  1:1338331/112‑1 (MQ=255)
                                                                               tCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTg                                                                                                                 >  3:808174/1‑107 (MQ=255)
                                                                               tCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTg                                                                                                                 <  4:808174/107‑1 (MQ=255)
                                                                               tCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCaa                                                                                      >  4:388785/1‑134 (MQ=255)
                                                                                cccGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACg                                                                      >  1:1052159/1‑149 (MQ=255)
                                                                                 ccGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGtt                                                                                                                                                  <  4:960514/72‑1 (MQ=255)
                                                                                 ccGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACcgcg                                                                                             >  2:1001499/1‑125 (MQ=255)
                                                                                 ccGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACcgcg                                                                                             <  1:1001499/125‑1 (MQ=255)
                                                                                    aTTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGa                                                                  >  2:313173/1‑149 (MQ=255)
                                                                                    aTTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGa                                                                  <  1:958313/149‑1 (MQ=255)
                                                                                       aGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCa                                                               <  4:161401/149‑1 (MQ=255)
                                                                                        gATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCaa                                                              >  4:758887/1‑149 (MQ=255)
                                                                                          tGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAatc                                                            <  2:1018945/149‑1 (MQ=255)
                                                                                            tCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCt                                                                                                                     <  1:621835/90‑1 (MQ=255)
                                                                                            tCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCt                                                                                                                     >  2:621835/1‑90 (MQ=255)
                                                                                            tCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGCCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAatcat                                                          >  1:251732/1‑149 (MQ=255)
                                                                                             cAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAatcatc                                                         >  3:601143/1‑149 (MQ=255)
                                                                                              aaGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCaa                                                                                                             <  1:1146884/96‑1 (MQ=255)
                                                                                              aaGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAatc                                                            <  2:376377/145‑1 (MQ=255)
                                                                                              aaGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCAtct                                                        >  2:96476/1‑149 (MQ=255)
                                                                                               aGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCaa                                                                                      >  4:743269/1‑118 (MQ=255)
                                                                                                 aaCTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCtt                                                     <  3:758887/149‑1 (MQ=255)
                                                                                                  aCTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGCTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTc                                                    <  4:965882/149‑1 (MQ=255)
                                                                                                   cTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCt                                                                                                                     >  1:348094/1‑83 (MQ=255)
                                                                                                    ttGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACa                                                                                                           <  3:1137333/92‑1 (MQ=255)
                                                                                                     tGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCaaa                                                 <  4:601143/149‑1 (MQ=255)
                                                                                                             cTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCtt                                         <  2:251732/149‑1 (MQ=255)
                                                                                                               gtgCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGt                                       <  2:1132181/149‑1 (MQ=255)
                                                                                                               gtgCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGt                                       <  3:27966/149‑1 (MQ=255)
                                                                                                                 gCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGttt                                     <  3:141047/149‑1 (MQ=255)
                                                                                                                    gcgTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCaa                                                                                      <  2:1128967/97‑1 (MQ=255)
                                                                                                                    gcgTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGAt                                  >  4:719091/1‑149 (MQ=255)
                                                                                                                    gcgTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGAt                                  <  4:1427991/149‑1 (MQ=255)
                                                                                                                     cgTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATa                                 <  2:1420518/149‑1 (MQ=255)
                                                                                                                      gTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAGCCGCGATTTCaa                                                                                      <  4:1163222/95‑1 (MQ=255)
                                                                                                                         gTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATaccac                             >  1:1407596/1‑149 (MQ=255)
                                                                                                                          ttGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGa                                                                  >  2:449886/1‑111 (MQ=255)
                                                                                                                            gtgCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCt                                                                             >  4:340955/1‑98 (MQ=255)
                                                                                                                             tgCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTg                                    <  3:20695/138‑1 (MQ=255)
                                                                                                                              gCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATACCACCGAtg                        >  3:275725/1‑149 (MQ=255)
                                                                                                                                tCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATACCACCGAtgtg                      >  4:545112/1‑149 (MQ=255)
                                                                                                                                       ttccgTAAACTGAAAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCa                                                   <  1:475042/110‑1 (MQ=255)
                                                                                                                                       ccGCGTAAACTTTCGGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATACCACCGATGTGGCCGgtg               >  4:329849/1‑149 (MQ=255)
                                                                                                                                       ccGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAat                                                             >  3:336257/1‑103 (MQ=255)
                                                                                                                                       ccGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGc                                           <  2:1240513/121‑1 (MQ=255)
                                                                                                                                       ccGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGc                                           >  1:1240513/1‑121 (MQ=255)
                                                                                                                                        cgcgTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGtt                                      >  3:434697/1‑125 (MQ=255)
                                                                                                                                         ccgTAAACTGAAAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCAtctc                                                       >  3:15766/2‑107 (MQ=255)
                                                                                                                                          cgTAAACTGAAAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGa                                                                  >  3:356927/1‑95 (MQ=255)
                                                                                                                                          cgTAAACTGAAAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATACCACCGAtg                        >  1:460256/1‑137 (MQ=255)
                                                                                                                                          cgTAAACTGAAAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATACCACCGAtg                        <  2:460256/137‑1 (MQ=255)
                                                                                                                                           gTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCaa                                                                                      >  4:61915/1‑74 (MQ=255)
                                                                                                                                           gTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCaa                                                                                      <  3:61915/74‑1 (MQ=255)
                                                                                                                                           gTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAg                                                                          <  4:1037727/86‑1 (MQ=255)
                                                                                                                                           gTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAg                                                                          >  3:1037727/1‑86 (MQ=255)
                                                                                                                                           gTAAACTGAAAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCt                                                                                                      <  4:949726/58‑1 (MQ=38)
                                                                                                                                           gTAAACTGAAAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCt                                                                                                      >  2:239327/1‑58 (MQ=38)
                                                                                                                                           gTAAACTGAAAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCt                                                                                                      <  1:239327/58‑1 (MQ=38)
                                                                                                                                           gTAAACTGAAAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCAtctc                                                       <  4:115866/105‑1 (MQ=255)
                                                                                                                                           gTAAACTGAAAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATACCACCGATGTGGCCGGTgtggt           >  2:1148335/1‑149 (MQ=255)
                                                                                                                                             aaaCTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATACCACCGATGTGGCCGGTGTGGTGa         >  3:851361/1‑149 (MQ=255)
                                                                                                                                              aaCTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGAt                                  >  1:692846/1‑123 (MQ=255)
                                                                                                                                              aaCTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGAt                                  <  2:692846/123‑1 (MQ=255)
                                                                                                                                               aCTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATACCACCGATGTGGCCGGTGTGGTGATa       >  4:1378875/1‑149 (MQ=255)
                                                                                                                                               aCTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATATCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGCTCGCTTGTTTGATACCACCGAGGTGGCCGGTGTGGTGATa       <  1:423774/149‑1 (MQ=255)
                                                                                                                                                cTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATACCACCGATGTGGCCGGTGTGGTGATAg      <  2:1407596/149‑1 (MQ=255)
                                                                                                                                                cTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATACCACCGATGTGGCCGGTGTGGTGATAg      <  2:744041/149‑1 (MQ=255)
                                                                                                                                                 tttCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAAt                                                                                     >  2:814280/1‑69 (MQ=255)
                                                                                                                                                 tttCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAAt                                                                                     <  1:814280/69‑1 (MQ=255)
                                                                                                                                                 tttCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGAt                                  >  3:959652/1‑120 (MQ=255)
                                                                                                                                                 tttCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATACCACCGATGTGGCCGGTGTGGTGATAGt     >  1:762299/1‑149 (MQ=255)
                                                                                                                                                  ttCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATACCACCGATGTGGCCGGTGTGGTGATAGTa    <  1:1130701/149‑1 (MQ=255)
                                                                                                                                                  ttCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATACCACCGATGTGGCCGGTGTGGGGATAGTa    >  1:440711/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCGGTTTGATGAAAACGCGAGCTGCGGAACTTTTGCGGCGACCAGTGCCGTAGTATTGATTTTCAGCCATTGCCTATAATCCCGATTAGATGTCAAGAACTTGCGGTTGCTGTGCCGCGTGGTTGTGCTCGTTACCCGCGTAAACTTTCAGTTTACGGAACATAGCACGACCCAGCGGGCCTTTTGGCAACATGCCTTTAACCGCGATTTCAATCACACGCTCAGGACGGCGAGCAATCATCTCTTCAAAGGTCGCTTGTTTGATACCACCGATGTGGCCGGTGTGGTGATAGTACA  >  NC_000913/3378138‑3378434

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