Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,459,293 (G)5→6 coding (955/1482 nt) gspE → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,459,2881.G96.5% 510.3 / 10.5 144coding (950/1482 nt)gspEgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base G (65/74);  total (67/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 155 total aligned reads displayed.

TGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCAT  >  NC_000913/3459189‑3459383
                                                                                                     |                                                                                                
tGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtggt                                                                                                   <  1:224902/100‑1 (MQ=255)
     ggCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTggggg                                                                                             >  2:4048256/1‑100 (MQ=255)
      gCATTGGGCATACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCTCCAGGACCCGGATGTCGTCATGGTgggggg                                                                                            <  2:2750354/100‑2 (MQ=255)
      gCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggggg                                                                                            <  2:2310099/100‑2 (MQ=255)
        aTTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTggggggaa                                                                                          >  1:1893719/1‑97 (MQ=255)
          tGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTggggggaaat                                                                                        >  2:698379/1‑95 (MQ=255)
          tGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTggggggaaat                                                                                        >  1:3038367/1‑95 (MQ=255)
          tGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGgaaat                                                                                       >  2:2278552/1‑98 (MQ=255)
           gggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAAtt                                                                                       >  1:415751/1‑100 (MQ=255)
           gggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTggggagaa                                                                                          >  1:1585423/1‑93 (MQ=255)
            ggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTc                                                                                      >  1:3978611/1‑100 (MQ=255)
              cAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGt                                                                                    <  2:1231755/100‑1 (MQ=255)
              cAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGt                                                                                    >  2:163199/1‑100 (MQ=255)
              cAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGt                                                                                    <  1:1252012/100‑1 (MQ=255)
               aaaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTg                                                                                   <  1:2132548/100‑1 (MQ=255)
               aaaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTg                                                                                   <  1:2581995/100‑1 (MQ=255)
                 aCGCAGGTGAATACCTGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGAt                                                                                 <  1:2429226/100‑1 (MQ=255)
                   gCAGGAGAATACCCGTGTGGATATGTCATTCGCTCGAGGCATGCGCGGCATACATCGCCAGGACCCGGATGTCGACATGGTGGGGGGCAATTCGTGATAc                                                                               <  2:616623/100‑1 (MQ=255)
                    cAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACa                                                                              <  2:2990364/100‑1 (MQ=255)
                      ggTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGa                                                                            <  1:2034633/100‑1 (MQ=255)
                      ggTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGa                                                                            >  2:1031031/1‑100 (MQ=255)
                      ggTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGa                                                                            <  1:3268291/100‑1 (MQ=255)
                      ggTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGa                                                                            >  1:1376699/1‑100 (MQ=255)
                      ggTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGa                                                                            >  2:3629317/1‑100 (MQ=255)
                      ggTGAATACCCGTGTGGATATGTCTTGCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGa                                                                            <  2:3890386/100‑1 (MQ=255)
                       gTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGaa                                                                           >  1:2537749/1‑100 (MQ=255)
                         gAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAAc                                                                         >  2:3408875/1‑100 (MQ=255)
                         gAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAAc                                                                         <  1:3718444/100‑1 (MQ=255)
                           aTACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACcg                                                                       <  2:2535629/100‑1 (MQ=255)
                            tACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACcgc                                                                      >  1:2128770/1‑100 (MQ=255)
                            tACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACcgc                                                                      >  2:3174418/1‑100 (MQ=255)
                             aCCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACcgcg                                                                     <  2:605111/100‑1 (MQ=255)
                              cccGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCgcgc                                                                    <  2:1292362/100‑1 (MQ=255)
                              cccGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCgcgc                                                                    >  1:1213291/1‑100 (MQ=255)
                              cccGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGCTGGGGGGAAATTCGTGATACAGAAACcgccc                                                                    <  1:2815393/100‑3 (MQ=255)
                               ccGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGGCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCa                                                                   <  2:1958640/100‑1 (MQ=255)
                                cGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAg                                                                  <  1:3082783/100‑1 (MQ=255)
                                    tgGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTg                                                              <  2:3040584/100‑1 (MQ=255)
                                      gATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCg                                                            >  2:807618/1‑100 (MQ=255)
                                       atatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCgg                                                           <  1:1440343/100‑1 (MQ=255)
                                        tatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGt                                                          >  2:1717523/1‑100 (MQ=255)
                                          tGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTc                                                        >  1:3667585/1‑100 (MQ=255)
                                           gTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCa                                                       >  1:3373305/1‑100 (MQ=255)
                                            tCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAg                                                      <  2:1893705/100‑1 (MQ=255)
                                            tCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAg                                                      <  1:3651263/100‑1 (MQ=255)
                                             cTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgg                                                     <  1:2957493/100‑1 (MQ=255)
                                               ttCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGcc                                                   >  2:1764604/1‑100 (MQ=255)
                                               ttCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGc                                                    >  2:4135351/1‑99 (MQ=255)
                                                tcgctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCt                                                  >  1:101213/1‑100 (MQ=255)
                                                tcgctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCt                                                  <  1:1415671/100‑1 (MQ=255)
                                                tcgctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCt                                                  >  2:2147695/1‑100 (MQ=255)
                                                  gctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCg                                                >  1:3480353/1‑100 (MQ=255)
                                                  gctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCg                                                >  1:1990337/1‑100 (MQ=255)
                                                   ctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGc                                               >  2:2588703/1‑100 (MQ=255)
                                                     cgcgGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCt                                              <  1:2848470/99‑1 (MQ=255)
                                                     cgcgGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTc                                             >  2:287276/1‑100 (MQ=255)
                                                     cgcgGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTc                                             >  2:3373518/1‑100 (MQ=255)
                                                      gcgGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCa                                            >  1:2319840/1‑100 (MQ=255)
                                                          ccTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCgg                                        >  1:320281/1‑100 (MQ=255)
                                                           cTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGc                                       <  1:2254387/100‑1 (MQ=255)
                                                            tGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGCTTGCGGTTCAGGCCTCGCTCACCg                                         >  2:1034626/1‑97 (MQ=255)
                                                            tGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGcc                                      <  1:3344464/100‑1 (MQ=255)
                                                             gcgcgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCa                                     >  1:967230/1‑100 (MQ=255)
                                                              cgcgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCAt                                    <  1:1031036/100‑1 (MQ=255)
                                                                cgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCt                                  >  1:1463874/1‑100 (MQ=255)
                                                                 gcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAGCCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTg                                 <  2:3581001/100‑1 (MQ=255)
                                                                  ccATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTgn                                <  2:1702266/100‑2 (MQ=255)
                                                                   cATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGt                               <  2:1826713/100‑1 (MQ=255)
                                                                   cATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGt                               <  1:3766092/100‑1 (MQ=255)
                                                                   cATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGt                               <  1:476214/100‑1 (MQ=255)
                                                                    aTACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTa                              >  2:1850673/1‑100 (MQ=255)
                                                                    aTACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTa                              >  1:1576819/1‑100 (MQ=255)
                                                                     tACTTGGCCAGGACCCGGATGTCGTCATGGTGGGGGGACATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTAc                             <  2:970145/100‑1 (MQ=255)
                                                                     tACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACAGGCCATCTGGTAc                             >  2:1965880/1‑100 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTagat                           >  1:14891/1‑97 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                           >  1:1964481/1‑100 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                           >  1:2157647/1‑100 (MQ=255)
                                                                        ttCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActct                          >  2:665485/1‑100 (MQ=255)
                                                                        ttCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActct                          <  1:1509222/100‑1 (MQ=255)
                                                                         tCGCCAGGACCCGGATGTCGTCATGGGGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGGGGGTCAGGCCTCGCTCACCGGCCATCTGGTACtctc                         >  2:2252605/1‑100 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGt                               <  2:14891/93‑1 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACtctc                         <  2:2211138/99‑1 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                        >  2:1801685/1‑100 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                        >  1:3597743/1‑100 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                        <  1:1717535/100‑1 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                        <  1:2278572/100‑1 (MQ=255)
                                                                            ccAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCaa                       <  2:1363302/99‑1 (MQ=255)
                                                                             cAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAaca                     <  1:769834/100‑1 (MQ=255)
                                                                                gACACGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTc                  <  2:394247/100‑1 (MQ=255)
                                                                                  cccGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCca                >  2:455368/1‑100 (MQ=255)
                                                                                   tccgATCTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCcac               <  1:757485/97‑1 (MQ=255)
                                                                                       aTGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTa           >  1:566420/1‑100 (MQ=255)
                                                                                         gtcgtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAAc         >  2:882748/1‑100 (MQ=255)
                                                                                          tcgtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACa        >  2:1091960/1‑100 (MQ=255)
                                                                                          tcgtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACa        >  2:146297/1‑100 (MQ=255)
                                                                                           cgtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAg       <  1:3416775/100‑1 (MQ=255)
                                                                                           cgtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAg       >  2:1819886/1‑100 (MQ=255)
                                                                                            gtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGt      >  2:3716158/1‑100 (MQ=255)
                                                                                              cATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGc    <  1:2373619/100‑1 (MQ=255)
                                                                                                tggtggGGGCAAATTCGTGATACAGAAACCGATCAGATTGCGGTTCAGGCCTCGCTCACCGGCACTCTCGTACTCTCAACACTCCCCACTAACAGTGCAt  >  2:3491954/1‑100 (MQ=255)
                                                                                                     |                                                                                                
TGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCAT  >  NC_000913/3459189‑3459383

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