Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,931,354 G→A G119G (GGC→GGT)  fucA ← L‑fuculose‑1‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,931,3540GA99.5% 711.0 / NA 212G119G (GGC→GGT) fucAL‑fuculose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (100/110);  total (101/111)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 231 total aligned reads displayed.

TGCCTTACGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAAC  >  NC_000913/2931264‑2931443
                                                                                          |                                                                                         
tGCCTTACGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTaccacc                                                                                         <  1:862520/93‑1 (MQ=255)
tGCCTTACGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTaccacc                                                                                         >  1:1695796/1‑93 (MQ=255)
tGCCTTACGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTaccacc                                                                                         >  1:3841081/1‑93 (MQ=255)
 gCCTTACGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTAccacca                                                                                        <  2:806822/93‑1 (MQ=255)
   cTTACGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAgc                                                                                      <  1:3408952/93‑1 (MQ=255)
    ttACGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAgcc                                                                                     >  1:1571862/1‑93 (MQ=255)
    ttACGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAgcc                                                                                     >  2:2755418/1‑93 (MQ=255)
     tACGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAgccg                                                                                    >  1:265872/1‑93 (MQ=255)
     tACGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAgccg                                                                                    >  1:700366/1‑93 (MQ=255)
        gATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGccgccg                                                                                 >  1:2538409/1‑93 (MQ=255)
        gATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGccgccg                                                                                 >  2:4062181/1‑93 (MQ=255)
              tGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCa                                                                           <  2:133527/93‑1 (MQ=255)
                 gagCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGt                                                                        <  1:3675063/93‑1 (MQ=255)
                 gagCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGt                                                                        <  2:3612633/93‑1 (MQ=255)
                 gagCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGt                                                                        <  1:4093840/93‑1 (MQ=255)
                  agCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTa                                                                       <  2:2795948/93‑1 (MQ=255)
                  agCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTa                                                                       <  2:3218948/93‑1 (MQ=255)
                    ccAGCGCAACATGTTCAGAAAGTTCGCGTGTTCAAACGGTCGCATAGGGAGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGt                                                                     <  2:4681282/93‑1 (MQ=255)
                     cAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTg                                                                    >  1:2116327/1‑93 (MQ=255)
                      aGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGa                                                                   <  2:1665811/93‑1 (MQ=255)
                       gcgcAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGaa                                                                  <  2:1223012/93‑1 (MQ=255)
                       gcgcAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGaa                                                                  >  1:3005668/1‑93 (MQ=255)
                        cgcAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAAt                                                                 <  2:3405201/93‑1 (MQ=255)
                        cgcAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAAt                                                                 <  2:3644754/93‑1 (MQ=255)
                          cAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAg                                                               <  2:1871472/93‑1 (MQ=255)
                           aaCATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGc                                                              >  1:277899/1‑93 (MQ=255)
                           aaCATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGc                                                              >  2:2910303/1‑93 (MQ=255)
                           aaCATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGc                                                              <  2:3324700/93‑1 (MQ=255)
                           aaCATGTTCAGAAAGTTCGCGTGTTCCAAAGGCCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGc                                                              <  1:2250796/93‑1 (MQ=255)
                            aCATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCg                                                             >  1:4183678/1‑93 (MQ=255)
                            aCATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCg                                                             >  1:4618470/1‑93 (MQ=255)
                            aCATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCg                                                             <  2:4448690/93‑1 (MQ=255)
                             cATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCgg                                                            >  1:1093644/1‑93 (MQ=255)
                              aTGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCggg                                                           <  1:1290534/93‑1 (MQ=255)
                              aTGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCggg                                                           >  2:4409831/1‑93 (MQ=255)
                                gTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGga                                                         <  1:2400724/93‑1 (MQ=255)
                                 ttCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGgat                                                        <  1:3755072/93‑1 (MQ=255)
                                 ttCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGgat                                                        <  2:2113328/93‑1 (MQ=255)
                                 ttCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGgat                                                        >  1:1649322/1‑93 (MQ=255)
                                 ttCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGgat                                                        <  1:4543800/93‑1 (MQ=255)
                                 tgcAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGgat                                                        <  2:1097915/91‑1 (MQ=255)
                                     gAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTCCCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGgatcgat                                                    <  2:1868158/93‑1 (MQ=255)
                                       aaGTTTGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGatcgatcg                                                  >  2:2600275/1‑93 (MQ=255)
                                       aaGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGatcgatcg                                                  >  1:3397449/1‑93 (MQ=255)
                                       aaGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGatcgatcg                                                  >  2:3962962/1‑93 (MQ=255)
                                         gTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGt                                                >  2:1920033/1‑93 (MQ=255)
                                         gTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGt                                                <  2:4066488/93‑1 (MQ=255)
                                           tCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTa                                              <  1:4421773/93‑1 (MQ=255)
                                           tCGCGTGTGCAAAAGGCCGCATAAGGCACGCAAGGAAGAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTa                                              <  2:4093002/93‑1 (MQ=255)
                                            cgcgTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTaa                                             <  2:3173739/93‑1 (MQ=255)
                                             gcgTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAg                                            >  1:1881475/1‑93 (MQ=255)
                                             gcgTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAg                                            >  1:3818607/1‑93 (MQ=255)
                                             gcgTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAg                                            >  2:2480215/1‑93 (MQ=255)
                                                tgtTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAAt                                         >  2:1276353/1‑93 (MQ=255)
                                                  ttCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATgg                                       <  1:4022452/93‑1 (MQ=255)
                                                    ccAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGaa                                     >  2:4218454/1‑93 (MQ=255)
                                                     cAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGaaa                                    <  2:1692094/93‑1 (MQ=255)
                                                      aaaGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAAc                                   <  2:4479424/93‑1 (MQ=255)
                                                      aaaGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAAc                                   >  1:3607519/1‑93 (MQ=255)
                                                      aaaGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGCGAATAGCGGGGATCGATCGGTTAAGAATGGAAAc                                   >  2:2254604/1‑93 (MQ=255)
                                                       aaGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACt                                  >  2:2967742/1‑93 (MQ=255)
                                                       aaGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACt                                  <  2:3025875/93‑1 (MQ=255)
                                                       aaGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACt                                  <  1:81533/93‑1 (MQ=255)
                                                       aaGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACt                                  >  1:3938212/1‑93 (MQ=255)
                                                       aaGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACt                                  >  1:475237/1‑93 (MQ=255)
                                                        aGGTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTg                                 <  2:3060182/93‑1 (MQ=255)
                                                         ggTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGc                                >  2:4241788/1‑93 (MQ=255)
                                                          gTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGcc                               >  2:2211552/1‑93 (MQ=255)
                                                          gTCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGcc                               >  2:4161783/1‑93 (MQ=255)
                                                           tCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCg                              >  1:3996488/1‑93 (MQ=255)
                                                           tCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCg                              <  1:1228715/93‑1 (MQ=255)
                                                           tCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCg                              <  1:112405/93‑1 (MQ=255)
                                                           tCGCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCg                              >  1:519636/1‑93 (MQ=255)
                                                            cGCATAAGGCGCGCAAGGAATAGAATTACCACCGGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGt                             <  2:3122055/93‑1 (MQ=255)
                                                             gCATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTg                            >  1:1564612/1‑93 (MQ=255)
                                                              cATAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAGTAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTga                           >  1:1600154/1‑92 (MQ=255)
                                                               aTAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCa                          >  1:4232821/1‑93 (MQ=255)
                                                               aTAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCa                          >  1:1020498/1‑93 (MQ=255)
                                                               aTAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCa                          <  2:1004892/93‑1 (MQ=255)
                                                                tAAGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCaa                         <  1:2743968/93‑1 (MQ=255)
                                                                 aaGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAAt                        <  1:1352368/93‑1 (MQ=255)
                                                                  aGGCGCGCAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATg                       >  1:1400861/1‑93 (MQ=255)
                                                                    gcgcgcAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGaa                     <  2:637447/93‑1 (MQ=255)
                                                                     cgcgcAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAAc                    >  1:3230318/1‑93 (MQ=255)
                                                                     cgcgcAAGAAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAAc                    <  1:3430491/93‑1 (MQ=255)
                                                                       cgcAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACgg                  <  1:4394580/93‑1 (MQ=255)
                                                                         cAAGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCa                <  2:3897596/93‑1 (MQ=255)
                                                                          aaGGAATAGAATTACCACCAGCCGCCGCAATCTTGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCAt               >  2:1575990/1‑93 (MQ=255)
                                                                          aaGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAAGAGCGGGGATCGATCGGCTAAGAATGGAAACTGCCGGGCAATTAATTGCAt               >  1:2515472/1‑93 (MQ=255)
                                                                           aGGAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATg              >  1:3103738/1‑93 (MQ=255)
                                                                            ggAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATgt             >  1:640831/1‑92 (MQ=255)
                                                                             gAATAGAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGAAATGAACGGCATGAt            >  2:4575092/1‑93 (MQ=255)
                                                                                  gAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGa       >  1:1982055/1‑93 (MQ=255)
                                                                                  gAATTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGa       >  2:3482887/1‑93 (MQ=255)
                                                                                   aaTTCCCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaa      <  1:219652/93‑1 (MQ=255)
                                                                                   aaTTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaa      <  2:870319/93‑1 (MQ=255)
                                                                                   aaTTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaa      <  2:4329663/93‑1 (MQ=255)
                                                                                    aTTACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaac     >  1:3193936/1‑93 (MQ=255)
                                                                                      tACCACCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaacaa   <  1:3644427/93‑1 (MQ=255)
                                                                                       accaccAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaacaac  <  2:2724861/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGCCTTACGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAAC  >  NC_000913/2931264‑2931443

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑