Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,8480AG53.9% 113.4 / inf 508F45F (TTT→TTCompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base A (131/101);  new base G (153/118);  total (285/223)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 518 total aligned reads displayed.

GTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGCTACTAACAGAGCA  >  NC_000913/986750‑986947
                                                                                                  |                                                                                                   
  cAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGagaaa                                                                                                  <  2:748873/100‑1 (MQ=255)
    gATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAATTGTTTGCACCGTTACCCTTGGAAAaata                                                                                                >  1:4697504/1‑100 (MQ=255)
    gATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAaata                                                                                                >  1:1324262/1‑100 (MQ=255)
        ggAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGc                                                                                            <  2:2718905/100‑1 (MQ=255)
          aaTTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAg                                                                                          <  2:3938071/100‑1 (MQ=255)
          aaTTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAg                                                                                          <  2:674779/100‑1 (MQ=255)
           aTTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGa                                                                                         >  1:180804/1‑100 (MQ=255)
                 ttGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACa                                                                                   >  1:2431244/1‑100 (MQ=255)
                   gAGTTTCCCCTTTAACGCCCAGAGGGGCATCGGTACTGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGACCAATACTGAAGACCACCAGc                                                                                 <  2:231447/100‑1 (MQ=255)
                       ttCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTa                                                                             >  1:2886941/1‑100 (MQ=255)
                       ttCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACCGCTTTa                                                                             >  2:4872504/1‑100 (MQ=255)
                       ttCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTa                                                                             <  2:1421850/100‑1 (MQ=255)
                       ttCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTa                                                                             >  1:886899/1‑100 (MQ=255)
                       ctCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTa                                                                             >  1:3580878/2‑100 (MQ=255)
                       ctCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTa                                                                             >  1:2542308/2‑100 (MQ=255)
                       atCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTa                                                                             >  1:3126190/2‑100 (MQ=255)
                        tCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAc                                                                            >  1:5030558/1‑100 (MQ=255)
                        tCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTCTAc                                                                            <  2:3128451/100‑1 (MQ=255)
                         ccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAcc                                                                           >  1:3724318/1‑100 (MQ=255)
                         ccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAcc                                                                           >  1:4222570/1‑100 (MQ=255)
                         ccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAcc                                                                           >  1:4323423/1‑100 (MQ=255)
                         ccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAcc                                                                           >  1:5138264/1‑100 (MQ=255)
                         ccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAcc                                                                           >  1:3388689/1‑100 (MQ=255)
                         ccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAcc                                                                           >  1:1115662/1‑100 (MQ=255)
                         ccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAcc                                                                           >  1:1254389/1‑100 (MQ=255)
                         ccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAcc                                                                           >  1:1460016/1‑100 (MQ=255)
                         ccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAcc                                                                           >  1:1702435/1‑100 (MQ=255)
                         ccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAcc                                                                           >  1:2346886/1‑100 (MQ=255)
                         ccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAcc                                                                           >  1:1756379/1‑100 (MQ=255)
                         ccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAcc                                                                           >  1:1889718/1‑100 (MQ=255)
                         ccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCGTTAcc                                                                           >  1:1489281/1‑100 (MQ=255)
                          cccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCg                                                                          >  1:2657692/1‑100 (MQ=255)
                          cccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCg                                                                          <  2:2808298/100‑1 (MQ=255)
                          cccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCg                                                                          >  2:3174522/1‑100 (MQ=255)
                          cccTTTAAAACCAAGACGGGCATAGGTAATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAGCAGCTTTACCg                                                                          >  1:4538440/1‑100 (MQ=255)
                            cTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCCTCACCGTa                                                                        >  1:3060772/1‑100 (MQ=255)
                             tttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTAc                                                                       <  2:2477445/100‑1 (MQ=255)
                             tttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTAc                                                                       <  2:255762/100‑1 (MQ=255)
                             tttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTAc                                                                       >  1:4445516/1‑100 (MQ=255)
                                aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGa                                                                    >  1:2152489/1‑100 (MQ=255)
                                 aaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGAt                                                                   >  1:1933196/1‑100 (MQ=255)
                                 aaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGAt                                                                   >  1:1201354/1‑100 (MQ=255)
                                    ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTa                                                                >  1:3816282/1‑100 (MQ=255)
                                     cAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTAc                                                               <  1:4267573/100‑1 (MQ=255)
                                      aaGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACt                                                              <  2:1178794/100‑1 (MQ=255)
                                       aGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACtt                                                             >  2:3456337/1‑100 (MQ=255)
                                          cGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTttgt                                                          >  1:4887497/1‑100 (MQ=255)
                                           gatcttAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTttgtt                                                         <  1:344532/95‑1 (MQ=255)
                                            ggCATAGGTCATGTCGCCTTTGCCACCGTAACTGTTTTCAACGTTACCCTTGGAGAACTAATGCAGTCCAACAGCTTTACCGTACAGATCTACTttcttn                                                        <  2:1042349/100‑2 (MQ=255)
                                            ggCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTttgttg                                                        <  2:2272432/100‑1 (MQ=255)
                                            ggCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTttgttg                                                        <  2:2351237/100‑1 (MQ=255)
                                            ggCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTttgttg                                                        <  2:1613520/100‑1 (MQ=255)
                                            ggCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTttgttg                                                        <  1:2588057/100‑1 (MQ=255)
                                            ggCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTttgttg                                                        <  2:4404474/100‑1 (MQ=255)
                                            ggCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTttgttg                                                        >  2:5012232/1‑100 (MQ=255)
                                             gCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGc                                                       >  1:1530006/1‑100 (MQ=255)
                                               aTAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCa                                                     >  1:2798496/1‑100 (MQ=255)
                                                     cATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTg                                               <  2:3067625/100‑1 (MQ=255)
                                                       tGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGtt                                             >  2:1731032/1‑100 (MQ=255)
                                                        gTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTta                                            >  1:783744/1‑100 (MQ=255)
                                                        gTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTta                                            >  1:1631729/1‑100 (MQ=255)
                                                        gTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTta                                            >  1:524006/1‑100 (MQ=255)
                                                         tCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTtat                                           <  2:3513305/100‑1 (MQ=255)
                                                         tCGCCATTGCCACCGTAACTGTTATCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTCTTGCCATCTTTGTtat                                           >  2:4248573/1‑100 (MQ=255)
                                                          cGCCATTGCCACCGTAAATGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCCACAGCTTTACCGTACAGACCTACTTTGTTGCCATCCTTGTtata                                          >  2:465379/1‑100 (MQ=255)
                                                           gCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAg                                         >  1:2208881/1‑100 (MQ=255)
                                                           gCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAg                                         <  2:2137033/100‑1 (MQ=255)
                                                           gCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAg                                         >  2:3739972/1‑100 (MQ=255)
                                                             cATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGAt                                       <  2:1892773/100‑1 (MQ=255)
                                                             cATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGAt                                       >  1:746589/1‑100 (MQ=255)
                                                                tGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTc                                    <  2:3795849/100‑1 (MQ=255)
                                                                tGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTc                                    <  2:4045490/100‑1 (MQ=255)
                                                                tGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTc                                    <  2:2936352/100‑1 (MQ=255)
                                                                tGCCACCGTAACTGGTGTCACTGTTTCCCTTGGAGAGATAATGCATACCGAGAGCTTTACCGTACACATCTACCTTGGTGCCATCTTTGTTATAGATTTc                                    <  2:4799005/100‑1 (MQ=255)
                                                                  ccACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTg                                  <  2:3688459/100‑1 (MQ=255)
                                                                   cACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGc                                 >  1:2958224/1‑100 (MQ=255)
                                                                   cACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGc                                 <  2:1989551/100‑1 (MQ=255)
                                                                    aCCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCa                                >  2:834502/1‑100 (MQ=255)
                                                                    aCCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCa                                >  1:102208/1‑100 (MQ=255)
                                                                           cTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCAACAGATTTACCGTACAGAACTACATTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTg                         <  2:937074/100‑1 (MQ=255)
                                                                            tGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGc                        <  2:1287798/100‑1 (MQ=255)
                                                                             gTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCa                       <  2:531849/100‑1 (MQ=255)
                                                                                 acACCGTTACCCTTGGAGAAATAATGCAGACCGAGAGCTTTACCGTACAGATCAACTTTGTTGCCATCTTTGTTGTAGATTTCTGCAGCGTTTGCAGTAc                   >  1:4789841/2‑100 (MQ=255)
                                                                                  cACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTAcc                  >  1:3772769/1‑100 (MQ=255)
                                                                                  cACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTAcc                  <  2:4952090/100‑1 (MQ=255)
                                                                                    ccGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTg                >  1:2074029/1‑100 (MQ=255)
                                                                                      gTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGct              <  2:4180545/100‑1 (MQ=255)
                                                                                      gTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGct              <  2:146913/100‑1 (MQ=255)
                                                                                      gTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGct              <  1:3210367/100‑1 (MQ=255)
                                                                                          cccTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGctacta          >  1:246398/1‑100 (MQ=255)
                                                                                          cccTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGctacta          >  1:4962272/1‑100 (MQ=255)
                                                                                          cccTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGctacta          >  1:4366440/1‑100 (MQ=255)
                                                                                           ncTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGCTACTaa         >  1:4110127/2‑100 (MQ=255)
                                                                                           ccTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGCTACTaa         >  1:3509003/1‑100 (MQ=255)
                                                                                           ccTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGCTACTaa         >  1:3914904/1‑100 (MQ=255)
                                                                                            cTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGCTACTAAc        >  1:636542/1‑100 (MQ=255)
                                                                                            cTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGCTACTAAc        <  2:447198/100‑1 (MQ=255)
                                                                                              tGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGCTACTAACag      <  2:2447982/100‑1 (MQ=255)
                                                                                              tGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACATTGTTGCCATATTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGCTACTAACag      >  1:574543/1‑100 (MQ=37)
                                                                                                  aaaaTAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGCTACTAACAGAGCa  <  2:3581501/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTATAGATTTCTGCAGCGTTTGCAGTACCTGCTACTAACAGAGCA  >  NC_000913/986750‑986947

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: