Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,246,800 A→G 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,8000AG96.2% 511.0 / 4.1 157intergenic (‑31/+263)cirA/lysPcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter
Reads supporting (aligned to +/- strand):  ref base A (5/1);  new base G (82/69);  total (87/70)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 167 total aligned reads displayed.

CAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCG  >  NC_000913/2246703‑2246896
                                                                                                 |                                                                                                
caGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaatgaa                                                                                                >  1:876014/1‑100 (MQ=255)
  ggCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaatgaaat                                                                                              >  2:3182431/1‑100 (MQ=255)
   gCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAaaatgaaatg                                                                                             >  2:1281938/1‑100 (MQ=255)
    ccATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATgg                                                                                            <  2:1972370/100‑1 (MQ=255)
      aggCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGt                                                                                          <  1:1679188/98‑1 (MQ=255)
         cacaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTAt                                                                                       >  1:3067976/1‑100 (MQ=255)
         cacaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTAt                                                                                       >  1:3027087/1‑100 (MQ=255)
          acaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCCTGAGGTAACTACGAAAAAGAAATGGGTTAt                                                                                       >  1:2184828/1‑99 (MQ=255)
          acaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATc                                                                                      >  2:1717448/1‑100 (MQ=255)
          acaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATc                                                                                      <  2:2804749/100‑1 (MQ=255)
           caAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCg                                                                                     >  2:2456811/1‑100 (MQ=255)
           caAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCg                                                                                     >  1:510738/1‑100 (MQ=255)
           caAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCg                                                                                     >  1:483807/1‑100 (MQ=255)
           caAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCg                                                                                     >  1:2152260/1‑100 (MQ=255)
            aaGAAATAGCTGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGc                                                                                    <  1:2730624/100‑1 (MQ=255)
             aGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCt                                                                                   >  1:2963026/1‑100 (MQ=255)
              gAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTc                                                                                  >  2:1582829/1‑100 (MQ=255)
              gAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTc                                                                                  >  2:1822837/1‑100 (MQ=255)
                 aTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTcaca                                                                               >  1:1785391/1‑100 (MQ=255)
                 aTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTcaca                                                                               >  2:1215852/1‑100 (MQ=255)
                   aGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtc                                                                             >  2:1148640/1‑100 (MQ=255)
                    gCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtct                                                                            >  1:1964095/1‑100 (MQ=255)
                    gCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACAtct                                                                            >  1:746086/1‑100 (MQ=255)
                      ggACAAACACAGCCCGACCCGTACGAAAGGGTTGAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtattc                                                                          >  1:592680/1‑100 (MQ=255)
                       gACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtcttct                                                                         <  2:2933712/100‑1 (MQ=255)
                       gACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtcttct                                                                         <  2:3068026/100‑1 (MQ=255)
                        aCAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATcttctt                                                                        <  2:48228/100‑1 (MQ=255)
                         cAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAGATGGGTTATCGCTCACATCttcttc                                                                       >  2:392817/1‑100 (MQ=255)
                          aaaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCa                                                                      <  2:1377860/100‑1 (MQ=255)
                           aaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCAc                                                                     >  1:2202194/1‑100 (MQ=255)
                           aaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCAc                                                                     <  1:709824/100‑1 (MQ=255)
                           aaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCAc                                                                     <  1:1613931/100‑1 (MQ=255)
                           aaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCAc                                                                     >  1:396185/1‑100 (MQ=255)
                             cacaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCAGAtcggaagagcac                                                                   >  2:758991/1‑90 (MQ=255)
                             cacaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCAc                                                                     >  1:3413045/1‑98 (MQ=255)
                             cacaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGt                                                                   >  2:2059749/1‑100 (MQ=255)
                              acaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTc                                                                  <  1:2035365/100‑1 (MQ=255)
                              acaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTc                                                                  <  1:1281921/100‑1 (MQ=255)
                               caGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTcc                                                                 <  1:1484495/100‑1 (MQ=255)
                                 gCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTcccc                                                               <  1:2688015/100‑1 (MQ=255)
                                   ccGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCtt                                                             >  2:2408114/1‑100 (MQ=255)
                                    cGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCTCTCACATCTTCTTCACGTCCCCtt                                                             >  2:146207/1‑99 (MQ=255)
                                     gACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTg                                                            <  2:984311/99‑1 (MQ=255)
                                       cccGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGt                                                         <  2:533000/100‑1 (MQ=255)
                                           tACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACgg                                                     >  1:858821/1‑100 (MQ=255)
                                            aCGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgc                                                    <  2:2724863/100‑1 (MQ=255)
                                            aCGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGCACGgc                                                    <  2:1008550/100‑1 (MQ=255)
                                             cGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgca                                                   <  2:137962/100‑1 (MQ=255)
                                             cGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgca                                                   >  1:297223/1‑100 (MQ=255)
                                              gAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgcag                                                  <  1:2014759/100‑1 (MQ=255)
                                               aaaGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgcagc                                                 >  2:2603244/1‑100 (MQ=255)
                                               aaaGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgcag                                                  >  1:2865697/1‑99 (MQ=255)
                                                    gTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATcgc                                            >  2:2015646/1‑100 (MQ=255)
                                                       cAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGt                                         <  2:3027136/100‑1 (MQ=255)
                                                       cAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGt                                         <  1:1101670/100‑1 (MQ=255)
                                                         aCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTgg                                       >  2:2213832/1‑100 (MQ=255)
                                                         aCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTgg                                       >  1:325773/1‑100 (MQ=255)
                                                         aCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTgg                                       <  2:2019624/100‑1 (MQ=255)
                                                          ccTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGa                                      >  1:3091466/1‑100 (MQ=255)
                                                          ccTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGa                                      >  2:2180259/1‑100 (MQ=255)
                                                          ccTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGa                                      <  2:831442/100‑1 (MQ=255)
                                                          ccTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGa                                      >  1:3200750/1‑100 (MQ=255)
                                                          ccTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGa                                      >  2:3121166/1‑100 (MQ=255)
                                                           cTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGAc                                     >  1:1110084/1‑100 (MQ=255)
                                                            tAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACt                                    <  1:1116465/100‑1 (MQ=255)
                                                            tAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACt                                    <  2:1490667/100‑1 (MQ=255)
                                                            tAAACATTCCATACCTCCATGAGGTAACTACGAAAATGAAATGGGTTATAGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACt                                    >  2:2333617/1‑100 (MQ=255)
                                                             aaaCATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACtt                                   >  2:107671/1‑100 (MQ=255)
                                                               aCATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAt                                 <  2:3373369/100‑1 (MQ=255)
                                                               aCATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAt                                 >  1:1020453/1‑100 (MQ=255)
                                                                 aTTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAttt                               >  2:829054/1‑100 (MQ=255)
                                                                   tCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAttttt                             <  2:3181787/100‑1 (MQ=255)
                                                                    ccATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAtttttt                            >  2:1661694/1‑100 (MQ=255)
                                                                    ccATATCTCCATGAGGTAACGACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAtttttt                            >  1:2126291/1‑100 (MQ=255)
                                                                     cATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTa                           <  2:855490/100‑1 (MQ=255)
                                                                     cATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTa                           <  1:1783357/100‑1 (MQ=255)
                                                                     cATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTa                           <  2:3463455/100‑1 (MQ=255)
                                                                      atatCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTAt                          <  1:244768/100‑1 (MQ=255)
                                                                      atatCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTAt                          <  1:2313180/100‑1 (MQ=255)
                                                                       tatCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATg                         <  1:3019436/100‑1 (MQ=255)
                                                                       tatCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTACCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATg                         >  1:2382478/1‑100 (MQ=255)
                                                                         tctcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCa                       >  1:3123792/1‑100 (MQ=255)
                                                                          ctcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAg                      >  1:315231/1‑100 (MQ=255)
                                                                          ctcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAg                      >  1:1267658/1‑100 (MQ=255)
                                                                               tGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGat                 >  1:2922159/1‑100 (MQ=255)
                                                                                gAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATGTTCTTCACGTCCCCTTGCGTACGGCGGCATCGCGGTGGACTTATTTTTTCTGCAGGTg                   >  1:1187270/1‑97 (MQ=255)
                                                                                 acgTAACTACGAAAATGAAATTGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTAGACTTATTTTTTATGCAGGTGatca               <  1:3335170/98‑1 (MQ=255)
                                                                                    tAACTACGAACATAAAATGGGTTATCGCACCCACCCTCCTCACGTCAACTAGAGTACTGCCGCATCCCGGTAGACTTATCTTTTATGCAGGTGACCATcc            >  1:2151174/1‑100 (MQ=255)
                                                                                     aaCTACGAAAGTGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGGGGAGGACTTATTTTTTATGCAGGTGATCATccc           <  1:2279869/100‑1 (MQ=255)
                                                                                     aaCTACGAAAATGAAATGGGTTAGCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATccc           <  1:708571/100‑1 (MQ=255)
                                                                                      aCTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCa          <  1:1550840/100‑1 (MQ=255)
                                                                                       cTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAg         <  2:3200802/100‑1 (MQ=255)
                                                                                         aCGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaa       >  2:1631280/1‑100 (MQ=255)
                                                                                          cGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaaa      >  2:1207255/1‑100 (MQ=255)
                                                                                          cGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaaa      >  2:1794050/1‑100 (MQ=255)
                                                                                          cGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGCTGGACTTATTTTTTATGCAGGTGATCATCCCAGaaa      >  1:1045341/1‑100 (MQ=255)
                                                                                           gAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAg     >  2:183573/1‑100 (MQ=255)
                                                                                            aaaatgaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGc    >  1:852781/1‑100 (MQ=255)
                                                                                             aaatgaaatgGGTTATCGCTCACATCTTCTTCACGTCCCCTAGCGTACGGCAGCATCGCGGCGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGcc   >  1:1895212/1‑100 (MQ=255)
                                                                                              aatgaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCg  >  1:3410027/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCG  >  NC_000913/2246703‑2246896

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: