Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,246,800 A→G 60.5% intergenic (‑31/+263) cirA ← / ← lysP colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,246,8000AG60.5% 80.6 / 121.7 114intergenic (‑31/+263)cirA/lysPcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter
Reads supporting (aligned to +/- strand):  ref base A (21/24);  new base G (38/31);  total (59/55)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

Only 100 of 127 total aligned reads displayed.

CACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAA  >  NC_000913/2246701‑2246899
                                                                                                   |                                                                                                   
cacaGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaata                                                                                                     >  1:1087474/1‑100 (MQ=255)
cacaGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaata                                                                                                     <  2:632198/100‑1 (MQ=255)
cacaGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaata                                                                                                     <  1:2953578/100‑1 (MQ=255)
   aGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaataaaa                                                                                                  >  1:3032747/1‑100 (MQ=255)
    ggCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaataaaat                                                                                                 >  1:486446/1‑100 (MQ=255)
         tGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGtt                                                                                            >  2:2100559/1‑100 (MQ=255)
          gCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTa                                                                                           <  1:3343997/100‑1 (MQ=255)
            acaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATc                                                                                         >  1:2421541/1‑100 (MQ=255)
            acaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATc                                                                                         >  2:1410189/1‑100 (MQ=255)
             caAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCg                                                                                        <  2:2151377/100‑1 (MQ=255)
               aGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCt                                                                                      <  2:1569432/100‑1 (MQ=255)
               aGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCt                                                                                      >  2:2133557/1‑100 (MQ=255)
                gAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTc                                                                                     <  2:156775/100‑1 (MQ=255)
                 aaaTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTca                                                                                    <  2:2983641/100‑1 (MQ=255)
                 aaaTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTca                                                                                    >  2:1628214/1‑100 (MQ=255)
                  aaTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGGAAATAAAATGGGTTATCGCTcac                                                                                   <  1:2669389/100‑1 (MQ=255)
                   aTAGCGGACAAACACAGCCCGAGCCCTACGAAAGGGTTCAACATAAACAGTCCAGATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTcaca                                                                                  <  1:1190897/100‑1 (MQ=255)
                     aGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATAGGCTATCGCTcaca                                                                                  >  1:3369674/1‑98 (MQ=255)
                      gCGGCCAAACACAGCCCGACCTGTACGAAAAGGTTCAACCTAAGCATTCTATATCTCCATGAGGAAACTACGAACATGAAATGATTGATCGCTTACAtct                                                                               >  1:2693031/1‑100 (MQ=255)
                       cGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACAtctt                                                                              >  2:3283156/1‑100 (MQ=255)
                         gACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtcttct                                                                            >  2:1979494/1‑100 (MQ=255)
                         gACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtcttct                                                                            <  1:454513/100‑1 (MQ=255)
                           cAAACACATCCCGACCCGTACGAAAGTGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGTGTTATCGCTCACATCttcttc                                                                          <  1:2208178/100‑1 (MQ=255)
                            aaaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCa                                                                         <  2:1944246/100‑1 (MQ=255)
                             aaCACAGCCCGACCCGTACGAAATGGTTCATCCTAAACATTCCCTATCTCCATGAGTGAACTACGAAAATAAAATGGGTTATCGCTTACATCTTCTTCAc                                                                        >  1:2019633/1‑100 (MQ=255)
                             aaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCAc                                                                        >  2:2966265/1‑100 (MQ=255)
                             aaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCAc                                                                        >  2:361835/1‑100 (MQ=255)
                              acacAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACg                                                                       <  1:701002/100‑1 (MQ=255)
                               cacaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGt                                                                      >  2:1806018/1‑100 (MQ=255)
                                 caGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTcc                                                                    <  2:137974/100‑1 (MQ=255)
                                  aGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTccc                                                                   <  2:1014534/100‑1 (MQ=255)
                                   gCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCGTATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTcccc                                                                  >  1:1811809/1‑100 (MQ=255)
                                    cccGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCt                                                                 <  2:3157603/100‑1 (MQ=255)
                                     ccGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCtt                                                                >  2:899337/1‑100 (MQ=255)
                                      cGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTg                                                               <  2:2390916/100‑1 (MQ=255)
                                       gACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGc                                                              >  2:2760760/1‑100 (MQ=255)
                                       gACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGc                                                              >  2:2669104/1‑100 (MQ=255)
                                         cccGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGt                                                            >  2:1099332/1‑100 (MQ=255)
                                           cGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTAc                                                          >  1:1947164/1‑100 (MQ=255)
                                           cGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTAc                                                          >  2:2084493/1‑100 (MQ=255)
                                           cGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTAc                                                          >  1:1587695/1‑100 (MQ=255)
                                           cGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTAc                                                          >  2:506918/1‑100 (MQ=255)
                                           cGTACAAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTAc                                                          <  1:1649762/100‑1 (MQ=255)
                                            gTACGAAAGGGTTCAACCTAAACATTCCATATATCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCAAATATTCTTCACGTCCCCTTGCGt                                                            >  1:1834859/1‑97 (MQ=255)
                                               cGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgca                                                      >  2:1896466/1‑100 (MQ=255)
                                               cGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgca                                                      >  2:1794158/1‑100 (MQ=255)
                                               cGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgca                                                      <  1:2133542/100‑1 (MQ=255)
                                                 aaaGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgcagc                                                    >  2:2844143/1‑100 (MQ=255)
                                                 aaaGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgcagc                                                    >  1:954016/1‑100 (MQ=255)
                                                  aaGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgcagca                                                   <  1:1760597/100‑1 (MQ=255)
                                                      gTTCAACCTAAACATTCCATATCTCCTTGATGTAACTACGAAAGTGAAATGGGTGATCGCTCACATCGTCTTCACGTTCCCTTGCGTACGGCAGCATcgc                                               <  1:1964013/100‑1 (MQ=255)
                                                      gTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATcgc                                               <  1:2682728/100‑1 (MQ=255)
                                                         cAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGt                                            >  1:1646062/1‑100 (MQ=255)
                                                         cAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGt                                            >  1:1346953/1‑100 (MQ=255)
                                                          aaCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTg                                           <  1:3053932/100‑1 (MQ=255)
                                                          aaCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTg                                           <  1:716374/100‑1 (MQ=255)
                                                            ccTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGa                                         >  2:3100915/1‑100 (MQ=255)
                                                            ccTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGa                                         >  1:1353936/1‑100 (MQ=255)
                                                            ccTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGa                                         >  1:1080234/1‑100 (MQ=255)
                                                             cTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGa                                         <  1:2934879/99‑1 (MQ=255)
                                                             cTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGAc                                        <  1:2356184/100‑1 (MQ=255)
                                                             cTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGAc                                        <  1:1641587/100‑1 (MQ=255)
                                                             cTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTGCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGAc                                        <  1:1263249/100‑1 (MQ=255)
                                                              tAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACt                                       >  1:177916/1‑100 (MQ=255)
                                                                aaCATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTa                                     >  1:65617/1‑100 (MQ=255)
                                                                  cATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAtt                                   <  1:1979481/100‑1 (MQ=255)
                                                                  cATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAtt                                   >  1:1744081/1‑100 (MQ=255)
                                                                  cATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAtt                                   >  1:1138497/1‑100 (MQ=255)
                                                                    ttCCATATCTCCATGAGGTAACTACGAAAATGAAATTGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAtttt                                 <  1:981359/100‑1 (MQ=255)
                                                                     tCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAttttt                                >  2:499772/1‑100 (MQ=255)
                                                                      ccATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAtttttt                               >  2:2113885/1‑100 (MQ=255)
                                                                       cATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTa                              <  1:1329572/100‑1 (MQ=255)
                                                                       cATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAtttttt                               <  1:2449594/99‑1 (MQ=255)
                                                                        atatCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTAt                             <  2:2921628/100‑1 (MQ=255)
                                                                         tatCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATg                            <  1:2882060/100‑1 (MQ=255)
                                                                         tatCTCCATGAGGTAAATACGAAAATGAAATGGGTTATCGCTCACATATTCTTCACGTCCCCTTGCGTATGGCAGCATCGCGGAGGACTTACTTTTTAt                             >  1:1977488/1‑99 (MQ=255)
                                                                          atCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTAt                             >  1:1474164/1‑98 (MQ=255)
                                                                            ctcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAg                         >  2:813697/1‑100 (MQ=255)
                                                                             tcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTGACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAgg                        >  1:2085714/1‑100 (MQ=255)
                                                                             tcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAgg                        >  1:2155722/1‑100 (MQ=255)
                                                                             tcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAgg                        >  1:3040095/1‑100 (MQ=255)
                                                                             tcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAgg                        >  1:3104813/1‑100 (MQ=255)
                                                                              ccATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGt                       <  2:2917268/100‑1 (MQ=255)
                                                                               caggAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTg                      <  1:11657/97‑1 (MQ=255)
                                                                                 tGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGat                    <  1:1099329/100‑1 (MQ=255)
                                                                                  gAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCCTCGCGGTGGACTTATTTTTTATGCAGGTGatc                   <  2:3217438/100‑1 (MQ=255)
                                                                                    ggTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGatcat                 >  2:1407807/1‑100 (MQ=255)
                                                                                     gTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGatcatc                <  2:1384728/100‑1 (MQ=255)
                                                                                      tAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATcc               >  2:886485/1‑100 (MQ=255)
                                                                                        aCTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCa             <  2:3348929/100‑1 (MQ=255)
                                                                                        aCTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCa             <  2:411596/100‑1 (MQ=255)
                                                                                        aCTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCa             <  2:709480/100‑1 (MQ=255)
                                                                                         cTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAg            <  2:708742/100‑1 (MQ=255)
                                                                                         cTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAg            <  2:1554676/100‑1 (MQ=255)
                                                                                           aCGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGa           >  2:3376372/1‑99 (MQ=255)
                                                                                            cGAAAATAAAATGGGTTATCGCTCACATCTGCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaaa         >  1:2057127/1‑100 (MQ=255)
                                                                                              aaaatgaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGc       >  1:2464259/1‑100 (MQ=255)
                                                                                               aaatgaaatgGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGcc      >  1:2668759/1‑100 (MQ=255)
                                                                                                  taaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGa   <  2:698255/100‑1 (MQ=255)
                                                                                                   aaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGaa  >  2:3203121/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAA  >  NC_000913/2246701‑2246899

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: