Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,459,293:1 (G)5→6 100% 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.G92.4% 383.6 / 24.8 119coding (950/1482 nt)gspEgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base . (4/5);  new base G (47/61);  total (52/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

GGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGC  >  NC_000913/3459190‑3459381
                                                                                                     |                                                                                              
ggAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtggtg                                                                                               <  2:519696/100‑1 (MQ=255)
 gAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtggtgg                                                                                              <  2:2213596/100‑1 (MQ=255)
   aGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggg                                                                                            <  1:2781824/100‑1 (MQ=255)
    ggCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTggggg                                                                                          <  1:1919987/100‑1 (MQ=255)
     gCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggg                                                                                            <  1:792668/98‑1 (MQ=255)
      cATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggggga                                                                                        <  2:1661033/100‑3 (MQ=255)
      cATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggggga                                                                                        >  1:2284526/1‑98 (MQ=255)
      cATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggggga                                                                                        >  2:2859998/1‑98 (MQ=255)
          gggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAAtt                                                                                     <  2:2304380/100‑1 (MQ=255)
          gggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTc                                                                                    >  2:1813208/1‑100 (MQ=255)
           ggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGCTGCCGTCATGGGGGGGGgaaa                                                                                      >  1:2476957/1‑96 (MQ=255)
           ggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTc                                                                                    >  2:1832688/1‑100 (MQ=255)
            gCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGt                                                                                  <  2:397259/100‑1 (MQ=255)
             cAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGt                                                                                  <  1:1126150/100‑1 (MQ=255)
             cAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGt                                                                                  <  2:291924/100‑1 (MQ=255)
              aaaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTg                                                                                 <  2:2563323/100‑1 (MQ=255)
              aaaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTg                                                                                 <  1:2207558/100‑1 (MQ=255)
              aaaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTg                                                                                 <  2:1472884/100‑1 (MQ=255)
              aaaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTg                                                                                 >  2:568762/1‑100 (MQ=255)
               aaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCTCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGa                                                                                <  1:40410/100‑1 (MQ=255)
               aaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGa                                                                                >  1:1738496/1‑100 (MQ=255)
               aaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGa                                                                                >  1:1737345/1‑100 (MQ=255)
                 cGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATa                                                                              >  1:2075501/1‑100 (MQ=255)
                 cGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATa                                                                              >  2:575514/1‑100 (MQ=255)
                  gCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATAc                                                                             >  2:1635748/1‑100 (MQ=255)
                   cAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACa                                                                            <  1:534546/100‑1 (MQ=255)
                   cAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACa                                                                            >  1:1679786/1‑100 (MQ=255)
                    aGGTTAATACCCGTGTGGATATGTCTTTCGCCCGTGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATAc                                                                             <  1:403040/98‑1 (MQ=255)
                    aGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGAGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATAc                                                                             <  1:2194126/98‑1 (MQ=255)
                     ggTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGa                                                                          >  1:2807313/1‑100 (MQ=255)
                      gTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGaa                                                                         <  2:2208699/100‑1 (MQ=255)
                      gTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTCGGGGGAAATTCGTGATACAGaa                                                                         <  1:1690592/100‑1 (MQ=255)
                        gAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAAc                                                                       >  2:187877/1‑100 (MQ=255)
                        gAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAAc                                                                       <  2:188815/100‑1 (MQ=255)
                          aTACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACcg                                                                     >  2:776586/1‑100 (MQ=255)
                           tACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACcgc                                                                    <  2:940720/100‑1 (MQ=255)
                            aCCCGTGTGGATATGTCTTTCGCTCGCGGGCTGCGCGCCATACTTCGCCGGGACCCGGATGTCGACATCGTAGGGGGAGATTCGTGATACAGAAACagcg                                                                   >  1:2305273/1‑100 (MQ=255)
                             cccGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCgcgc                                                                  >  1:1159470/1‑100 (MQ=255)
                             cccGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCgcgc                                                                  >  2:186600/1‑100 (MQ=255)
                              ccGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCa                                                                 <  2:510531/100‑1 (MQ=255)
                               cGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATAGTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCa                                                                 <  1:1870763/99‑1 (MQ=255)
                               cGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAg                                                                >  2:2500458/1‑100 (MQ=255)
                               cGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAg                                                                <  2:113152/100‑1 (MQ=255)
                                   tgGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTg                                                            <  2:1760389/100‑1 (MQ=255)
                                    ggATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGc                                                           <  2:2284536/100‑1 (MQ=255)
                                    ggATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGc                                                           <  1:1369451/100‑1 (MQ=255)
                                    ggATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGc                                                           >  1:2311215/1‑100 (MQ=255)
                                      atatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCgg                                                         <  2:2340209/100‑1 (MQ=255)
                                      atatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCgg                                                         <  2:2807323/100‑1 (MQ=255)
                                      atatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCgg                                                         <  2:2075508/100‑1 (MQ=255)
                                      atatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGt                                                        <  2:1924386/100‑1 (MQ=255)
                                         tGTCTTTCGCTCGCGGCCTGCGCGCCATAGTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTc                                                      <  1:1088787/100‑1 (MQ=255)
                                          gTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGCTGTCGTCATGGTGGGGGGAAATTCCTGATACAGAAACCGCGCAGATTGCGGTTCa                                                     >  1:829412/1‑100 (MQ=255)
                                          gTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCa                                                     <  1:2769484/100‑1 (MQ=255)
                                          gTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAg                                                    <  1:149377/100‑1 (MQ=255)
                                           tCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAg                                                    <  2:2765709/100‑1 (MQ=255)
                                            cTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgg                                                   >  2:151658/1‑100 (MQ=255)
                                            cTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgg                                                   >  1:1982363/1‑100 (MQ=255)
                                             tttCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTAGGCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGc                                                  <  2:1091686/100‑1 (MQ=255)
                                                cgctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTc                                               >  1:255814/1‑100 (MQ=255)
                                                   tcgcGGCCTGCGCGCCATACTTACCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTc                                           <  1:1302919/100‑1 (MQ=255)
                                                     gcgGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCa                                          >  2:2503/1‑100 (MQ=255)
                                                       ggCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCAcc                                        >  2:1381749/1‑100 (MQ=255)
                                                        gCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCg                                       >  2:169952/1‑100 (MQ=255)
                                                           tGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGcc                                    <  2:952677/100‑1 (MQ=255)
                                                           tGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGcc                                    <  1:2803946/100‑1 (MQ=255)
                                                             cgcgcCATACTTCGCCAGGACCCGGATGTCGTCATGGGGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCa                                   <  1:1582038/99‑1 (MQ=255)
                                                               cgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCt                                >  2:2601906/1‑100 (MQ=255)
                                                                 ccATACTTCGGCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTgg                              <  1:2356089/100‑1 (MQ=255)
                                                                   aTACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTa                            <  2:320193/100‑1 (MQ=255)
                                                                   aTACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTa                            <  2:136719/100‑1 (MQ=255)
                                                                    tACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTAc                           <  2:255816/100‑1 (MQ=255)
                                                                    tACTTCGCCAGGACCCGGATGTCGCCATGGTTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGGTCAGACCTCGCTCACCGGCCATCTGGTAc                           >  1:221752/1‑100 (MQ=255)
                                                                      cTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                         <  1:516065/100‑1 (MQ=255)
                                                                      cTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                         >  1:1901584/1‑100 (MQ=255)
                                                                       ttCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActct                        <  2:154611/100‑1 (MQ=255)
                                                                          gCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCaa                     >  2:2319611/1‑100 (MQ=255)
                                                                           ccAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCaan                    <  2:808232/100‑2 (MQ=255)
                                                                           ccAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                    >  1:456156/1‑100 (MQ=255)
                                                                           ccAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                    <  2:281710/100‑1 (MQ=255)
                                                                           ccAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                    <  2:935957/100‑1 (MQ=255)
                                                                           ccAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                    <  2:1906923/100‑1 (MQ=255)
                                                                              ggacgtggATGTTGTCATGGTGGGGTGAGATTCGTGATACAGAAACCGTGTAGATTTCGGTTCAGCCCTCGCTCACCGGCCATCTGGTACTCTCAacac                  <  1:1745485/93‑1 (MQ=255)
                                                                              ggACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGTGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACt                 <  1:745279/100‑1 (MQ=255)
                                                                                aCCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTcc               <  2:1240657/100‑1 (MQ=255)
                                                                                  ccGGGTGTCGTCATGGTGGGGGGAAAGTCGTGATACAGGACGCGCGCAGATTGCGGTTCAGGCCTCGCTCGCCGGCCATCTGGTACTCTCAACACTCcac             <  1:1204805/100‑1 (MQ=255)
                                                                                  ccGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCcac             >  1:613017/1‑100 (MQ=255)
                                                                                     gATGTCGTCAGGGTGGGGGGAAATGCGTGATACCGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACt          <  1:2324382/100‑1 (MQ=255)
                                                                                      aTGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTa         <  2:2851687/100‑1 (MQ=255)
                                                                                       tGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTaa        >  2:2416951/1‑100 (MQ=255)
                                                                                        gtcgtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAAc       <  2:646262/100‑1 (MQ=255)
                                                                                         tcgtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACa      >  1:2338386/1‑100 (MQ=255)
                                                                                          cgtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAg     >  1:313081/1‑100 (MQ=255)
                                                                                          cgtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAg     >  2:1701991/1‑100 (MQ=255)
                                                                                          cgtcATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGt    >  2:475739/1‑100 (MQ=255)
                                                                                           gtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGt    >  1:2101003/1‑100 (MQ=255)
                                                                                            tcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTg   >  1:2176835/1‑100 (MQ=255)
                                                                                            tcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTg   <  1:2503/100‑1 (MQ=255)
                                                                                             cATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGc  >  2:2733971/1‑100 (MQ=255)
                                                                                              aTGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGc  <  1:1684481/99‑1 (MQ=255)
                                                                                                     |                                                                                              
GGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGC  >  NC_000913/3459190‑3459381

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