Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,602,188 T→C 78.3% I233I (ATT→ATC lsrA → autoinducer 2 import ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,602,1880TC78.3% 202.3 / 41.7 104I233I (ATT→ATClsrAautoinducer 2 import ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (11/11);  new base C (34/45);  total (47/57)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

GCCGGAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAAC  >  NC_000913/1602089‑1602286
                                                                                                   |                                                                                                  
 ccGGAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTc                                                                                                   <  2:3011276/100‑1 (MQ=255)
   ggAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAg                                                                                                 <  2:1804817/100‑1 (MQ=255)
      aactCGCCAGATTGCCGATCGAATTAGCGTGAGGCGCGACGGAACCATCGGCTTACGCGGTTCAGCCAGCGAACTGTCTATCGATGACATTATGCAGGcc                                                                                              <  2:593503/97‑1 (MQ=255)
      aaTTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGcc                                                                                              <  1:650339/100‑1 (MQ=255)
          cgccacacagCCGATCGAATTAGCGTGATGAGCGACGGAACCATCGCATTAATCGGCAAAACCAGCGAACTGTCTACCGACGACATTATGCAGGCCAt                                                                                            >  2:1284997/10‑98 (MQ=255)
          cGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCa                                                                                          <  1:2282529/100‑1 (MQ=255)
          cGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCa                                                                                          <  1:410092/100‑1 (MQ=255)
          cGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCa                                                                                          >  1:2076749/1‑100 (MQ=255)
           gCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCAc                                                                                         <  2:3168853/100‑1 (MQ=255)
           gCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCAc                                                                                         <  1:2598654/100‑1 (MQ=255)
            ccAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCAcc                                                                                        <  1:812536/100‑1 (MQ=255)
             cAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCAccc                                                                                       <  1:2219108/100‑1 (MQ=255)
              aGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCAcccc                                                                                      <  2:2016643/100‑1 (MQ=255)
              aGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCAcccc                                                                                      >  2:21440/1‑100 (MQ=255)
               gATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAACCCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCa                                                                                     <  2:224061/100‑1 (MQ=255)
                 ttGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGc                                                                                   <  2:2149207/100‑1 (MQ=255)
                 ttGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCACCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAg                                                                                    >  2:897782/1‑99 (MQ=255)
                  tGCCGATTGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCg                                                                                  >  2:1646320/1‑100 (MQ=255)
                  tGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCg                                                                                  <  1:2080563/100‑1 (MQ=255)
                  tGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCg                                                                                  >  1:2080996/1‑100 (MQ=255)
                   gCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCg                                                                                  >  2:2891221/1‑99 (MQ=255)
                   gCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCgg                                                                                 >  2:918128/1‑100 (MQ=255)
                   gCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCgg                                                                                 <  2:1882020/100‑1 (MQ=255)
                    ccGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGt                                                                                >  1:3837984/1‑100 (MQ=255)
                     cGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTa                                                                               >  2:466079/1‑100 (MQ=255)
                      gctCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAACGCATCGAACTGTCTACCGACGACATTATTCACGCCATCACCCCAGCGGt                                                                                <  2:2834487/96‑1 (MQ=255)
                      gATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTAc                                                                              >  1:1889437/1‑100 (MQ=255)
                           aaTTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGaa                                                                         >  2:4001476/1‑100 (MQ=255)
                           aaTTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGaa                                                                         <  1:3246082/100‑1 (MQ=255)
                            aTTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGaaa                                                                        <  1:2154940/100‑1 (MQ=255)
                             ttAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGaaaa                                                                       <  2:3419621/100‑1 (MQ=255)
                              tAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGaaaaa                                                                      >  1:3363967/1‑100 (MQ=255)
                               aGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAAt                                                                     >  1:2305304/1‑100 (MQ=255)
                                gCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCTCCCCAGCGGTACGGGAAAAATc                                                                    <  2:2418501/100‑1 (MQ=255)
                                 cGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCg                                                                   >  1:1309398/1‑100 (MQ=255)
                                  gTGCTGCGCGACGGAACCATCGCCTTAAGCGGCAAACCCAGCGAACTGTCGACCGACGACATTATCCGGGCCATCACCCCAGCGGTACGGGAAAAATCGc                                                                  <  1:1191736/100‑1 (MQ=255)
                                  gTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGc                                                                  >  2:2884323/1‑100 (MQ=255)
                                   tGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACAGGAAAAAtacct                                                                 >  2:2535430/1‑96 (MQ=255)
                                    gATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGctc                                                                <  1:881678/100‑1 (MQ=255)
                                      tGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGctctc                                                              <  2:921165/100‑1 (MQ=255)
                                        cgtgACGGAACCATCGCCGTAATCGGCAAAACTAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGTAAAAATCGctctct                                                             >  2:1869216/4‑99 (MQ=255)
                                            aCGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTgccag                                                        <  1:1504464/100‑1 (MQ=255)
                                                aaCCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaa                                                    >  1:406472/1‑100 (MQ=255)
                                                 aCCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaaa                                                   <  1:1808085/100‑1 (MQ=255)
                                                    aTCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAAt                                                >  1:1508141/1‑100 (MQ=255)
                                                    aTCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAAt                                                >  1:2893231/1‑100 (MQ=255)
                                                      cGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTa                                              <  1:2646164/100‑1 (MQ=255)
                                                       gCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCCAAAATTAt                                             >  1:2256994/1‑100 (MQ=255)
                                                       gCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTAt                                             <  2:3134500/100‑1 (MQ=255)
                                                       gCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTAt                                             >  2:3542502/1‑100 (MQ=255)
                                                        ccTTAAGCGGCAAAACCAGCGAACTGTCTACCGAAGACATTATCCAGGTCATCACCCAAACAATAAGAGAAAAATCGCTCTCTGCCAGCCAAAAATTAt                                             <  2:3549709/99‑1 (MQ=255)
                                                          ttAAGCGGCAAAGCCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCAAGCCAAAAATTATGGc                                          <  2:1674180/100‑1 (MQ=255)
                                                          ttAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGc                                          <  2:3233930/100‑1 (MQ=255)
                                                          ttAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGc                                          <  1:2792498/100‑1 (MQ=255)
                                                           tAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCt                                         <  2:972313/100‑1 (MQ=255)
                                                            aaGCGGCAAAACCAGCGATCTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTg                                        <  2:2946042/100‑1 (MQ=255)
                                                            aaGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTg                                        >  1:1554062/1‑100 (MQ=255)
                                                             aGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTgg                                       >  2:2041083/1‑100 (MQ=255)
                                                             aGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTgg                                       <  1:2570816/100‑1 (MQ=255)
                                                             aGCGGAAAAACCAGCGAACTGTCTAACGAAGACATTATTCAGGCCATCACCCCAGAAAAAAAGAAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCt                                         <  2:972471/98‑1 (MQ=255)
                                                              gCGGTAAAACCCGCGAACGGTGTACCGACGACAGTATCCAGGCCAGCGCCCCAGCGGTACGGGAAAAATCGCTCTCGGCCAGCCATAATTTTTGGCTggc                                      <  2:1632650/100‑2 (MQ=255)
                                                              gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCGGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                      <  1:963528/100‑2 (MQ=255)
                                                              gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                      <  1:1582448/100‑2 (MQ=255)
                                                              gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCCCTCTGCCAGCCAAAAATTATGGCTGGa                                      >  1:1155732/1‑100 (MQ=255)
                                                               cGGCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAg                                     <  1:4579/100‑1 (MQ=255)
                                                                ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGt                                    >  2:1613823/1‑100 (MQ=255)
                                                                ggCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGa                                      >  2:3729419/1‑98 (MQ=255)
                                                                 gCAAAACCAGCGAACTGTCTACCGACGACATTATCCATGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGtt                                   <  2:1041891/100‑1 (MQ=255)
                                                                  cAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcagag                                  <  1:2844728/100‑6 (MQ=38)
                                                                  cAAAACCAGCGAACTGTCTACCGACGACATGATGCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATCGCTGGAGTTa                                  >  2:2740111/1‑100 (MQ=255)
                                                                    aaaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTAcc                                >  1:1589743/1‑100 (MQ=255)
                                                                     aaCCAGCGAACTGTCTCCCGACGACATTATCCAGCCCATCACCCCAGCGGTACCGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTAcc                                <  1:2921017/99‑1 (MQ=255)
                                                                     aaCCAGCGAACTGTCTACCGACGACATTATCTAGGCCATCACCCCAGCGGTACGGGAATAATCGCTCACTGCCAGCCAAAAATTATGGCTGGAGTTAACt                               >  2:2888257/1‑100 (MQ=255)
                                                                     aaCCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTAcc                                <  2:1836482/99‑1 (MQ=255)
                                                                     aaCCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCt                               <  2:1601020/100‑1 (MQ=255)
                                                                       ccAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTgg                             >  2:306290/1‑100 (MQ=255)
                                                                       ccAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTg                              <  1:1001453/99‑1 (MQ=255)
                                                                        cAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGt                            >  1:2009227/1‑100 (MQ=255)
                                                                        cAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGt                            <  1:2093520/100‑1 (MQ=255)
                                                                          gCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAACAATTATGGCTGGAGTTACCTGGTaa                          >  2:1549695/1‑100 (MQ=255)
                                                                          gCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTaa                          >  1:1373894/1‑100 (MQ=255)
                                                                            gAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAAcc                        >  2:2204297/1‑100 (MQ=255)
                                                                            gAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAAcc                        <  2:3110665/100‑1 (MQ=255)
                                                                            gAACTGTCTACCGACGACATTATCAAGGCCATAAACACAGGGGTACGGGAAAAATCGAAAAAAAAAAAAAAAAAATTATGGCTGGAGTTACCTGGTaa                          <  2:2037209/98‑1 (MQ=255)
                                                                               cTGTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGcc                     >  1:950460/1‑100 (MQ=255)
                                                                                 gTCTACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCca                   <  1:306293/100‑1 (MQ=255)
                                                                                    tACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaa                <  2:1300896/100‑1 (MQ=255)
                                                                                    tACCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaa                >  1:254902/1‑100 (MQ=255)
                                                                                     aCCGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaac               <  1:1732549/100‑1 (MQ=255)
                                                                                      ccGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaaca              <  1:2041093/100‑1 (MQ=255)
                                                                                      ccGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACAGCCCacaac               <  2:3095613/99‑1 (MQ=255)
                                                                                      ccGACGACATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaaca              >  2:2254804/1‑100 (MQ=255)
                                                                                        gacgacATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATg            <  1:3209899/100‑1 (MQ=255)
                                                                                        gacgacATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATg            >  2:833109/1‑100 (MQ=255)
                                                                                        gacgacATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATg            >  1:1653878/1‑100 (MQ=255)
                                                                                        gacgacATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCGCTCAGCAAGCCAAAAAGTATGGCTGGAGTTACCTGGTACCCGCCAAAAACAt             >  2:3201601/1‑99 (MQ=255)
                                                                                          cgacATTATCCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgcc          >  1:322274/1‑100 (MQ=255)
                                                                                            acATTATTCAGACCCTCACCACACGGGTACAACAAACTTCGCTCTCTCACAGCCAATAATTATGGCTGGAGTTATCTGTTAACCGCCCAAAACATgccgc        <  2:661145/100‑1 (MQ=255)
                                                                                             nattatcCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccgcc       >  1:1119107/2‑100 (MQ=255)
                                                                                                  ttCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACGGCCCACAACATGCCGCCGGAAc  >  2:1760658/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GCCGGAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAAC  >  NC_000913/1602089‑1602286

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