Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,967,307 A→G D222G (GAC→GGC)  wzzE → putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,967,3070AG98.6% 438.6 / ‑3.4 140D222G (GAC→GGC) wzzEputative transport protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (63/74);  total (65/75)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 146 total aligned reads displayed.

ATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTA  >  NC_000913/3967217‑3967398
                                                                                          |                                                                                           
atgatgAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGcc                                                                                           >  2:2331252/1‑93 (MQ=255)
atgatgAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGcc                                                                                           <  1:3198630/93‑1 (MQ=255)
 tgatgaGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGccg                                                                                          >  2:3078451/1‑93 (MQ=255)
  gatgaGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGccgc                                                                                         >  2:2003768/1‑93 (MQ=255)
  gatgaGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGccgc                                                                                         <  1:4039695/93‑1 (MQ=255)
  gatgaGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGccgc                                                                                         >  1:2693439/1‑93 (MQ=255)
  gatgaGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGccgc                                                                                         <  1:2014855/93‑1 (MQ=255)
      aGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGa                                                                                     <  2:1061701/93‑1 (MQ=255)
       gCTGAAAGGCGCATGGGCGGCTCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGAt                                                                                    <  2:3802881/93‑1 (MQ=255)
       gCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGAt                                                                                    <  1:1087281/93‑1 (MQ=255)
       gCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGAt                                                                                    <  1:1087280/93‑1 (MQ=255)
       gCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGAt                                                                                    <  1:3050611/93‑1 (MQ=255)
        cTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATg                                                                                   <  2:955838/93‑1 (MQ=255)
        cTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATg                                                                                   <  2:2530867/93‑1 (MQ=255)
         tGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGa                                                                                  >  1:4630281/1‑93 (MQ=255)
          gAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGaa                                                                                 <  1:3716141/93‑1 (MQ=255)
          gAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGaa                                                                                 <  1:2080817/93‑1 (MQ=255)
          gAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGaa                                                                                 >  2:2902519/1‑93 (MQ=255)
          gAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGaa                                                                                 >  1:2770868/1‑93 (MQ=255)
           aaaGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAAc                                                                                >  2:1984632/1‑93 (MQ=255)
           aaaGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAAc                                                                                <  1:3439270/93‑1 (MQ=255)
            aaGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACa                                                                               <  2:4289130/93‑1 (MQ=255)
            aaGGCGCATGGGCGGAGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACa                                                                               <  1:4279969/93‑1 (MQ=255)
             aGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAg                                                                              >  1:3336135/1‑93 (MQ=255)
             aGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAg                                                                              <  1:1576566/93‑1 (MQ=255)
                cgcATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCAt                                                                           <  2:3324174/93‑1 (MQ=255)
                 gcATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCAtt                                                                          <  1:414607/93‑1 (MQ=255)
                   aTGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGa                                                                        >  2:3325671/1‑93 (MQ=255)
                   aTGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGa                                                                        <  1:1193897/93‑1 (MQ=255)
                    tGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAg                                                                       >  1:4275161/1‑93 (MQ=255)
                    tGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAg                                                                       >  1:4540451/1‑93 (MQ=255)
                      ggcggcGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCa                                                                     <  1:151304/93‑1 (MQ=255)
                      ggcggcGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCa                                                                     <  1:3393881/93‑1 (MQ=255)
                      ggcggcGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCa                                                                     <  2:720566/93‑1 (MQ=255)
                         ggcgCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGGCGGGTGAACAGCATTGAGCAggg                                                                  >  1:1730111/1‑92 (MQ=255)
                         ggcgCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGgc                                                                  <  2:3268529/93‑1 (MQ=255)
                          gcgcGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGgcg                                                                 >  2:2901364/1‑93 (MQ=255)
                            gcgTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCt                                                               >  2:1316323/1‑93 (MQ=255)
                               tACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGaa                                                            >  1:1946039/1‑93 (MQ=255)
                               tACCATCAGGATGAAAGCTCAGGTGAGGCGTCAGAAGGAGGTGGCGAAAGCCATCTACGGCCGCCGGTTGACCAGCATTAGGAAGGCGCTGaa                                                            <  2:4211589/93‑1 (MQ=255)
                                aCCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGaaa                                                           >  1:762636/1‑93 (MQ=255)
                                aCCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGaaa                                                           >  2:2355733/1‑93 (MQ=255)
                                 ccATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGaaaa                                                          <  1:97573/93‑1 (MQ=255)
                                 ccATCCAGATGAAAGCTCAGGTGAAGCATCAGGAAGAGGTGGCGAAAGCCATCTACGCCCGCCGGATGAACAGCATTGAGCAGGCGCTGaaaa                                                          <  2:1357904/93‑1 (MQ=255)
                                  cATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAACAt                                                         >  1:59049/1‑93 (MQ=255)
                                   aTCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAAtt                                                        >  2:2821328/1‑93 (MQ=255)
                                   aTCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAAtt                                                        >  2:3175025/1‑93 (MQ=255)
                                     ccAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGc                                                      <  2:3275916/93‑1 (MQ=255)
                                        gATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGa                                                   >  2:2163755/1‑93 (MQ=255)
                                        gATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGa                                                   <  1:4718954/93‑1 (MQ=255)
                                        gATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGa                                                   <  2:3882144/93‑1 (MQ=255)
                                        gATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGa                                                   <  2:347040/93‑1 (MQ=255)
                                          tGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGagc                                                 >  2:4617613/1‑93 (MQ=255)
                                          tGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGagc                                                 <  2:1563225/93‑1 (MQ=255)
                                           gAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGagca                                                >  2:892056/1‑93 (MQ=255)
                                           gAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGagca                                                <  2:1610612/93‑1 (MQ=255)
                                             aaGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGagcagc                                              >  2:350387/1‑93 (MQ=255)
                                              aGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGagcagaa                                             >  1:1588630/1‑91 (MQ=255)
                                              aGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAgcagca                                             >  1:2302235/1‑93 (MQ=255)
                                               gCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAgcagcaa                                            >  2:1542548/1‑92 (MQ=255)
                                                 tCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCataa                                          <  1:2575682/93‑1 (MQ=255)
                                                  cAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCataat                                         <  2:4029677/93‑1 (MQ=255)
                                                    ggTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAatat                                       >  1:3526373/1‑93 (MQ=255)
                                                      tGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATAttt                                     <  2:3116732/93‑1 (MQ=255)
                                                       gAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTc                                    >  1:1636367/1‑93 (MQ=255)
                                                       gAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTc                                    <  1:3715633/93‑1 (MQ=255)
                                                       gAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTc                                    <  1:2835644/93‑1 (MQ=255)
                                                       gAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTc                                    >  1:2349159/1‑93 (MQ=255)
                                                         aGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTcgc                                  <  2:565460/93‑1 (MQ=255)
                                                          gCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTcgcg                                 >  1:301138/1‑93 (MQ=255)
                                                          gCGTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAAAATTTcgcg                                 >  1:262939/1‑93 (MQ=255)
                                                            gTCAGGAAGCGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCa                               >  2:2556485/1‑93 (MQ=255)
                                                            gTCAGGAAGAGGTTGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCa                               <  2:2743230/93‑1 (MQ=255)
                                                            gTCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCa                               <  1:759127/93‑1 (MQ=255)
                                                             tCAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAg                              >  1:4541548/1‑93 (MQ=255)
                                                              cAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGt                             <  1:3858999/93‑1 (MQ=255)
                                                              cAGGAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGt                             >  1:82342/1‑93 (MQ=255)
                                                                ggAAGAGGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGc                           >  1:1859096/1‑93 (MQ=255)
                                                                   agagGTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGAc                        >  2:4445530/1‑93 (MQ=255)
                                                                      ggTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGa                     <  2:671159/93‑1 (MQ=255)
                                                                       gTGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGAt                    <  1:1580623/93‑1 (MQ=255)
                                                                        tGGCGAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATg                   <  2:1918102/93‑1 (MQ=255)
                                                                            gAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTGcc               <  1:3824645/93‑1 (MQ=255)
                                                                            gAAAGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTAcc               >  1:76065/1‑93 (MQ=255)
                                                                               aGCCATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGc            <  2:4451837/93‑1 (MQ=255)
                                                                                  cATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGa         <  1:1696066/93‑1 (MQ=255)
                                                                                  cATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGa         <  1:4113343/93‑1 (MQ=255)
                                                                                  cATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGa         <  2:1684511/93‑1 (MQ=255)
                                                                                  cATCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGa         >  2:3910189/1‑93 (MQ=255)
                                                                                   aTCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAg        >  1:1495772/1‑93 (MQ=255)
                                                                                    tCTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAgg       >  2:4373242/1‑93 (MQ=255)
                                                                                     cTACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGa      <  2:766890/93‑1 (MQ=255)
                                                                                      tACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGaa     <  1:2561902/93‑1 (MQ=255)
                                                                                      tACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGaa     >  1:753247/1‑93 (MQ=255)
                                                                                      tACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGaa     <  1:677470/93‑1 (MQ=255)
                                                                                      tACGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGaa     <  2:3827184/93‑1 (MQ=255)
                                                                                       aCGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAAt    <  2:4323836/93‑1 (MQ=255)
                                                                                       aCGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAAt    >  1:2175059/1‑93 (MQ=255)
                                                                                       aCGGCCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCCAGGAAt    <  2:4607200/93‑1 (MQ=255)
                                                                                         gACCGCCAGATGAAGAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACGGATGTACCTGCCGCGGAATTa  >  2:2092732/1‑93 (MQ=255)
                                                                                          |                                                                                           
ATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTA  >  NC_000913/3967217‑3967398

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