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,0500AG99.4% 543.1 / ‑5.1 175Y31C (TAC→TGC) yfhDputative periplasmic binding transport protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (107/67);  total (107/68)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 179 total aligned reads displayed.

TGAGCACCATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTACCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACCTGCTGGCGG  >  NC_000913/2693960‑2694139
                                                                                          |                                                                                         
tgAGCACCATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGcc                                                                                         <  2:1091446/93‑1 (MQ=255)
     aCCATTCATTCTCCCCTGACTTATAACGAAATCAGCGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGc                                                                                    <  1:76918/93‑1 (MQ=255)
     aCCATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGc                                                                                    <  2:2399264/93‑1 (MQ=255)
     aCCATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGc                                                                                    <  1:4552266/93‑1 (MQ=255)
     aCCATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGc                                                                                    >  1:1913081/1‑93 (MQ=255)
     aCCATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGc                                                                                    >  1:41322/1‑93 (MQ=255)
       cATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGt                                                                                  >  2:2879805/1‑93 (MQ=255)
       cATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGt                                                                                  >  1:2086045/1‑93 (MQ=255)
         ttCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTaa                                                                                >  2:2564419/1‑93 (MQ=255)
          tCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTaaa                                                                               >  2:773768/1‑93 (MQ=255)
          tCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTaaa                                                                               >  1:2436066/1‑93 (MQ=255)
           cATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTaaaa                                                                              <  1:4347976/93‑1 (MQ=255)
           cATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTaaaa                                                                              >  2:3994819/1‑93 (MQ=255)
            aTACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAAc                                                                             <  2:1825434/93‑1 (MQ=255)
             tACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACt                                                                            <  2:3589990/93‑1 (MQ=255)
             tACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACt                                                                            >  2:4528184/1‑93 (MQ=255)
               cTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGa                                                                          >  2:3040294/1‑93 (MQ=255)
                 ccccTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGaaa                                                                        >  2:1295500/1‑93 (MQ=255)
                 ccccTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGaaa                                                                        <  1:61241/93‑1 (MQ=255)
                    cTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTg                                                                     >  1:498308/1‑93 (MQ=255)
                    cTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTg                                                                     <  2:1147311/93‑1 (MQ=255)
                     tGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGa                                                                    <  1:3810119/93‑1 (MQ=255)
                      gACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGAc                                                                   <  2:1769053/93‑1 (MQ=255)
                      gACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGAc                                                                   <  1:1533145/93‑1 (MQ=255)
                       aCTTATAACGAAATCAACGGGAAGCCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGAcc                                                                  <  1:1177702/93‑1 (MQ=255)
                        cTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGATCg                                                                 >  2:3290705/1‑93 (MQ=255)
                        cTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCg                                                                 >  1:353164/1‑93 (MQ=255)
                           ataACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGc                                                              >  1:3251360/1‑93 (MQ=255)
                           ataACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGc                                                              <  1:4078519/93‑1 (MQ=255)
                            taACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCg                                                             >  1:395465/1‑93 (MQ=255)
                            taACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCg                                                             >  1:4269305/1‑93 (MQ=255)
                             aaCGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCgg                                                            >  2:4120529/1‑93 (MQ=255)
                                gAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAg                                                         >  2:2745549/1‑93 (MQ=255)
                                gAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAg                                                         <  1:3000206/93‑1 (MQ=255)
                                  aaTCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTACAACTGAAAGTGACCGTGCGGCAGaa                                                       >  1:2156479/1‑93 (MQ=255)
                                  aaTCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGaa                                                       >  1:3675751/1‑93 (MQ=255)
                                  aaTCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGaa                                                       <  1:848600/93‑1 (MQ=255)
                                   aTCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAat                                                      >  1:4482185/1‑93 (MQ=255)
                                    tCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAata                                                     <  2:2959566/93‑1 (MQ=255)
                                       aCGGGAACCCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGCGACCGTGCGGCAGAgtataa                                                  >  1:3016182/1‑91 (MQ=255)
                                       aCGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATca                                                  >  2:2423610/1‑93 (MQ=255)
                                       aCGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATca                                                  <  2:4329269/93‑1 (MQ=255)
                                        cGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcag                                                 >  2:4056776/1‑93 (MQ=255)
                                         gggAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcagc                                                >  1:1112304/1‑93 (MQ=255)
                                          ggAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcagcc                                               >  2:2116700/1‑93 (MQ=255)
                                            aaaCCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcagccag                                             >  1:107558/1‑93 (MQ=255)
                                             aaCCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcagccagc                                            <  2:3785159/93‑1 (MQ=255)
                                             aaCCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcagccagc                                            >  1:1861113/1‑93 (MQ=255)
                                             aaCCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcagccagc                                            >  2:4476581/1‑93 (MQ=255)
                                              aCCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCt                                           >  1:949011/1‑93 (MQ=255)
                                              aCCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCt                                           >  1:700354/1‑93 (MQ=255)
                                              aCCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCt                                           >  1:139501/1‑93 (MQ=255)
                                              aCCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCt                                           >  1:3989826/1‑93 (MQ=255)
                                              aCCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGACAGCt                                           >  1:1634403/1‑93 (MQ=255)
                                               ccTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATcagccagacg                                          >  1:2972647/1‑90 (MQ=255)
                                                cTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGt                                         >  2:953640/1‑93 (MQ=255)
                                                cTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGt                                         >  1:1934240/1‑93 (MQ=255)
                                                cTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGt                                         >  1:3366581/1‑93 (MQ=255)
                                                 ttttGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGtt                                        <  1:4280665/93‑1 (MQ=255)
                                                    tGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTga                                     <  1:3293301/93‑1 (MQ=255)
                                                     ggCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGATGTTTgac                                    >  2:3890700/1‑93 (MQ=255)
                                                      gCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTgacg                                   >  2:4444361/1‑93 (MQ=255)
                                                      gCCTGGATTACGAACTGGCGAAACAGCTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTgacg                                   <  2:3796791/93‑1 (MQ=255)
                                                        cTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTgacgac                                 >  2:3703519/1‑93 (MQ=255)
                                                        cTGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTgacgac                                 >  1:371657/1‑93 (MQ=255)
                                                         tGGATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGAcc                                <  1:650192/93‑1 (MQ=255)
                                                           gATTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCtt                              >  2:1252972/1‑93 (MQ=255)
                                                            aTTACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTg                             >  1:3445135/1‑93 (MQ=255)
                                                             ttACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAGAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGa                            <  1:2017427/93‑1 (MQ=255)
                                                              tACGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGAt                           >  1:265162/1‑93 (MQ=255)
                                                               aCGAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATa                          <  2:1528431/93‑1 (MQ=255)
                                                                 gAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAAt                        <  1:4172757/93‑1 (MQ=255)
                                                                 gAACTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAAt                        >  1:3581225/1‑93 (MQ=255)
                                                                   aCTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATgg                      >  2:2672543/1‑93 (MQ=255)
                                                                    cTGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGt                     >  1:3254749/1‑93 (MQ=255)
                                                                     tGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTa                    <  1:2641102/93‑1 (MQ=255)
                                                                     tGGCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGGGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTa                    >  2:1406164/1‑93 (MQ=255)
                                                                      ggCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTaa                   <  2:4176124/93‑1 (MQ=255)
                                                                       gCGAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAAc                  <  2:519876/93‑1 (MQ=255)
                                                                         gAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGc                >  1:3311811/1‑93 (MQ=255)
                                                                         gAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGc                >  1:2399206/1‑93 (MQ=255)
                                                                         gAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGc                >  1:2526845/1‑93 (MQ=255)
                                                                         gAAACAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGc                <  1:1406273/93‑1 (MQ=255)
                                                                          aaaCAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGcc               >  2:2149803/1‑93 (MQ=255)
                                                                           aaCAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCg              >  2:4572462/1‑93 (MQ=255)
                                                                            aCAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGa             >  2:343167/1‑93 (MQ=255)
                                                                             cAGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGAc            <  2:2117048/93‑1 (MQ=255)
                                                                              aGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGAcc           <  2:1410875/93‑1 (MQ=255)
                                                                              aGTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGAcc           <  1:1224233/93‑1 (MQ=255)
                                                                               gTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACct          >  1:1454851/1‑93 (MQ=255)
                                                                               gTTTGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACct          >  2:1474915/1‑93 (MQ=255)
                                                                                tttGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACctg         >  1:2854961/1‑93 (MQ=255)
                                                                                 ttGCCGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACctgc        <  2:1443679/93‑1 (MQ=255)
                                                                                    ccGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACCTGCTgg     >  1:3761000/1‑93 (MQ=255)
                                                                                    ccGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACCTGCTgg     >  2:2340318/1‑93 (MQ=255)
                                                                                    ccGATCGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACCTGCTgg     <  2:4245891/93‑1 (MQ=255)
                                                                                     cGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACCTGCTggc    <  2:3983700/93‑1 (MQ=255)
                                                                                     cGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACCTGCTggc    <  1:2369545/93‑1 (MQ=255)
                                                                                     cGATTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACCTGCTggc    <  2:1922317/93‑1 (MQ=255)
                                                                                       aTTGCCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGATGTTTGACGACCTTGATAATGGTAACGCCGACCTGCTggcgg  >  2:415711/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGAGCACCATTCATACTCCCCTGACTTATAACGAAATCAACGGGAAACCTTTTGGCCTGGATTACGAACTGGCGAAACAGTTTGCCGATTACCTCGGCGTAAAACTGAAAGTGACCGTGCGGCAGAATATCAGCCAGCTGTTTGACGACCTTGATAATGGTAACGCCGACCTGCTGGCGG  >  NC_000913/2693960‑2694139

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