Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,246,800 A→G 20.1% 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,8000AG20.1% 303.0 / 68.8 149intergenic (‑31/+263)cirA/lysPcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter
Reads supporting (aligned to +/- strand):  ref base A (70/49);  new base G (20/10);  total (90/59)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

CACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAA  >  NC_000913/2246701‑2246899
                                                                                                   |                                                                                                   
cacaGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaata                                                                                                     >  2:1968341/1‑100 (MQ=255)
cacaGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaata                                                                                                     >  2:2325267/1‑100 (MQ=255)
  caGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaataaa                                                                                                   >  2:3140393/1‑100 (MQ=255)
   aGGCCATGCACAAGAACTAGCGAACAAACACAGCCCGACCCCTACGAAAGGGATCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaataaa                                                                                                   <  2:1373500/99‑1 (MQ=255)
    ggCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaatgaaat                                                                                                 <  1:3723723/100‑1 (MQ=255)
    ggCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaatgaaat                                                                                                 <  1:98454/100‑1 (MQ=255)
    ggCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaataaaat                                                                                                 <  1:3006127/100‑1 (MQ=255)
       cATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATggg                                                                                              <  2:2826006/100‑1 (MQ=255)
       cATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATggg                                                                                              >  2:1709052/1‑100 (MQ=255)
         tGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGtt                                                                                            <  1:2693864/100‑1 (MQ=255)
           cacaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTAt                                                                                          >  2:3624989/1‑100 (MQ=255)
           cacaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTAt                                                                                          >  2:2991575/1‑100 (MQ=255)
            atctgAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAGTGAAATGGGTTATc                                                                                         <  2:3871796/96‑1 (MQ=255)
            acaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATc                                                                                         >  1:86857/1‑100 (MQ=255)
            acaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATgggcccca                                                                                         >  1:311338/1‑95 (MQ=255)
             caAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCg                                                                                        >  2:2776058/1‑100 (MQ=255)
              aaGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGc                                                                                       >  1:3803335/1‑100 (MQ=255)
               aGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCt                                                                                      >  2:3114166/1‑100 (MQ=255)
                gAAATAGCGGACAAACACAGTCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTc                                                                                     <  2:1770812/100‑1 (MQ=255)
                gAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATcagat                                                                                     >  1:3871821/1‑96 (MQ=255)
                gAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTc                                                                                     <  1:1709064/100‑1 (MQ=255)
                  aaTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTcac                                                                                   >  1:115094/1‑100 (MQ=255)
                   aTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTCACTACGAAAATAAAATGGGTTATCGCTcac                                                                                   >  1:766998/1‑99 (MQ=255)
                   aTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTcaca                                                                                  <  2:2084355/100‑1 (MQ=255)
                      gCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACAtct                                                                               <  2:2146213/100‑1 (MQ=255)
                      gCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACAtct                                                                               <  1:3412842/100‑1 (MQ=255)
                       cGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtctt                                                                              >  2:2386043/1‑100 (MQ=255)
                        ggACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACAtct                                                                               >  2:355327/1‑98 (MQ=255)
                         gACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtcttct                                                                            <  1:3625012/100‑1 (MQ=255)
                         gACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACAtcttct                                                                            <  1:720275/100‑1 (MQ=255)
                         gACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACAtcttct                                                                            <  1:1856400/100‑1 (MQ=255)
                           cAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCttcttc                                                                          <  1:3962897/100‑1 (MQ=255)
                           cAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAAGAAAATGGGTTATCGCTCACATCttcttc                                                                          >  1:3704537/1‑100 (MQ=255)
                            aaaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCa                                                                         <  1:149041/100‑1 (MQ=255)
                            aaaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCa                                                                         <  1:2757553/100‑1 (MQ=255)
                             aaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCAc                                                                        >  1:314802/1‑100 (MQ=255)
                             aaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAAAATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTcccacactattcac                                                                        >  1:2406718/1‑90 (MQ=255)
                              acacAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACg                                                                       <  1:582862/100‑1 (MQ=255)
                               cacaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGt                                                                      >  1:2691098/1‑100 (MQ=255)
                               cacaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGt                                                                      >  1:458088/1‑100 (MQ=255)
                               cacaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGt                                                                      >  1:2542646/1‑100 (MQ=255)
                                acaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTc                                                                     <  2:971416/100‑1 (MQ=255)
                                 caGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTcc                                                                    >  1:1229204/1‑100 (MQ=255)
                                   gCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTcccc                                                                  >  1:2883759/1‑100 (MQ=255)
                                    cccGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCt                                                                 <  1:2235774/100‑1 (MQ=255)
                                    cccGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCt                                                                 <  1:3057130/100‑1 (MQ=255)
                                     ccGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTccc                                                                   >  2:3823049/1‑97 (MQ=255)
                                     ccGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCtt                                                                <  1:940555/100‑1 (MQ=255)
                                     ccGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGGAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCtt                                                                >  2:268116/1‑100 (MQ=255)
                                          ccGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTa                                                           >  2:2436063/1‑100 (MQ=255)
                                           cGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTAc                                                          >  2:2758457/1‑100 (MQ=255)
                                            gTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTAAGAAAATAAAATGGATTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACg                                                         >  1:1007872/1‑100 (MQ=255)
                                             tACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACgg                                                        <  1:3103474/100‑1 (MQ=255)
                                              aCGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgc                                                       >  1:1038343/1‑100 (MQ=255)
                                              aCGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgc                                                       <  2:1997011/100‑1 (MQ=255)
                                              aCGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgc                                                       >  2:3316746/1‑100 (MQ=255)
                                              aCGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgc                                                       <  1:1464566/100‑1 (MQ=255)
                                               cGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgca                                                      >  2:224714/1‑100 (MQ=255)
                                                gAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgcag                                                     >  1:407403/1‑100 (MQ=255)
                                                gAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgcag                                                     >  1:1007474/1‑100 (MQ=255)
                                                    gggTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCCGCATc                                                 >  2:4053896/1‑100 (MQ=255)
                                                     ggTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATcg                                                >  1:129925/1‑100 (MQ=255)
                                                      gTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATcgc                                               <  2:694311/100‑1 (MQ=255)
                                                       ttCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATcgcg                                              <  2:1202466/100‑1 (MQ=255)
                                                        tCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCgg                                             <  2:3221121/100‑1 (MQ=255)
                                                        tCAAACTAAACATTCCATATCTCCATGAGGAAACAAAAAAAAAAAAATGGGTTATCGCTCACATCTTCTTCACGTCCTCTTGCGTACGGCAGCATcgc                                               <  2:2519035/98‑1 (MQ=255)
                                                          aaCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTg                                           >  2:526989/1‑100 (MQ=255)
                                                           aCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTgg                                          <  2:2424894/100‑1 (MQ=255)
                                                           aCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTgg                                          >  2:2383346/1‑100 (MQ=255)
                                                           aCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTgg                                          <  1:3292387/100‑1 (MQ=255)
                                                             cTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGAc                                        <  1:2872488/100‑1 (MQ=255)
                                                             cTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGAc                                        <  2:3423140/100‑1 (MQ=255)
                                                                 aCATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAt                                    >  1:2286780/1‑100 (MQ=255)
                                                                   aTTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAttt                                  >  1:418164/1‑100 (MQ=255)
                                                                     tCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAttttt                                <  1:1063901/100‑1 (MQ=255)
                                                                       cATATCTCCATGAGGTAACTACGAAAATGAAATGGATTATCGCTCACATCTTCTTCACGTCCTCTTGCGTACGGCAGCATCGCGGTGGACTTAtttttt                               >  2:4052905/1‑99 (MQ=255)
                                                                       cATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAttttttn                              <  1:3714483/100‑2 (MQ=255)
                                                                       cATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTa                              <  2:3834668/100‑1 (MQ=255)
                                                                        atatCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTAt                             <  1:2753317/100‑1 (MQ=255)
                                                                         tatCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATg                            >  1:1984603/1‑100 (MQ=255)
                                                                          atCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGc                           >  1:738302/1‑100 (MQ=255)
                                                                           tctcCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATTCa                          >  2:859998/1‑100 (MQ=255)
                                                                           tctcCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCa                          >  2:1450418/1‑100 (MQ=255)
                                                                             tcCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAgg                        <  1:1295237/100‑1 (MQ=255)
                                                                               ccgatcttAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTg                      <  1:972543/93‑1 (MQ=255)
                                                                                aTGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGa                     >  1:2425305/1‑100 (MQ=255)
                                                                                aTGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGa                     <  1:3218883/100‑1 (MQ=255)
                                                                                  gAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGatc                   >  2:1229006/1‑100 (MQ=255)
                                                                                  gAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGatc                   >  2:1357304/1‑100 (MQ=255)
                                                                                    ggTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGatcat                 >  1:3510850/1‑100 (MQ=255)
                                                                                      tAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATcc               <  1:2436077/100‑1 (MQ=255)
                                                                                      tAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATCCAGGTGatc                   >  2:3052834/1‑96 (MQ=255)
                                                                                       aaCTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATccc              <  1:3211628/100‑1 (MQ=255)
                                                                                       aaCTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATccc              <  2:1542530/100‑1 (MQ=255)
                                                                                         cTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAg            >  1:827289/1‑100 (MQ=255)
                                                                                         cTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAg            <  2:3659697/100‑1 (MQ=255)
                                                                                         cTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAg            >  2:4011533/1‑100 (MQ=255)
                                                                                         cTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGACCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTATTTATGCAGGTGATCATcc               >  2:2207724/1‑97 (MQ=255)
                                                                                               aaataaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGcc      <  1:3952431/100‑1 (MQ=255)
                                                                                                   aaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGaa  >  1:2269321/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 < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: