Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 966,368 A→G 100% A683A (GCA→GCG ycaI → ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913966,3680AG98.4% 422.7 / ‑3.4 125A683A (GCA→GCGycaIComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (64/59);  total (65/60)
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 125 total aligned reads displayed.

GCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCA  >  NC_000913/966272‑966464
                                                                                                |                                                                                                
gcagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTc                                                                                               >  2:702607/1‑100 (MQ=255)
gcagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTc                                                                                               >  2:1470884/1‑100 (MQ=255)
 cagcaAGTACAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCg                                                                                              <  2:110399/100‑1 (MQ=255)
     aaGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatc                                                                                          >  2:1073751/1‑100 (MQ=255)
      aGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcg                                                                                         <  1:421389/100‑1 (MQ=255)
       gTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcgg                                                                                        <  1:2099745/100‑1 (MQ=255)
       gTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATAATCGTTGCAATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcgg                                                                                        >  2:665832/1‑100 (MQ=255)
         gCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggca                                                                                      >  1:1279800/1‑100 (MQ=255)
          cAGGCAACATTGCTTCAGGTAACAAACAAAGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAATGAATCGAAAAGTGGCGCTccgttaattaa                                                                                     >  2:1300642/1‑90 (MQ=255)
            ggCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggcatcg                                                                                   >  2:2850109/1‑100 (MQ=255)
            ggCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggcatcg                                                                                   >  1:2785929/1‑100 (MQ=255)
              cAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGACAAGTGGCGATCGCATCGGCATcgc                                                                                  >  2:693203/1‑99 (MQ=255)
               aaCATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCgcgc                                                                                >  1:2390205/1‑100 (MQ=255)
                 cATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCta                                                                              <  1:2828717/100‑1 (MQ=255)
                 cATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCta                                                                              >  1:1496942/1‑100 (MQ=255)
                  aTTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCtat                                                                             <  2:1393263/100‑1 (MQ=255)
                  aTTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCtat                                                                             <  1:923746/100‑1 (MQ=255)
                   ttGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGGAAAGTGGCGCTCGCATCGGCATCGCGCtata                                                                            >  2:2046179/1‑100 (MQ=255)
                     gCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAAc                                                                          >  2:2535627/1‑100 (MQ=255)
                      cTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACg                                                                         >  1:300407/1‑100 (MQ=255)
                      cTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACg                                                                         <  1:387323/100‑1 (MQ=255)
                        tCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCa                                                                       >  1:1939173/1‑100 (MQ=255)
                           ggTGCCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGAATCGCGCTATAACGCATgg                                                                    <  2:1113630/100‑1 (MQ=255)
                           ggTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATgg                                                                    >  1:2113715/1‑100 (MQ=255)
                           ggTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATgg                                                                    >  2:1837794/1‑100 (MQ=255)
                            gTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGc                                                                   <  2:1334205/100‑1 (MQ=255)
                             tACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCg                                                                  <  2:2855862/100‑1 (MQ=255)
                             tACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCg                                                                  <  1:2022584/100‑1 (MQ=255)
                               ccTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGAc                                                                <  1:1470899/100‑1 (MQ=255)
                               ccTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGAc                                                                <  2:2094996/100‑1 (MQ=255)
                               ccTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGAc                                                                >  2:308990/1‑100 (MQ=255)
                               ccTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGAc                                                                <  2:1908625/100‑1 (MQ=255)
                               ccTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGAc                                                                <  1:1793918/100‑1 (MQ=255)
                                cTCACCATTGCAGTAATACCTCATCATCGATTCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCTCATCGGCATCGCGCTATAACGCATGGCGAc                                                                <  2:2270501/99‑1 (MQ=255)
                                cTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACt                                                               >  1:718883/1‑100 (MQ=255)
                                 tCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTg                                                              >  1:989143/1‑100 (MQ=255)
                                 tCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTg                                                              <  2:1686227/100‑1 (MQ=255)
                                 tCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTg                                                              <  1:1639176/100‑1 (MQ=255)
                                  cACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAATGCATGGCGACttc                                                             <  1:2046202/100‑3 (MQ=255)
                                  cACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGc                                                             <  2:483098/100‑1 (MQ=255)
                                   aCCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGcc                                                            <  1:968357/100‑1 (MQ=255)
                                    ccATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGccc                                                           <  1:1010736/100‑1 (MQ=255)
                                    ccATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGccc                                                           <  1:271301/100‑1 (MQ=255)
                                    ccATGGCAGTAATACCTCATCATCGTTGCCATGAATGCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGccc                                                           <  2:331955/100‑1 (MQ=255)
                                     cATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCct                                                          <  2:36993/100‑1 (MQ=255)
                                     cATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCct                                                          >  2:539391/1‑100 (MQ=255)
                                      aTGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctn                                                         <  2:27908/100‑2 (MQ=255)
                                      aTGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctc                                                         <  2:690233/100‑1 (MQ=255)
                                      aTGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctc                                                         <  2:2257596/100‑1 (MQ=255)
                                      aTGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctc                                                         <  1:2233138/100‑1 (MQ=255)
                                        ggCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTa                                                       >  2:2732426/1‑100 (MQ=255)
                                         gtaGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTa                                                       <  1:2535658/97‑1 (MQ=255)
                                         gCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTa                                                       <  1:117208/99‑1 (MQ=255)
                                          cAGTAATACCTCATCATCGTTGCTATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAAc                                                     <  1:925311/100‑1 (MQ=255)
                                          cAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAAc                                                     >  1:257447/1‑100 (MQ=255)
                                          cAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAAc                                                     <  1:1773580/100‑1 (MQ=255)
                                          cAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAAc                                                     >  2:1773560/1‑100 (MQ=255)
                                          cAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAAc                                                     >  1:1710144/1‑100 (MQ=255)
                                          cAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAAc                                                     >  2:2516113/1‑100 (MQ=255)
                                            gTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGTGCTCGCATCGGAATCGTGCTATAACGCATGGCGCGTGCCATCTAACa                                                    >  2:2126879/1‑99 (MQ=255)
                                             tAAGACCTCATCATCGTTGGCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGTGCTATAACTAATGGTGATTGCCCTCGAACaaa                                                  <  2:1098334/100‑1 (MQ=255)
                                              aaTACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAg                                                 <  1:428671/100‑1 (MQ=255)
                                              aaTACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAg                                                 <  1:580289/100‑1 (MQ=255)
                                              aaTACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAg                                                 <  1:1705128/100‑1 (MQ=255)
                                                tACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGtt                                               <  1:702628/100‑1 (MQ=255)
                                                  ccTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTaa                                             >  1:1477197/1‑100 (MQ=255)
                                                   cTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTaa                                             <  2:325553/99‑1 (MQ=255)
                                                   cTCATAATCTTTGTTATTAATTCTTCGAGTGAATTGAAAAGTGGAACTCGCATCGGCATCGCGCTATCACGCATGGCGACTGCCCTCTAACAAAGTTAAg                                            <  1:1796576/100‑1 (MQ=255)
                                                    tcatcaTCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGc                                           <  1:2036391/100‑1 (MQ=255)
                                                    tcatcaTCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGc                                           <  1:2370417/100‑1 (MQ=255)
                                                    tcatcaTCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGc                                           <  2:18766/100‑1 (MQ=255)
                                                     catcatCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCa                                          >  2:2485588/1‑100 (MQ=255)
                                                         atcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATCACGCATGGCGACTGCCCTCTAACAGAGTTAAGCATCGc                                      <  2:1874872/100‑1 (MQ=255)
                                                             ttGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATc                                  <  2:2261517/100‑1 (MQ=255)
                                                               gCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCaa                                >  1:2178080/1‑100 (MQ=255)
                                                                ccATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAAc                               >  1:3021516/1‑100 (MQ=255)
                                                                 cATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACt                              >  2:2530609/1‑100 (MQ=255)
                                                                 cATTAATTCAGCGAGTGAATGGAAAAATGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATGACTATCaa                                >  2:1998249/1‑98 (MQ=255)
                                                                    tAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCa                           >  1:2084767/1‑100 (MQ=255)
                                                                     aaTTCAGCGAGCGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCaa                          >  1:2988547/1‑100 (MQ=255)
                                                                       ttCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAgg                        <  1:2617745/100‑1 (MQ=255)
                                                                         cAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGat                      >  1:265688/1‑100 (MQ=255)
                                                                         cAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGat                      >  1:22487/1‑100 (MQ=255)
                                                                         cAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCACTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGat                      >  1:2448661/1‑100 (MQ=255)
                                                                          aGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGata                     >  1:1633585/1‑100 (MQ=255)
                                                                            cGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATc                   <  2:1731700/100‑1 (MQ=255)
                                                                             gAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCa                  >  1:819576/1‑100 (MQ=255)
                                                                                 gAATTGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATgg              >  1:2256636/1‑100 (MQ=255)
                                                                                 gAATTGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATgg              >  1:2419641/1‑100 (MQ=255)
                                                                                 gAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATgg              >  1:1614925/1‑100 (MQ=255)
                                                                                 gAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATgg              >  2:889101/1‑100 (MQ=255)
                                                                                   aTGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAt            >  2:1321414/1‑100 (MQ=255)
                                                                                   aTGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGAATCGCTATCAACTGCAAGGATATAAATGGAt            >  1:2085963/1‑100 (MQ=255)
                                                                                    tGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAtt           >  1:2998244/1‑100 (MQ=255)
                                                                                     ggAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTg          >  1:1287151/1‑100 (MQ=255)
                                                                                      gAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGa         >  1:1838077/1‑100 (MQ=255)
                                                                                       aaaaGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGAt        <  1:702611/100‑1 (MQ=255)
                                                                                         aaGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATAc      <  2:2081585/100‑1 (MQ=255)
                                                                                            tGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTcc   <  2:2988514/100‑1 (MQ=255)
                                                                                             ggCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCca  >  1:2740215/1‑100 (MQ=255)
                                                                                                |                                                                                                
GCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCA  >  NC_000913/966272‑966464

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: