Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,694,050 A→G Y31C (TAC→TGC)  yfhD → putative periplasmic binding transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,694,0500AG98.8% 510.6 / ‑2.7 161Y31C (TAC→TGC) yfhDputative periplasmic binding transport protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (88/71);  total (88/73)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

Only 100 of 167 total aligned reads displayed.

GAGCACCATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTACCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACCTGCTGGCGG  >  NC_000913/2693961‑2694139
                                                                                         |                                                                                         
gAGCACCATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCt                                                                                        >  2:744810/1‑93 (MQ=255)
gAGCACCATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCt                                                                                        <  1:1513140/93‑1 (MQ=255)
 aGCACCATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTc                                                                                       >  2:3047226/1‑93 (MQ=255)
    aCCATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGc                                                                                    >  1:653331/1‑93 (MQ=255)
    aCCATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGc                                                                                    <  1:3641660/93‑1 (MQ=255)
     ccATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGTCGATTGCCTCGGCg                                                                                   >  2:3686833/1‑93 (MQ=255)
      cATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGCAACAGTTTGCCGATTGCCTCGGCGt                                                                                  >  2:2243524/1‑93 (MQ=255)
       aTTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTa                                                                                 >  2:2111397/1‑93 (MQ=255)
              cTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGa                                                                          >  1:4526207/1‑93 (MQ=255)
              cTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTACCTCGGCGTAAAACTGa                                                                          <  1:2326473/93‑1 (MQ=255)
              cTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGAGTGCCTCGGCGTAAAACTGa                                                                          >  2:4350270/1‑93 (MQ=255)
               tCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGaa                                                                         >  2:161447/1‑93 (MQ=255)
                ccccTGACTTATAACGCAATCAACGGGAAACCTTTTGGCCTGGATACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAg                                                                       >  2:2254750/1‑93 (MQ=255)
                ccccTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGaaa                                                                        >  2:2254751/1‑93 (MQ=255)
                 cccTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAg                                                                       >  1:946060/1‑93 (MQ=255)
                 cccTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAg                                                                       <  1:1859679/93‑1 (MQ=255)
                  ccTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAgg                                                                      >  2:3470369/1‑92 (MQ=255)
                  ccTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGt                                                                      >  1:3856406/1‑93 (MQ=255)
                  ccTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGt                                                                      <  1:353013/93‑1 (MQ=255)
                  ccTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGt                                                                      <  2:4395199/93‑1 (MQ=255)
                  catGACTTATAACGAAATCAACGGGAATCCTTTTGGCCTGGGTTACGAACTGGAGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGt                                                                      <  1:353012/91‑1 (MQ=255)
                   cTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTg                                                                     >  1:2238645/1‑93 (MQ=255)
                    tGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGa                                                                    <  2:3410465/93‑1 (MQ=255)
                    tGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGa                                                                    >  1:362236/1‑93 (MQ=255)
                    tGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGa                                                                    <  1:1051836/93‑1 (MQ=255)
                     gACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGAc                                                                   >  1:1738113/1‑93 (MQ=255)
                     gACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGAc                                                                   >  1:430376/1‑93 (MQ=255)
                     gACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGAc                                                                   >  1:546082/1‑93 (MQ=255)
                     gACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGAc                                                                   <  1:77269/93‑1 (MQ=255)
                      aCTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGAcc                                                                  <  2:3434225/93‑1 (MQ=255)
                       cTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCg                                                                 <  2:2336967/93‑1 (MQ=255)
                        ttATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGt                                                                <  1:2653910/93‑1 (MQ=255)
                        ttATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGt                                                                >  2:3248027/1‑93 (MQ=255)
                         tataACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTg                                                               >  2:638854/1‑93 (MQ=255)
                          ataACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAAGAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGc                                                              >  2:2169846/1‑93 (MQ=255)
                              cGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCa                                                          >  1:1018900/1‑93 (MQ=255)
                                aaaTCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGa                                                        >  2:3555023/1‑93 (MQ=255)
                                  aTCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAat                                                      >  2:274651/1‑93 (MQ=255)
                                  aTCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAat                                                      >  1:2045755/1‑93 (MQ=255)
                                  aTCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAat                                                      <  2:3993090/93‑1 (MQ=255)
                                    cAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAatat                                                    >  2:2489118/1‑93 (MQ=255)
                                    cAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAatat                                                    >  2:1051006/1‑93 (MQ=255)
                                    cAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAatat                                                    >  1:3282775/1‑93 (MQ=255)
                                     aaCGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATc                                                   <  2:2170500/93‑1 (MQ=255)
                                       cGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcag                                                 <  2:399771/93‑1 (MQ=255)
                                        gggAAACCTTTTTGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcagc                                                >  1:1189178/1‑93 (MQ=255)
                                         ggAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAACACTGAAAGTGACCGTGCGGCAGAATATcagcc                                               >  1:3519023/1‑93 (MQ=255)
                                         ggAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcagcc                                               >  1:2587223/1‑93 (MQ=255)
                                         ggAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcagcc                                               <  1:4374874/93‑1 (MQ=255)
                                         ggAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcagcc                                               <  2:4291011/93‑1 (MQ=255)
                                          gAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcagcca                                              <  2:1117200/93‑1 (MQ=255)
                                            aaCCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcagccagc                                            <  1:3035152/93‑1 (MQ=255)
                                            aaCCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcagccagc                                            >  2:2926920/1‑93 (MQ=255)
                                                ttttGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGtt                                        >  2:3844770/1‑93 (MQ=255)
                                                ttttGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGtt                                        >  1:3776998/1‑93 (MQ=255)
                                                 tttGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCGGCCAGCTGttt                                       <  1:3679454/93‑1 (MQ=255)
                                                    ggCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTgac                                    >  1:1907273/1‑93 (MQ=255)
                                                     gCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTgacg                                   >  2:418395/1‑93 (MQ=255)
                                                      ccTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTgacgc                                  >  2:1637076/1‑92 (MQ=255)
                                                       cTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTgacgac                                 >  2:4613801/1‑93 (MQ=255)
                                                       cTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTgacgac                                 >  2:500961/1‑93 (MQ=255)
                                                         ggATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCt                               <  1:1700999/93‑1 (MQ=255)
                                                          gATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCtt                              >  1:3437415/1‑93 (MQ=255)
                                                           aTTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTg                             >  1:4114331/1‑93 (MQ=255)
                                                            ttACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTTACGACCTTGa                            <  1:3871490/93‑1 (MQ=255)
                                                            ttACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGa                            <  2:2450392/93‑1 (MQ=255)
                                                             tACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGAt                           <  2:623841/93‑1 (MQ=255)
                                                             tACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGAt                           <  2:2582933/93‑1 (MQ=255)
                                                              aCGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATa                          >  2:1597018/1‑93 (MQ=255)
                                                               cGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATaa                         <  1:4326558/93‑1 (MQ=255)
                                                               cGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATaa                         >  1:76227/1‑93 (MQ=255)
                                                                gAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAAt                        <  2:2392374/93‑1 (MQ=255)
                                                                gAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAAt                        <  1:2659293/93‑1 (MQ=255)
                                                                gAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAAt                        <  1:1672478/93‑1 (MQ=255)
                                                                 aaCTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATg                       >  2:1775769/1‑93 (MQ=255)
                                                                  aCTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATgg                      <  1:602839/93‑1 (MQ=255)
                                                                   cTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGt                     >  1:321131/1‑93 (MQ=255)
                                                                    tGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTa                    >  1:4570049/1‑93 (MQ=255)
                                                                     ggCGAAATAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTaa                   >  2:3744734/1‑93 (MQ=255)
                                                                       cGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACg                 >  1:2532279/1‑93 (MQ=255)
                                                                        gAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCGGCCAGCTGTTTGACGACCTTGATAATGGTAACGc                >  2:2353886/1‑93 (MQ=255)
                                                                        gAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGc                >  2:4552674/1‑93 (MQ=255)
                                                                        gAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGc                >  2:1841955/1‑93 (MQ=255)
                                                                         aaaCAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGcc               <  1:2349864/93‑1 (MQ=255)
                                                                          aaCAGTTTGCTGATTGCCTCGGCGTAAAACTGAAAGTGAACGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCg              <  1:2680986/93‑1 (MQ=255)
                                                                            cAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGAc            <  1:10985/93‑1 (MQ=255)
                                                                            cAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGAc            <  1:2598466/93‑1 (MQ=255)
                                                                             aGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGAcc           <  2:1399189/93‑1 (MQ=255)
                                                                             aGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGAcc           <  2:3942617/93‑1 (MQ=255)
                                                                             aGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGAcc           <  1:335123/93‑1 (MQ=255)
                                                                              gTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACct          >  2:2613579/1‑93 (MQ=255)
                                                                                ttGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACctgc        <  2:3113417/93‑1 (MQ=255)
                                                                                ttGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACctgc        <  1:3490267/93‑1 (MQ=255)
                                                                                ttGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACctgc        <  1:4438928/93‑1 (MQ=255)
                                                                                 tGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACctgct       <  1:1018191/93‑1 (MQ=255)
                                                                                 tGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACctgct       >  2:2866850/1‑93 (MQ=255)
                                                                                   ccGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACCTGCTgg     <  1:1466630/93‑1 (MQ=255)
                                                                                    cGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACCTGCTggc    >  2:3851199/1‑93 (MQ=255)
                                                                                      aTTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACCTGCTggcgg  >  2:4136399/1‑93 (MQ=255)
                                                                                      aTTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACCTGCTggcgg  >  1:1630633/1‑93 (MQ=255)
                                                                                         |                                                                                         
GAGCACCATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTACCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACCTGCTGGCGG  >  NC_000913/2693961‑2694139

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