Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,943,843 C→A 13.9% noncoding (140/2904 nt) rrlC → 23S ribosomal RNA of rrnC operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,943,8430CA13.9% 199.9 / 20.6 65noncoding (140/2904 nt)rrlC23S ribosomal RNA of rrnC operon
Reads supporting (aligned to +/- strand):  ref base C (38/18);  new base A (3/6);  total (41/24)
Fisher's exact test for biased strand distribution p-value = 6.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

GAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGT  >  NC_000913/3943702‑3943987
                                                                                                                                             |                                                                                                                                                
attcaggctctgggctgctccccgttcgctcgccgctactgggggaatctcggttgatttcttttcctcggggtacttagatgtttcagttcccccggttcgcctcattaacctatggattcagttaatgatAGTGTGTTTCGacaca                                                                                                                                            >  2:396528‑M1/133‑148 (MQ=255)
  tcaggctctgggctgctccccgttcgctcgccgctactgggggaatctcggttgatttcttttcctcggggtacttagatgtttcagttcccccggttcgcctcattaacctatggattcagttaatgatAGTGTGTTTCGACACACt                                                                                                                                          >  1:639783‑M1/131‑148 (MQ=255)
     ggctctgggctgctccccgttcgctcgccgctactgggggaatctcggttgatttcttttcctcggggtacttagatgtttcagttcccccggttcgcctcattaacctatggattcagttaatgatAGTGTGTTTCGACACACTATCa                                                                                                                                      >  3:821905‑M1/128‑149 (MQ=255)
                                                    gttgatttcttttcctcggggtacttagatgtttcagttcccccggttcgcctcattaacctatggattcagttaatgatAGTGTGTTTCGacaca                                                                                                                                            <  1:393235‑M1/16‑1 (MQ=255)
                                                    gttgatttcttttcctcggggtacttagatgtttcagttcccccggttcgcctcattaacctatggattcagttaatgatAGTGTGTTTCGacaca                                                                                                                                            >  2:393235‑M1/81‑96 (MQ=255)
                                                                                                        tcattaacctatggattcagttaatgatAGTGTGTCGAAACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAgcggcg                                    <  1:778808‑M1/120‑1 (MQ=255)
                                                                                                                   tggattcagttaatgatAGTGTGTCGAAACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACgggg                         >  2:778808‑M1/18‑148 (MQ=255)
                                                                                                                       ttcagttaatgatAGTGTGTCGAAACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGgcgagcga                               <  2:417189‑M1/125‑1 (MQ=255)
                                                                                                                       ttcagttaatgatAGTGTGTCGAAACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGgcgagcga                               >  1:417189‑M1/14‑138 (MQ=255)
                                                                                                                        tcagttaatgatAGTGTGTCGAAACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGaa                                                                <  1:565331‑M1/92‑1 (MQ=255)
                                                                                                                        tcagttaatgatAGTGTGTCGAAACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGaa                                                                >  2:565331‑M1/13‑104 (MQ=255)
                                                                                                                            ttaatgatAGTGTGTCGAAACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCag               <  4:637431‑M1/141‑1 (MQ=255)
                                                                                                                               atgatAGTGTGTCGAAACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGc            <  3:205454‑M1/144‑1 (MQ=255)
                                                                                                                               atgatAGTGTGTCGAAACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGc            <  4:826076‑M1/144‑1 (MQ=255)
                                                                                                                                    AGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGATGCGAACCGAGGGAACTGTAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGaa       >  4:775683‑M1/1‑149 (MQ=255)
                                                                                                                                    AGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGaa       <  3:393838‑M1/149‑1 (MQ=255)
                                                                                                                                    AGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGaa       >  3:893177‑M1/1‑149 (MQ=255)
                                                                                                                                    AGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGaa       <  4:519283‑M1/149‑1 (MQ=255)
                                                                                                                                    AGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGaa       >  4:792168‑M1/1‑149 (MQ=255)
                                                                                                                                    AGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGaa       <  4:912198‑M1/149‑1 (MQ=255)
                                                                                                                                    AGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGa        <  3:903986‑M1/148‑1 (MQ=255)
                                                                                                                                    AGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGa        >  2:232465‑M1/1‑148 (MQ=255)
                                                                                                                                    AGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGa        <  1:289002‑M1/148‑1 (MQ=255)
                                                                                                                                    AGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGa        >  4:903986‑M1/1‑148 (MQ=255)
                                                                                                                                    AGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGa        <  2:64624‑M1/148‑1 (MQ=255)
                                                                                                                                     GTGTGTTTCGACACTCTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGaa       <  1:590804‑M1/148‑1 (MQ=255)
                                                                                                                                     GTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGaa       >  2:90812‑M1/1‑148 (MQ=255)
                                                                                                                                     GTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAAt      <  3:342629‑M1/149‑1 (MQ=255)
                                                                                                                                     GTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAAt      >  4:458694‑M1/1‑149 (MQ=255)
                                                                                                                                      TGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTc                                                  <  3:711772‑M1/104‑1 (MQ=255)
                                                                                                                                      TGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTc                                                  >  4:711772‑M1/1‑104 (MQ=255)
                                                                                                                                      TGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAAt      <  1:476105‑M1/148‑1 (MQ=255)
                                                                                                                                      TGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAAt      >  2:947600‑M1/1‑148 (MQ=255)
                                                                                                                                      TGTGTCTCGACACACTATCATTAACTGAATCCATAGGGTAATGAGGCGAACCGGGGGAACTGAAACCTCTAAGTACCCCGAGGCAAAGAAATCAACCGAGATTCCCACAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGGCTGAAt      <  2:908529‑M1/148‑1 (MQ=255)
                                                                                                                                       GTGTTTCGACACACTATCATTAACTGAATTCATAGTGTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATc     >  1:553523‑M1/1‑148 (MQ=255)
                                                                                                                                       GTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGaa                                                                                          >  1:885837‑M1/1‑63 (MQ=255)
                                                                                                                                       GTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGaa                                                                                          <  2:885837‑M1/63‑1 (MQ=255)
                                                                                                                                       GTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATc     >  1:287846‑M1/1‑148 (MQ=255)
                                                                                                                                       GTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATc     <  1:75259‑M1/148‑1 (MQ=255)
                                                                                                                                       GTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATc     <  2:64906‑M1/148‑1 (MQ=255)
                                                                                                                                       GTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATc     <  1:232683‑M1/148‑1 (MQ=255)
                                                                                                                                       GTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCa    >  3:57899‑M1/1‑149 (MQ=255)
                                                                                                                                       GTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCa    >  4:816682‑M1/1‑149 (MQ=255)
                                                                                                                                       GTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCa    >  4:469198‑M1/1‑149 (MQ=255)
                                                                                                                                       GTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCa    >  4:760830‑M1/1‑149 (MQ=255)
                                                                                                                                       GTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCa    >  3:276821‑M1/1‑149 (MQ=255)
                                                                                                                                        TGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAATGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCa    >  1:736010‑M1/1‑148 (MQ=255)
                                                                                                                                        TGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCa    >  2:655841‑M1/1‑148 (MQ=255)
                                                                                                                                        TGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAg   >  3:548738‑M1/1‑149 (MQ=255)
                                                                                                                                        TGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAg   <  3:354006‑M1/149‑1 (MQ=255)
                                                                                                                                        TGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAg   >  4:807745‑M1/1‑149 (MQ=255)
                                                                                                                                         GTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAATCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATGAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGt  >  4:803353‑M1/1‑149 (MQ=255)
                                                                                                                                         GTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGATAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGt  >  4:68623‑M1/1‑149 (MQ=255)
                                                                                                                                         GTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAg   >  1:464245‑M1/1‑148 (MQ=255)
                                                                                                                                         GTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAg   >  1:235367‑M1/1‑148 (MQ=255)
                                                                                                                                         GTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAg   >  1:789045‑M1/1‑148 (MQ=255)
                                                                                                                                         GTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAg   >  2:114675‑M1/1‑148 (MQ=255)
                                                                                                                                         GTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAg   >  2:134952‑M1/1‑148 (MQ=255)
                                                                                                                                         GTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAg   >  2:371991‑M1/1‑148 (MQ=255)
                                                                                                                                         GTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAg   >  2:528977‑M1/1‑148 (MQ=255)
                                                                                                                                         GTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGt  >  4:354006‑M1/1‑149 (MQ=255)
                                                                                                                                         GTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGt  >  3:762133‑M1/1‑149 (MQ=255)
                                                                                                                                         GTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCGAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGt  >  4:516187‑M1/1‑149 (MQ=255)
                                                                                                                                          TTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGgcgag                                  >  4:287978‑M1/1‑116 (MQ=255)
                                                                                                                                          TTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGgcgag                                  <  3:287978‑M1/116‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
GAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGT  >  NC_000913/3943702‑3943987

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