Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,921,384 G→A A388T (GCC→ACC)  rsmF → 16S rRNA m(5)C1407 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,921,3840GA98.5% 215.2 / NA 67A388T (GCC→ACC) rsmF16S rRNA m(5)C1407 methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (29/35);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACC  >  NC_000913/1921253‑1921516
                                                                                                                                   |                                                                                                                                    
cTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACg                                                                                                                              <  2:1460590/140‑1 (MQ=255)
cTACAGGTGTTGGCTGAAACTGAGATGAATACCTGCGTCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACg                                                                                                                              <  1:597407/140‑1 (MQ=255)
 tACAGGTGTTGGTTTAACCTGGGATGAAATCCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCACGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGAGATTAAACTTGCCGAGACGc                                                                                                                             <  2:207581/140‑1 (MQ=255)
      gtgtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCacaac                                                                                                                        <  1:360563/140‑1 (MQ=255)
       tggtGGCTTAATCTGGGATGAAAACCTTCGCCTTTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTCCCGAGACGCacaaca                                                                                                                       <  2:1358346/137‑1 (MQ=255)
        gtTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCAcaacaa                                                                                                                      <  1:1363349/140‑1 (MQ=255)
          tGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCGCGGTTGGCGATAAACCTTACCTAGACGCACAACAAAg                                                                                                                    >  1:685695/1‑140 (MQ=255)
               tAAACTGGGATGAAATCCTGCGTCTTTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTAt                                                                                                               <  1:403129/140‑1 (MQ=255)
               tAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTAt                                                                                                               <  1:206176/140‑1 (MQ=255)
                aaaCTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATc                                                                                                              >  2:1334845/1‑140 (MQ=255)
                aaaCTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATc                                                                                                              >  1:146463/1‑140 (MQ=255)
                      ggATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGgc                                                                                                        >  1:1094961/1‑140 (MQ=255)
                      ggATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGgc                                                                                                        >  1:157043/1‑140 (MQ=255)
                       gATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGgca                                                                                                       >  2:1505614/1‑140 (MQ=255)
                       gATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGgca                                                                                                       >  1:1533404/1‑140 (MQ=255)
                             aaCCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGa                                                                                                 <  2:157043/140‑1 (MQ=255)
                               ccTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAg                                                                                               <  1:1172880/140‑1 (MQ=255)
                                  gcgcCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCgg                                                                                            >  1:239182/1‑140 (MQ=255)
                                  gcgcCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCgg                                                                                            >  1:1138813/1‑140 (MQ=255)
                                  gcgcCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCgg                                                                                            >  1:808555/1‑140 (MQ=255)
                                     ccTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTa                                                                                         <  1:1524834/140‑1 (MQ=255)
                                     ccTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTa                                                                                         <  2:91412/140‑1 (MQ=255)
                                     ccTCTGGCAGCGTGACAAAGAACTGTGGTTGATCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTa                                                                                         <  2:548867/140‑1 (MQ=255)
                                     ccTCTGGCAGCGTGACAAAGAACTGAGGTGGTTCCCGTTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTa                                                                                         <  1:728921/140‑1 (MQ=255)
                                      ctctGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTAt                                                                                        <  1:779946/140‑1 (MQ=255)
                                        ctGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTTCCGAGACGAACATCAACAGTTATCGATGGCAACATGAAGGGGTTATTg                                                                                      >  2:728921/1‑140 (MQ=255)
                                          ggCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAg                                                                                               >  1:529180/1‑129 (MQ=255)
                                          ggCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAg                                                                                               <  2:529180/129‑1 (MQ=255)
                                               cGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAg                                                                                               <  1:1439198/124‑1 (MQ=255)
                                               cGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAg                                                                                               >  2:1439198/1‑124 (MQ=255)
                                                gTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGcc                                                                              <  1:1334845/140‑1 (MQ=255)
                                                gTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGcc                                                                              <  2:1533404/140‑1 (MQ=255)
                                                   aCAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCa                                                                           >  2:913010/1‑140 (MQ=255)
                                                    cAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCAc                                                                          <  2:1138813/140‑1 (MQ=255)
                                                    cAAAGAACCGTGGTTGTTCCCGTTGGGCATTCACGACCTGATCGGTAATGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCAc                                                                          <  2:685695/140‑1 (MQ=255)
                                                        gAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGa                                                                      <  1:913010/140‑1 (MQ=255)
                                                         aaCTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGat                                                                     <  2:1094961/140‑1 (MQ=255)
                                                         aaCTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGat                                                                     >  2:1114159/1‑140 (MQ=255)
                                                              tgGTTTTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCa                                                                           >  1:1097200/1‑129 (MQ=255)
                                                              tgGTTTTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCa                                                                           <  2:1097200/129‑1 (MQ=255)
                                                                 ttgttCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGaa                                                             >  2:353859/1‑140 (MQ=255)
                                                                       ccGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCAtt                                                       <  1:353859/140‑1 (MQ=255)
                                                                        cGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAAAGCATTc                                                      >  2:922328/1‑140 (MQ=255)
                                                                           tGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCAAGAAGCGGTTATTACTCATGCCTCACCCGATAATATGAATGCATTCGAg                                                   >  2:555306/1‑140 (MQ=255)
                                                                               cATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGa                                               >  2:1149744/1‑140 (MQ=255)
                                                                               cATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGa                                               >  2:163494/1‑140 (MQ=255)
                                                                               cATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTACTCATGCCTCACCCGACAATACGAATGCATTCGAGCTGa                                               >  2:529281/1‑140 (MQ=255)
                                                                                aTTGCAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTTCCGAGACACACAACAAAGGTTATCGCTGGCAGTATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTggc                                              <  2:783736/140‑3 (MQ=255)
                                                                                    aaGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCg                                          >  1:48983/1‑140 (MQ=255)
                                                                                        ccTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAgg                                      <  2:808555/140‑1 (MQ=255)
                                                                                        ccTGATCGGTAAAGTCCGATTTCCTCGGTTGGGGAATAATCTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAgg                                      <  2:605300/140‑1 (MQ=255)
                                                                                         cTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGa                                     >  1:322956/1‑140 (MQ=255)
                                                                                          tGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGaa                                    <  2:48983/140‑1 (MQ=255)
                                                                                           gATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAgg                                      >  2:1324738/1‑137 (MQ=255)
                                                                                           gATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAgg                                      <  1:1324738/137‑1 (MQ=255)
                                                                                           gATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAg                                   <  2:33826/140‑1 (MQ=255)
                                                                                            aTCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGc                                  >  2:778730/1‑140 (MQ=255)
                                                                                              cGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCgg                                <  2:146463/140‑1 (MQ=255)
                                                                                              cGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCACTTGCCTAACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCgg                                <  1:1411540/140‑1 (MQ=255)
                                                                                                gTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCggag                              <  2:1051889/140‑1 (MQ=255)
                                                                                                         gATTTTCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTAt                     >  2:959101/1‑140 (MQ=255)
                                                                                                           ttttCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGAATTTGAGTTGACACCGCAGGAGGCGGAGGAGTGGTATcg                   >  2:612888/1‑140 (MQ=255)
                                                                                                           ttttCTCGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCCTGCCTCACCCGATAGTATGAATGCATTTGAGTTGACACCGGAGGACGCGGAGGAGCGGTATcg                   >  2:573725/1‑140 (MQ=255)
                                                                                                                 cGGTTGGGGATTAAACTTACCGCGACGCACAACAAAGTTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGgcg             <  1:398003/140‑1 (MQ=255)
                                                                                                                 cGGTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGgcg             <  2:239182/140‑1 (MQ=255)
                                                                                                                  ggTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAAAGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGgcgc            <  1:922328/140‑1 (MQ=255)
                                                                                                                   gTTGGGGATTAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGcgcg           <  1:1505614/140‑1 (MQ=255)
                                                                                                                            tAAACTTACCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGGACCCAAGGAGTGCTATCGGGGGCGCGCTGTTTAcc  >  2:728459/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                                    
CTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAACTGTGGTTGTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACC  >  NC_000913/1921253‑1921516

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