Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,381,705 G→A 6.2% P321P (CCG→CCA degQ → serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,381,7050GA 6.2% 201.4 / 2.8 64P321P (CCG→CCAdegQserine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base A (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

GGCACCGAGATGAGTGCCGATATCGCCAAAGCCTTCAACCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGG  >  NC_000913/3381559‑3381838
                                                                                                                                                  |                                                                                                                                     
ggCACCGAGATGAGTGCCGATATCGCCAAAGCCTTCAACCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGc                                                                                                                                      >  6:321735/1‑148 (MQ=255)
  cACCGAGATGAGTGCCGATATCGCCAAAGCCTTCAACCTTGACGTGCAGCGTGGCGCGTTTGTGAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTg                                                                                                                                    <  7:373774/148‑1 (MQ=255)
     cGAGATGAGTGCCGATATCGCCAAAGCCTTCAACCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGaa                                                                                                                                  >  5:287172/1‑147 (MQ=255)
      gagaTGAGTGCCGATATCGCCAAAGCCTTCAACCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATa                                                                                                                                <  7:318222/148‑1 (MQ=255)
         aTGAGTGCCGATATCGCCAAAGCCTTCAACCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCt                                                                                                                             <  1:376399/148‑1 (MQ=255)
          tGAGTGCCGATATCGCCAAAGCCTTCAACCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCtt                                                                                                                            <  4:197906/148‑1 (MQ=255)
            aGTGCCGATATCGCCAAAGCCTTCAACCTTGACGTGCATCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCATCGGGAAACCGCTGAATAGCTTTg                                                                                                                          <  5:368127/148‑1 (MQ=255)
            aGTGCCGATATCGCCAAAGCCTTCAACCTTGACGTGCAGCGTGGCGCGTTGGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTg                                                                                                                          <  4:105033/148‑1 (MQ=255)
                     atCGCCAAAGCCTTCAACCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTgc                                                                                                                 >  1:43996/1‑148 (MQ=255)
                        gCCAAAGCCTTCAACCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTgcgc                                                                                                               <  5:321735/147‑1 (MQ=255)
                          cAAAGCCTTCAACCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGctct                                                                                                            >  8:333978/1‑148 (MQ=255)
                            aaGCCTTCAACCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGTTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGa                                                                                                                      >  1:2924/1‑136 (MQ=255)
                            aaGCCTTCAACCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGa                                                                                                                      <  2:2924/136‑1 (MQ=255)
                              gCCTTCAACCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTg                                                                                                                       <  7:151545/133‑1 (MQ=255)
                              gCCTTCAACCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTg                                                                                                                       >  8:151545/1‑133 (MQ=255)
                                 ttCAACCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATcg                                                                                                     >  3:387523/1‑148 (MQ=255)
                                    aaCCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTTCCAGGTTCTGGCTCGGCAAATGCGGGCGTCAAAGCGGGGGATATTATTACGAGCCTCAAGGGCAAACCGCTGATTAGCTTTGCTGAGTTGCGCTCTCGTATcgcg                                                                                                   <  3:317318/147‑1 (MQ=255)
                                    aaCCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGa                                                                                                  <  7:333978/148‑1 (MQ=255)
                                      ccTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGacc                                                                                                >  2:360934/1‑148 (MQ=255)
                                       cTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGc                                                                                                                                      <  7:370003/109‑1 (MQ=255)
                                       cTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGc                                                                                                                                      >  8:370003/1‑109 (MQ=255)
                                       cTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGt                                                                                                                    <  1:10092/127‑1 (MQ=255)
                                       cTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGt                                                                                                                    >  2:10092/1‑127 (MQ=255)
                                                       gcgcGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGaa                                                                               <  4:387523/148‑1 (MQ=255)
                                                       gcgcGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGaa                                                                               <  7:189702/148‑1 (MQ=255)
                                                        cgcgTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGcc                                                                                        <  1:208532/138‑1 (MQ=255)
                                                        cgcgTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGcc                                                                                        >  2:208532/1‑138 (MQ=255)
                                                        cgcgTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGa                                                                                >  8:189702/1‑146 (MQ=255)
                                                          cgTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAATGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGt                                                                            <  2:43996/148‑1 (MQ=255)
                                                               gTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTTAAg                                                                        <  6:287172/147‑1 (MQ=255)
                                                                tCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCt                                                                      <  1:360934/148‑1 (MQ=255)
                                                                   gCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTgg                                                                   <  3:244484/148‑1 (MQ=255)
                                                                          gtTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCctgctg                                                            >  7:338192/1‑148 (MQ=255)
                                                                           ttGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCtgctgc                                                           >  5:37989/1‑148 (MQ=255)
                                                                             gCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCg                                                          >  5:373247/1‑147 (MQ=255)
                                                                                   ttCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACgg                                                    >  1:27098/1‑147 (MQ=255)
                                                                                    tCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAg                                                                                                                     >  5:367873/1‑81 (MQ=255)
                                                                                    tCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAg                                                                                                                     <  6:367873/81‑1 (MQ=255)
                                                                                    tCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGc                                                   >  5:209071/1‑147 (MQ=255)
                                                                                        gCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCaa                                                 <  6:373247/145‑1 (MQ=255)
                                                                                              cAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGa                                        <  6:209071/148‑1 (MQ=255)
                                                                                                 aaGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTa                                                        <  8:14155/129‑1 (MQ=255)
                                                                                                 aaGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTa                                                        >  7:14155/1‑129 (MQ=255)
                                                                                                  aGCGGGCGTCAAAGAGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGATGAGTTGCGATCGCGTATCGCGAACACTGAGCCGGGCGCGAAAGTGAGGCTTGGACTGCTGCGTAACTGCAAACAACTGCAAGTa                                    >  8:233632/1‑148 (MQ=255)
                                                                                                           cAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGaa                                                                               <  4:5055/96‑1 (MQ=255)
                                                                                                           cAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGaa                                                                               >  3:5055/1‑96 (MQ=255)
                                                                                                           cAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCACTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGc                                                                       <  1:14413/104‑1 (MQ=255)
                                                                                                           cAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCACTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGc                                                                       <  1:385511/104‑1 (MQ=255)
                                                                                                           cAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCACTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGc                                                                       >  2:14413/1‑104 (MQ=255)
                                                                                                           cAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCACTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGc                                                                       >  2:385511/1‑104 (MQ=255)
                                                                                                              aGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCctgc                                                              >  6:19779/1‑110 (MQ=255)
                                                                                                              aGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCctgc                                                              <  5:19779/110‑1 (MQ=255)
                                                                                                              aGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTc                        <  6:37989/148‑1 (MQ=255)
                                                                                                               gCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTa                                                                                                        >  1:237224/1‑67 (MQ=255)
                                                                                                               gCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTa                                                                                                        <  2:237224/67‑1 (MQ=255)
                                                                                                                cGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGaa                                       <  2:80991/131‑1 (MQ=255)
                                                                                                                cGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGaa                                       >  1:80991/1‑131 (MQ=255)
                                                                                                                      atatTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCTCGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCACACCACTGGAAGTAGAAGTGCCGCTCGATACCAg                >  3:358183/1‑148 (MQ=255)
                                                                                                                               ccAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTtc       >  5:421760/1‑148 (MQ=255)
                                                                                                                               ccAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACAAGCACCTCTtc       >  5:358802/1‑148 (MQ=255)
                                                                                                                                            cAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCgg  <  5:188259/140‑1 (MQ=255)
                                                                                                                                            cAAACAGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCgg  >  6:188259/1‑140 (MQ=255)
                                                                                                                                                ccGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCg                       <  8:73076/115‑1 (MQ=255)
                                                                                                                                                ccGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGAGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCg                       >  7:73076/1‑115 (MQ=255)
                                                                                                                                                  |                                                                                                                                     
GGCACCGAGATGAGTGCCGATATCGCCAAAGCCTTCAACCTTGACGTGCAGCGTGGCGCGTTTGTCAGCGAAGTGTTGCCAGGTTCTGGCTCGGCAAAAGCGGGCGTCAAAGCGGGCGATATTATTACCAGCCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGG  >  NC_000913/3381559‑3381838

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