Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 966,368 A→G 69.7% 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,3680AG69.7% 194.9 / 98.7 136A683A (GCA→GCGycaIComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (18/23);  new base G (46/48);  total (64/72)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 137 total aligned reads displayed.

CTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCA  >  NC_000913/966269‑966464
                                                                                                   |                                                                                                
cTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCa                                                                                                  >  1:2412509/1‑100 (MQ=255)
  ggCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAGTGGAAAAGTGGCACt                                                                                                >  2:3343966/1‑100 (MQ=255)
    cagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTc                                                                                               >  1:419069/1‑99 (MQ=255)
    cagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCg                                                                                              <  1:1975360/100‑1 (MQ=255)
    cagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCg                                                                                              <  2:476303/100‑1 (MQ=255)
     agcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgc                                                                                             >  2:3190809/1‑100 (MQ=255)
     agcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgc                                                                                             <  1:1903560/100‑1 (MQ=255)
      gcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgca                                                                                            >  2:1758169/1‑100 (MQ=255)
        aaGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCgcatc                                                                                          <  1:1549268/100‑1 (MQ=255)
         aGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcg                                                                                         >  2:797563/1‑100 (MQ=255)
          gTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcgg                                                                                        >  1:1054756/1‑100 (MQ=255)
          gTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcgg                                                                                        >  1:338300/1‑100 (MQ=255)
          gTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCgcatcgg                                                                                        <  1:1694913/100‑1 (MQ=255)
           tGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggc                                                                                       <  1:2483910/100‑1 (MQ=255)
             naGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCgcatcggcat                                                                                     >  1:1782990/2‑100 (MQ=255)
              aGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCgcatcggcatc                                                                                    >  1:2370583/1‑100 (MQ=255)
                gCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATcgc                                                                                  <  1:2625210/100‑1 (MQ=255)
                 cAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATcgcg                                                                                 >  2:1740498/1‑100 (MQ=255)
                 cAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATcgcg                                                                                 <  1:2781364/100‑1 (MQ=255)
                  aaCATTGCTTCAGGTACCTCCCCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCgcgc                                                                                <  2:2516831/100‑1 (MQ=255)
                  aaCATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCgcgc                                                                                <  2:1964011/100‑1 (MQ=255)
                  aaCATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCgcgc                                                                                <  2:1117966/100‑1 (MQ=255)
                  aaCATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGAATCGGCATcgcg                                                                                 >  1:23057/1‑99 (MQ=255)
                   aCATTGCTTCAGGTACCTAACCATGACAGTAATACCTAATCATCGTTGCCATTAATTCAGCGAGTGAATCAGAAAATGGCGCTCGCATCGGCATCgagc                                                                                >  1:919927/1‑99 (MQ=255)
                    cATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCta                                                                              <  1:560834/100‑1 (MQ=255)
                        gCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGTATCGGCATCGCTCTATAAc                                                                          <  1:2826806/100‑1 (MQ=255)
                        gCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAAc                                                                          >  2:847552/1‑100 (MQ=255)
                         cTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACg                                                                         <  2:13050/100‑1 (MQ=255)
                         cTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACg                                                                         <  2:1770063/100‑1 (MQ=255)
                          ttCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGc                                                                        <  1:2458671/100‑1 (MQ=255)
                            cAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCAt                                                                      >  2:2464976/1‑100 (MQ=255)
                               gTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGc                                                                   <  1:1901601/100‑1 (MQ=255)
                               gTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGc                                                                   <  1:1488625/100‑1 (MQ=255)
                                tACCTCACCATGGCAGTAATACCTCATCGTCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCg                                                                  >  1:3166515/1‑100 (MQ=255)
                                tACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCg                                                                  >  2:2919386/1‑100 (MQ=255)
                                 aCCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGa                                                                 >  2:804084/1‑100 (MQ=255)
                                  ccTCACCATGGCCGTAATACCTCATCATCGTTGCCATTACTTCAGCGAGTGAATGGAAAACTGGCGCTCGCATCGGCATCGCGCTATAACGTATGGCGAc                                                                >  1:1405682/1‑100 (MQ=255)
                                  ccTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGAc                                                                >  1:3029909/1‑100 (MQ=255)
                                    tCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTg                                                              >  2:3246283/1‑100 (MQ=255)
                                    tCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTg                                                              >  2:2612439/1‑100 (MQ=255)
                                      aCCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGcc                                                            >  1:949182/1‑100 (MQ=255)
                                       ccATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGccc                                                           >  2:3325100/1‑100 (MQ=255)
                                       ccATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGccc                                                           <  2:1712606/100‑1 (MQ=255)
                                        cATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGc                                                             <  1:1618135/97‑1 (MQ=255)
                                         aTGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctn                                                         <  1:3108509/100‑2 (MQ=255)
                                          tGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctct                                                        <  1:1216278/100‑1 (MQ=255)
                                          tGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctct                                                        <  2:866992/100‑1 (MQ=255)
                                           ggCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTa                                                       >  2:918943/1‑100 (MQ=255)
                                           ggCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTa                                                       <  1:505220/100‑1 (MQ=255)
                                              aGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACa                                                    >  1:188038/1‑100 (MQ=255)
                                               gTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACaa                                                   >  2:2783658/1‑100 (MQ=255)
                                                 aaTACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAg                                                 <  1:3253405/100‑1 (MQ=255)
                                                  aTACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGt                                                <  2:3115682/100‑1 (MQ=255)
                                                    aCCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGTGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTa                                              >  2:2374742/1‑100 (MQ=255)
                                                     ccTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTaa                                             <  2:3065639/100‑1 (MQ=255)
                                                     ccTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTCACAAAGTTaa                                             >  1:1150881/1‑100 (MQ=255)
                                                      cTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAg                                            <  2:1660289/100‑1 (MQ=255)
                                                       tcatcaTCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGc                                           <  2:549093/100‑1 (MQ=255)
                                                        catcatCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCa                                          >  2:1485740/1‑100 (MQ=255)
                                                           catcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCTCTCTAACAAAGTTAAGCATCg                                       <  1:2953599/100‑1 (MQ=255)
                                                           catcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCg                                       >  2:1635679/1‑100 (MQ=255)
                                                            atcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGc                                      <  2:1413758/100‑1 (MQ=255)
                                                            atcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGc                                      <  1:406710/100‑1 (MQ=255)
                                                            atcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGc                                      >  1:867312/1‑100 (MQ=255)
                                                            atcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGc                                      <  2:1011414/100‑1 (MQ=255)
                                                               gTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTAt                                   <  2:2112707/100‑1 (MQ=255)
                                                               gTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTAt                                   <  1:2116819/100‑1 (MQ=255)
                                                                  gCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCaa                                >  2:337209/1‑100 (MQ=255)
                                                                    cATTAATTCAGCGAGTGAATGGAAAAGTGTCACTCGCATCGGCATCGCGCTATAACGCATTGCGACTGCCCTCTAACAAATTTAAGCATCTCTATCAACt                              >  2:1339947/1‑100 (MQ=255)
                                                                    cATTAAATCAGCGAGGGAATGGAAGAGTGGCGCTGGCATCGGCATCGCGCTATAACACATGGCGACTGCCCTCTAACAACGTTAAGCATCGCTATGAACt                              >  1:596009/1‑100 (MQ=255)
                                                                      ttAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGc                            <  2:2172152/100‑1 (MQ=255)
                                                                      ttAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGc                            <  2:1174261/100‑1 (MQ=255)
                                                                        aaTTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATTGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCaa                          >  2:2988315/1‑100 (MQ=255)
                                                                         aTTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAg                         <  2:3389440/100‑1 (MQ=255)
                                                                         aTTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCACGGCGACTGCCCTCTAACAAAGTTCAGCATCGCTATCAACTGCa                           >  1:3174155/1‑98 (MQ=255)
                                                                          ttCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAgg                        >  1:2455336/1‑100 (MQ=255)
                                                                          ttCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAgg                        <  1:1368406/100‑1 (MQ=255)
                                                                           tCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGa                       <  1:739246/100‑1 (MQ=255)
                                                                           tCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGa                       <  2:2853534/100‑1 (MQ=255)
                                                                              gCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGatat                    >  2:1081988/1‑100 (MQ=255)
                                                                              gCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGatat                    <  1:2708798/100‑1 (MQ=255)
                                                                               cGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATc                   >  2:221575/1‑100 (MQ=255)
                                                                                 aGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCaa                 <  2:742900/100‑1 (MQ=255)
                                                                                  gTGAATGGAAAAGTGGCCCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATc                   <  1:1640876/97‑1 (MQ=255)
                                                                                   tGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATg               <  2:2106703/100‑1 (MQ=255)
                                                                                   tGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATg               <  2:1875904/100‑1 (MQ=255)
                                                                                      aTGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAt            >  2:1335623/1‑100 (MQ=255)
                                                                                      aTGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGa             <  2:493565/99‑1 (MQ=255)
                                                                                       tGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAtt           <  1:1117659/100‑1 (MQ=255)
                                                                                        ggAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTg          >  1:489990/1‑100 (MQ=255)
                                                                                        ggAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTg          >  1:1766646/1‑100 (MQ=255)
                                                                                        ggAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTg          >  1:2676999/1‑100 (MQ=255)
                                                                                         gAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGa         >  1:727598/1‑100 (MQ=255)
                                                                                         gAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGa         <  1:1668049/100‑1 (MQ=255)
                                                                                           aaaGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATa       <  2:1695980/100‑1 (MQ=255)
                                                                                               tGGCGCTCGCATTGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTcc   <  1:2988289/100‑1 (MQ=255)
                                                                                               tGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATCCATTGATatcct   <  2:2993613/100‑5 (MQ=255)
                                                                                                ggCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCca  >  2:3192957/1‑100 (MQ=255)
                                                                                                ggCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCca  <  2:2667128/100‑1 (MQ=255)
                                                                                                ggCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCca  >  2:1767167/1‑100 (MQ=255)
                                                                                                   |                                                                                                
CTGGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCA  >  NC_000913/966269‑966464

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