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% 270.1 / NA 82*51* (TAG→TAAhokEtoxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/46);  total (36/46)

CTTGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT  >  NC_000913/607854‑608112
                                                                                                                                      |                                                                                                                            
cTTGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGt                                                                                                                          <  2:1086103/139‑1 (MQ=255)
cTTGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGt                                                                                                                          <  2:597320/139‑1 (MQ=255)
   gCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTtgt                                                                                                                       >  2:1533356/1‑139 (MQ=255)
   gCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTtgt                                                                                                                       >  2:136929/1‑139 (MQ=255)
   gCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTtgt                                                                                                                       >  2:97942/1‑139 (MQ=255)
      gCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                                    <  1:1681656/139‑1 (MQ=255)
        aGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGATTTTACGGTGAAGGAACGTAATATTGAGTTTACGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCggg                                                                                                                  >  1:1222504/1‑139 (MQ=255)
        aGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCggg                                                                                                                  <  1:97942/139‑1 (MQ=255)
        aGTCATAGTGCTGGGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCggg                                                                                                                  <  1:50854/139‑1 (MQ=255)
          tCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGa                                                                                                                <  1:1301312/139‑1 (MQ=255)
                         tAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcg                                                                                                 <  2:1560728/139‑1 (MQ=255)
                          aaCGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcgc                                                                                                <  1:1174661/139‑1 (MQ=255)
                          aaCGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcgc                                                                                                <  1:136929/139‑1 (MQ=255)
                           aCGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                               <  2:44789/139‑1 (MQ=255)
                                 cTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGt                                                                                                                          <  1:537429/106‑1 (MQ=255)
                                 cTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGt                                                                                                                          >  2:537429/1‑106 (MQ=255)
                                     gATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                     <  1:1533356/139‑1 (MQ=255)
                                         tACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTaa                                                                                        <  1:677667/132‑1 (MQ=255)
                                         tACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTaa                                                                                        >  2:677667/1‑132 (MQ=255)
                                                cTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGt                                                                                                                          >  2:947827/1‑91 (MQ=255)
                                                cTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGt                                                                                                                          <  1:947827/91‑1 (MQ=255)
                                                   gTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                                                       <  1:896402/139‑1 (MQ=255)
                                                                  tgtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                                    <  2:1022374/79‑1 (MQ=255)
                                                                  tgtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                                    >  1:1022374/1‑79 (MQ=255)
                                                                  tgtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtgt                                                                                 <  2:1517531/114‑1 (MQ=255)
                                                                  tgtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtgt                                                                                 >  1:1517531/1‑114 (MQ=255)
                                                                   gtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                                                       <  2:1652624/139‑1 (MQ=255)
                                                                          tACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGa                                                                                                                <  1:580144/75‑1 (MQ=255)
                                                                          tACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGa                                                                                                                >  2:580144/1‑75 (MQ=255)
                                                                          tACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                >  2:1285270/1‑139 (MQ=255)
                                                                            cGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTa                                                                                                             <  1:130457/76‑1 (MQ=255)
                                                                            cGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTa                                                                                                             >  2:130457/1‑76 (MQ=255)
                                                                             ggTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                                                        >  2:1696999/1‑112 (MQ=255)
                                                                             ggTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                                                        <  1:1696999/112‑1 (MQ=255)
                                                                              gTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt                                            <  2:1684879/139‑1 (MQ=255)
                                                                                 aaGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa                                          <  2:1547550/138‑1 (MQ=255)
                                                                                 aaGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa                                          >  1:1547550/1‑138 (MQ=255)
                                                                                  aGGAACGTAATCTTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                        <  2:1222504/139‑1 (MQ=255)
                                                                                  aGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCgggg                                                                                                                 <  1:1080796/66‑1 (MQ=255)
                                                                                  aGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCgggg                                                                                                                 >  2:1080796/1‑66 (MQ=255)
                                                                                  aGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                        <  1:1485892/139‑1 (MQ=255)
                                                                                  aGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCAAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                        >  1:143686/1‑139 (MQ=255)
                                                                                        gTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAAt                                                                                                           >  2:1032691/1‑66 (MQ=255)
                                                                                        gTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAAt                                                                                                           <  1:1032691/66‑1 (MQ=255)
                                                                                           atatTGAGTTTAAGGCTGTTCTCGCTTACGGACCGAAGAAGTAACCGTTGTGCg                                                                                                                    >  1:1153563/1‑54 (MQ=255)
                                                                                           atatTGAGTTTAAGGCTGTTCTCGCTTACGGACCGAAGAAGTAACCGTTGTGCg                                                                                                                    <  2:1153563/54‑1 (MQ=255)
                                                                                           atatTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                                    <  2:1544926/54‑1 (MQ=255)
                                                                                           atatTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                                    >  1:1544926/1‑54 (MQ=255)
                                                                                            tatTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATaa                              <  1:1285270/139‑1 (MQ=255)
                                                                                                gAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTgg                          >  1:558950/1‑139 (MQ=255)
                                                                                                   tttAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa                                          >  1:1251619/1‑120 (MQ=255)
                                                                                                   tttAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa                                          <  2:1251619/120‑1 (MQ=255)
                                                                                                   tttAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCa                       >  1:1001156/1‑139 (MQ=255)
                                                                                                    ttAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                                    >  2:1824988/1‑45 (MQ=255)
                                                                                                    ttAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                                    <  1:1824988/45‑1 (MQ=255)
                                                                                                      aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                               <  2:25663/64‑1 (MQ=255)
                                                                                                      aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtg                                                                                  <  1:836226/77‑1 (MQ=255)
                                                                                                      aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtg                                                                                  >  2:836226/1‑77 (MQ=255)
                                                                                                      aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGAACGGGGAGTAATCCCATAAGCGCt                                                                                               >  1:25663/1‑64 (MQ=255)
                                                                                                         gCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                     <  2:994738/71‑1 (MQ=255)
                                                                                                         gCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                     >  1:994738/1‑71 (MQ=255)
                                                                                                         gCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAAtt                 <  2:1001156/139‑1 (MQ=255)
                                                                                                             ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                                    <  2:1486164/36‑1 (MQ=255)
                                                                                                             ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                                    >  1:1486164/1‑36 (MQ=255)
                                                                                                             ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAAt                                                                                                           >  2:1525682/1‑45 (MQ=255)
                                                                                                             ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAAt                                                                                                           <  1:1525682/45‑1 (MQ=255)
                                                                                                             ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcgc                                                                                                <  2:417754/56‑1 (MQ=255)
                                                                                                             ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcgc                                                                                                >  1:417754/1‑56 (MQ=255)
                                                                                                             ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                               >  2:1516578/1‑57 (MQ=255)
                                                                                                             ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                               <  1:1516578/57‑1 (MQ=255)
                                                                                                             ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                                                            >  2:1616374/1‑76 (MQ=255)
                                                                                                             ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                                                            <  1:1616374/76‑1 (MQ=255)
                                                                                                             ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                                                       <  1:841119/97‑1 (MQ=255)
                                                                                                             ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                                                       >  2:841119/1‑97 (MQ=255)
                                                                                                             ttCTCGCTTACGAACCAAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgg             <  2:143686/139‑1 (MQ=255)
                                                                                                                       cGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                                                                   <  2:428198/75‑1 (MQ=255)
                                                                                                                       cGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                                                                   >  1:428198/1‑75 (MQ=255)
                                                                                                                       cGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg                      <  1:225713/120‑1 (MQ=255)
                                                                                                                       cGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg                      >  2:225713/1‑120 (MQ=255)
                                                                                                                       cGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTgg   >  1:36126/1‑139 (MQ=255)
                                                                                                                       cGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTgg   >  2:1023716/1‑139 (MQ=255)
                                                                                                                        gAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGt  <  1:1113715/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                            
CTTGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT  >  NC_000913/607854‑608112

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