Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 607,988 G→A *51* (TAG→TAA hokE → toxic polypeptide, small

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913607,9880GA100.0% 308.8 / NA 95*51* (TAG→TAAhokEtoxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (46/49);  total (46/49)

CCCTTGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT  >  NC_000913/607852‑608112
                                                                                                                                        |                                                                                                                            
cccTTGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAAcc                                                                                                                            >  1:260846/1‑139 (MQ=255)
  cTTGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGt                                                                                                                          <  2:1544731/139‑1 (MQ=255)
    tGCGGCAGTCATAGTGCTGTGTTTCACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATAGTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTtg                                                                                                                        >  2:755038/1‑139 (MQ=255)
    tGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTtg                                                                                                                        >  2:334653/1‑139 (MQ=255)
    tGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTtg                                                                                                                        >  2:781196/1‑139 (MQ=255)
        gCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                                    <  2:377424/139‑1 (MQ=255)
        gCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                                    <  2:1125466/139‑1 (MQ=255)
          aGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCggg                                                                                                                  <  2:66803/139‑1 (MQ=255)
          aGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCggg                                                                                                                  <  1:486774/139‑1 (MQ=255)
                aGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTaa                                                                                                            <  2:260846/139‑1 (MQ=255)
                       tgtTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATa                                                                                                     >  1:1012306/1‑139 (MQ=255)
                       tgtTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATa                                                                                                     >  2:1368421/1‑139 (MQ=255)
                       tgtTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATa                                                                                                     >  1:1445134/1‑139 (MQ=255)
                         tttAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAg                                                                                                   >  2:518375/1‑139 (MQ=255)
                         tttAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAg                                                                                                   >  2:660441/1‑139 (MQ=255)
                         tttAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAg                                                                                                   >  2:1418482/1‑139 (MQ=255)
                            aaCGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcgc                                                                                                >  1:1621710/1‑139 (MQ=255)
                             aCGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                               <  1:518375/139‑1 (MQ=255)
                             aCGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                               <  1:334653/139‑1 (MQ=255)
                               ggTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTaa                                                                                             >  1:1654760/1‑139 (MQ=255)
                                 tGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACt                                                                                           >  2:559026/1‑139 (MQ=255)
                                   cTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGt                                                                                                              >  1:185232/1‑118 (MQ=255)
                                   cTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGt                                                                                                              <  2:185232/118‑1 (MQ=255)
                                       gATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                     <  1:660441/139‑1 (MQ=255)
                                           tACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCAt                                                                                                      >  2:1247477/1‑118 (MQ=255)
                                           tACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCAt                                                                                                      <  1:1247477/118‑1 (MQ=255)
                                           tACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                     >  1:1061656/1‑135 (MQ=255)
                                           tACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                     <  2:1061656/135‑1 (MQ=255)
                                           tACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtgt                                                                                 <  2:1654760/139‑1 (MQ=255)
                                           tACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtgt                                                                                 >  2:556539/1‑139 (MQ=255)
                                           tACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtgt                                                                                 >  2:436476/1‑139 (MQ=255)
                                                 tCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                                                           <  2:1368443/139‑1 (MQ=255)
                                                    ggTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                                                        <  1:556539/139‑1 (MQ=255)
                                                     gTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                                                       <  2:1621710/139‑1 (MQ=255)
                                                          agaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                                                  <  2:1012306/139‑1 (MQ=255)
                                                          agaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                                                  <  1:559026/139‑1 (MQ=255)
                                                          agaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                                                  <  1:781196/139‑1 (MQ=255)
                                                                          tttACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGt                                                                                                              <  1:445934/79‑1 (MQ=255)
                                                                          tttACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGt                                                                                                              >  2:445934/1‑79 (MQ=255)
                                                                          tttACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcgc                                                                                                <  1:574253/93‑1 (MQ=255)
                                                                          tttACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcgc                                                                                                >  2:574253/1‑93 (MQ=255)
                                                                              cGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtg                                                                                  <  1:1227601/103‑1 (MQ=255)
                                                                              cGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtg                                                                                  >  2:1227601/1‑103 (MQ=255)
                                                                                 tGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt                                           <  2:1445134/139‑1 (MQ=255)
                                                                                  gAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa                                          >  1:1124882/1‑139 (MQ=255)
                                                                                   aaGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg                                               >  2:358407/1‑133 (MQ=255)
                                                                                   aaGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg                                               <  1:358407/133‑1 (MQ=255)
                                                                                    aGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                               >  2:539656/1‑84 (MQ=255)
                                                                                    aGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                               <  1:539656/84‑1 (MQ=255)
                                                                                    aGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTaa                                                                                        >  1:1237275/1‑91 (MQ=255)
                                                                                    aGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTaa                                                                                        <  2:1237275/91‑1 (MQ=255)
                                                                                         cGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAc                                   <  2:1124882/139‑1 (MQ=255)
                                                                                             atatTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                                                               <  1:1047255/107‑1 (MQ=255)
                                                                                             atatTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                                                               >  2:1047255/1‑107 (MQ=255)
                                                                                                 tGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                                    <  1:263869/50‑1 (MQ=255)
                                                                                                 tGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                                    >  2:263869/1‑50 (MQ=255)
                                                                                                 tGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGt                                                                                                              <  2:621167/56‑1 (MQ=255)
                                                                                                 tGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGt                                                                                                              >  1:621167/1‑56 (MQ=255)
                                                                                                 tGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtgt                                                                                 <  1:329699/85‑1 (MQ=255)
                                                                                                 tGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtgt                                                                                 >  2:329699/1‑85 (MQ=255)
                                                                                                 tGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTg                           >  2:167161/1‑139 (MQ=255)
                                                                                                   aGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCAt                                                                                                      <  1:704589/62‑1 (MQ=255)
                                                                                                   aGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCAt                                                                                                      >  2:704589/1‑62 (MQ=255)
                                                                                                    gTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGa                                                                                                                >  2:533685/1‑51 (MQ=255)
                                                                                                    gTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGa                                                                                                                <  1:533685/51‑1 (MQ=255)
                                                                                                     tttAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGtt                                                                                                                         <  1:1581731/41‑1 (MQ=255)
                                                                                                     tttAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGtt                                                                                                                         >  2:1581731/1‑41 (MQ=255)
                                                                                                      ttAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg                      <  1:333490/139‑1 (MQ=255)
                                                                                                        aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                        >  2:1309831/1‑119 (MQ=255)
                                                                                                        aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                        <  1:1309831/119‑1 (MQ=255)
                                                                                                             tGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCt               >  1:1230994/1‑139 (MQ=255)
                                                                                                              gTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                               >  2:283960/1‑58 (MQ=255)
                                                                                                              gTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                               <  1:283960/58‑1 (MQ=255)
                                                                                                              gTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTg              <  1:436476/139‑1 (MQ=255)
                                                                                                               ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcgc                                                                                                >  1:849760/1‑56 (MQ=255)
                                                                                                               ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcgc                                                                                                <  2:849760/56‑1 (MQ=255)
                                                                                                               ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACtt                                                                                          <  2:233879/62‑1 (MQ=255)
                                                                                                               ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACtt                                                                                          >  1:233879/1‑62 (MQ=255)
                                                                                                                tctcGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc            <  2:1230994/139‑1 (MQ=255)
                                                                                                                 ctcGCTTACGAACCGAAGAAGTAACCGTTGTGCggg                                                                                                                  >  2:991351/1‑36 (MQ=255)
                                                                                                                 ctcACTTACGAACCGAAGAAGTAACCGTTGTGCggg                                                                                                                  <  1:991351/36‑1 (MQ=255)
                                                                                                                    gCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc            >  1:824740/1‑135 (MQ=255)
                                                                                                                    gCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc            <  2:824740/135‑1 (MQ=255)
                                                                                                                       tACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATAtt     >  2:1060904/1‑139 (MQ=255)
                                                                                                                         cGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                                                     <  1:884726/89‑1 (MQ=255)
                                                                                                                         cGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                                                     >  2:884726/1‑89 (MQ=255)
                                                                                                                         cGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAGCTGGTCAGCAATTCTGGCCCATATTgg   <  1:1368421/139‑1 (MQ=255)
                                                                                                                         cGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTgg   >  2:815651/1‑139 (MQ=255)
                                                                                                                          gAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGt  <  1:1418482/139‑1 (MQ=255)
                                                                                                                              cGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                     <  2:56273/52‑1 (MQ=255)
                                                                                                                              cGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                     >  1:56273/1‑52 (MQ=255)
                                                                                                                                aagaagTAACCGTTGTGCGGGGAGTCATCCCATAAGCGCTATCTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc            <  2:1362107/123‑1 (MQ=255)
                                                                                                                                aagaagTAACCGTTGTGCGGGGAGTCATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc            >  1:1362107/1‑123 (MQ=255)
                                                                                                                                aagaagTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                               <  2:1001226/40‑1 (MQ=39)
                                                                                                                                aagaagTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                               >  1:1001226/1‑40 (MQ=255)
                                                                                                                                        |                                                                                                                            
CCCTTGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT  >  NC_000913/607852‑608112

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: