Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,458,811 C→T A51V (GCG→GTG)  hofG → putative general protein secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,458,8110CT99.3% 1010.6 / NA 290A51V (GCG→GTG) hofGputative general protein secretion protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (147/140);  minor base G (2/0);  total (150/140)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 309 total aligned reads displayed.

TTATCGGCGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAGCGCCGA  >  NC_000913/3458721‑3458900
                                                                                          |                                                                                         
ttATCGGCGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGgtgc                                                                                         >  2:8704484/1‑93 (MQ=255)
ttATCGGCGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGgtgc                                                                                         <  1:572315/93‑1 (MQ=255)
ttATCGGCGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGgtgc                                                                                         <  2:1099054/93‑1 (MQ=255)
   tCGGCGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTgg                                                                                      <  1:5508075/93‑1 (MQ=255)
      gCGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGaaa                                                                                   >  1:2477841/1‑93 (MQ=255)
       cGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGaaaa                                                                                  <  2:1932001/93‑1 (MQ=255)
        gTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAAc                                                                                 >  1:1456383/1‑93 (MQ=255)
        gTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAAc                                                                                 >  2:6223324/1‑93 (MQ=255)
        gTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAAc                                                                                 >  1:909046/1‑93 (MQ=255)
        gTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAAc                                                                                 >  1:3211270/1‑93 (MQ=255)
         tatttGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACg                                                                                >  2:7825034/4‑93 (MQ=255)
         tACTTGCCAGTCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACg                                                                                >  1:1218076/1‑93 (MQ=255)
         tACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACg                                                                                >  2:6105156/1‑93 (MQ=255)
         tACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACg                                                                                <  1:926850/93‑1 (MQ=255)
          aCTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGc                                                                               >  2:1711352/1‑93 (MQ=255)
            ttGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGccc                                                                             <  1:5365632/93‑1 (MQ=255)
            ttGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGccc                                                                             <  1:3111219/93‑1 (MQ=255)
             tGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCt                                                                            >  1:6165218/1‑93 (MQ=255)
               ccAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTg                                                                          <  2:6753421/93‑1 (MQ=255)
                 aGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGat                                                                        <  2:4145229/93‑1 (MQ=255)
                 aGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGat                                                                        <  2:5251384/93‑1 (MQ=255)
                 aGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGat                                                                        >  2:6514540/1‑93 (MQ=255)
                 aGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGat                                                                        >  2:1262190/1‑93 (MQ=255)
                  gCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGata                                                                       >  1:1611176/1‑93 (MQ=255)
                  gCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGata                                                                       >  2:6470323/1‑93 (MQ=255)
                  gCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGata                                                                       >  2:264463/1‑93 (MQ=255)
                   ccTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGatat                                                                      <  1:5856238/93‑1 (MQ=255)
                    cTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATg                                                                     >  2:8798345/1‑93 (MQ=255)
                     tggtggTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGt                                                                    <  1:3597528/93‑1 (MQ=255)
                      ggtggtTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTa                                                                   >  2:3506765/1‑93 (MQ=255)
                       gtggtTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTAc                                                                  >  2:762283/1‑93 (MQ=255)
                        tggtTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACa                                                                 <  2:2106661/93‑1 (MQ=255)
                        tggtTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACa                                                                 <  2:3689836/93‑1 (MQ=255)
                           ttCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAAc                                                              >  2:3866069/1‑93 (MQ=255)
                            tCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACt                                                             >  2:5510149/1‑93 (MQ=255)
                            tCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACt                                                             >  1:1503137/1‑93 (MQ=255)
                             cctaacctaaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTc                                                            <  2:593989/93‑1 (MQ=255)
                             cctaacctaaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTc                                                            >  2:188896/1‑93 (MQ=255)
                               taacctaaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGa                                                          <  2:6343930/93‑1 (MQ=255)
                                 acctaaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACa                                                        >  2:3350683/1‑93 (MQ=255)
                                 acctaaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACa                                                        >  2:6924773/1‑93 (MQ=255)
                                 acctaaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACa                                                        >  2:4956180/1‑93 (MQ=255)
                                  cctaaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACaa                                                       <  1:3900126/93‑1 (MQ=255)
                                   ctaaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAat                                                      >  2:311109/1‑93 (MQ=255)
                                   ctaaTGGACAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAat                                                      <  2:4901069/93‑1 (MQ=255)
                                     aaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAatca                                                    <  2:5314608/93‑1 (MQ=255)
                                      aTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAatcat                                                   <  1:8557760/93‑1 (MQ=255)
                                      aTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAatcat                                                   <  2:7685013/93‑1 (MQ=255)
                                       tGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAatcatc                                                  <  1:696693/93‑1 (MQ=255)
                                        gggCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAAtcatca                                                 >  2:5388121/1‑93 (MQ=255)
                                        gggCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAAtcatca                                                 <  1:1040911/93‑1 (MQ=255)
                                         ggCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCAc                                                <  2:4664038/93‑1 (MQ=255)
                                         ggCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCAc                                                >  1:6763911/1‑93 (MQ=255)
                                          gCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACt                                               >  1:7774180/1‑93 (MQ=255)
                                          gCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCGTGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACt                                               <  2:7519901/93‑1 (MQ=255)
                                           cAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTa                                              <  1:5268539/93‑1 (MQ=255)
                                            aaTAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTAt                                             >  2:605880/1‑93 (MQ=255)
                                             aTAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATc                                            >  1:407959/1‑93 (MQ=255)
                                              tAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATcc                                           <  2:626886/93‑1 (MQ=255)
                                               aaaGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACACACTCGACAATCATCACTATCCa                                          >  2:1169240/1‑93 (MQ=255)
                                                 aGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAAc                                        <  2:5518821/93‑1 (MQ=255)
                                                 aGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAAc                                        <  1:3638770/93‑1 (MQ=255)
                                                  gAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATCTTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAAcc                                       >  1:1816520/1‑93 (MQ=255)
                                                  gAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAAcc                                       >  1:1729925/1‑93 (MQ=255)
                                                   aaaaGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACca                                      >  2:6635894/1‑93 (MQ=255)
                                                    aaaGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACcac                                     >  1:3845680/1‑93 (MQ=255)
                                                       ggCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACaaa                                  >  2:7161584/1‑93 (MQ=255)
                                                       ggCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACaaa                                  >  1:1956431/1‑93 (MQ=255)
                                                        gCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAAt                                 >  1:5332662/1‑93 (MQ=255)
                                                        gCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAAt                                 <  2:8565849/93‑1 (MQ=255)
                                                            aTAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAgg                             >  2:6670359/1‑93 (MQ=255)
                                                            aTAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAgg                             <  2:3608229/93‑1 (MQ=255)
                                                            aTAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAgg                             <  2:8118801/93‑1 (MQ=255)
                                                             tAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAggg                            >  1:8052296/1‑93 (MQ=255)
                                                              aaGTAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAgggg                           <  2:1526240/93‑1 (MQ=255)
                                                               aGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACTACAAATCAgtgga                          >  2:4095534/1‑92 (MQ=255)
                                                               aGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGc                          >  1:8411357/1‑93 (MQ=255)
                                                               aGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGc                          >  2:2477000/1‑93 (MQ=255)
                                                                gCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCt                         <  1:5670575/93‑1 (MQ=255)
                                                                gCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCt                         <  1:6270970/93‑1 (MQ=255)
                                                                gCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCt                         >  1:4751175/1‑93 (MQ=255)
                                                                 cAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCtt                        <  1:593890/93‑1 (MQ=255)
                                                                 cAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCtt                        <  2:103922/93‑1 (MQ=255)
                                                                    aaaGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCATTATCCAACCACAAATCAGGGGCTTGaa                     >  2:5664973/1‑93 (MQ=255)
                                                                    aaaGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGaa                     <  1:727793/93‑1 (MQ=255)
                                                                     aaGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAAt                    <  1:4858551/93‑1 (MQ=255)
                                                                      aGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATc                   >  1:6818997/1‑93 (MQ=255)
                                                                       gCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATcc                  <  1:2807628/93‑1 (MQ=255)
                                                                          gTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTta               >  1:4208310/1‑93 (MQ=255)
                                                                          gTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTta               <  1:2125859/93‑1 (MQ=255)
                                                                          ctCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTta               <  1:6469669/92‑1 (MQ=255)
                                                                               cGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGa          <  1:3435026/93‑1 (MQ=255)
                                                                               cGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGa          <  1:4416111/93‑1 (MQ=255)
                                                                               cGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGa          <  1:6030411/93‑1 (MQ=255)
                                                                                gATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGaa         <  1:1471540/93‑1 (MQ=255)
                                                                                  tatTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAgc       <  2:7319715/93‑1 (MQ=255)
                                                                                  tatTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAgc       <  2:4252706/93‑1 (MQ=255)
                                                                                    ttGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAgcgc     >  1:3332064/1‑93 (MQ=255)
                                                                                     tgtgGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAGCGcc    <  2:5080411/93‑1 (MQ=255)
                                                                                      gtggtgCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAGCGCCg   <  2:27384/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTATCGGCGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAGCGCCGA  >  NC_000913/3458721‑3458900

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