Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 857,578 C→G 16.7% intergenic (+23/+218) ybiT → / ← ybiU ABC‑F family putative regulatory ATPase/DUF1479 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913857,5780CG16.7% 417.3 / 69.4 150intergenic (+23/+218)ybiT/ybiUABC‑F family putative regulatory ATPase/DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base C (56/69);  new base G (20/5);  total (76/74)
Fisher's exact test for biased strand distribution p-value = 1.76e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

Only 100 of 152 total aligned reads displayed.

TGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACT  >  NC_000913/857479‑857677
                                                                                                   |                                                                                                   
tGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGAGTAGTAAAGGGCTCGAGTAAGTGATTATCCGTCACCCCGATAc                                                                                                     >  2:2523085/1‑100 (MQ=255)
 ggAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACt                                                                                                    >  1:3482832/1‑100 (MQ=255)
     aTCTCCCCTGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTc                                                                                                 <  1:1267569/99‑1 (MQ=255)
     aTCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCg                                                                                                <  2:3372378/100‑1 (MQ=255)
       cACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCggg                                                                                              <  2:2144690/100‑1 (MQ=255)
       cACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCggg                                                                                              <  1:3507420/100‑1 (MQ=255)
        aCCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgggg                                                                                             <  1:1241209/100‑1 (MQ=255)
         ccccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGt                                                                                            <  2:2808827/100‑1 (MQ=255)
            cGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAc                                                                                         >  1:2431185/1‑100 (MQ=255)
             ggAACGCGTGCTCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTCATTATCCGTCACCCCGATACTGTCGGGGTGAcg                                                                                        <  2:2601061/100‑1 (MQ=255)
             ggAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcg                                                                                        <  2:756393/100‑1 (MQ=255)
              gAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAcgg                                                                                       <  1:668059/100‑1 (MQ=255)
                  gcgTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACggcggc                                                                                   >  1:341463/1‑100 (MQ=255)
                       aTCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgatg                                                                              >  1:896188/1‑100 (MQ=255)
                       aTCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgatg                                                                              <  2:2002616/100‑1 (MQ=255)
                       aTCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgatg                                                                              >  2:485806/1‑100 (MQ=255)
                        tCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgatga                                                                             <  2:1605341/100‑1 (MQ=255)
                        tCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgatga                                                                             >  2:3251104/1‑100 (MQ=255)
                         cGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGAc                                                                            <  1:1708806/100‑1 (MQ=255)
                           aCTTTAGCGGTAATTACGAAGATTTCCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACaa                                                                          <  2:1235860/100‑1 (MQ=255)
                           aCTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACaa                                                                          <  2:530629/100‑1 (MQ=255)
                             tttAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAAt                                                                        <  1:3611523/100‑1 (MQ=255)
                              ttAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATg                                                                       <  2:1092282/100‑1 (MQ=255)
                               tAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGc                                                                      >  2:3499778/1‑100 (MQ=255)
                                aGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCa                                                                     <  1:1676000/100‑1 (MQ=255)
                                  cGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaaa                                                                   <  1:2220292/100‑1 (MQ=255)
                                  cGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaaa                                                                   >  1:3146755/1‑100 (MQ=255)
                                   ggTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaaaa                                                                  >  2:3400975/1‑100 (MQ=255)
                                         tACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCt                                                            >  1:3461842/1‑100 (MQ=255)
                                         tACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCt                                                            <  1:1759109/100‑1 (MQ=255)
                                          acaaaGATTACCTACGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGCTACTGTCGGGGTGACGGCGGCAGATGACAAATGCAAAAGAGCCt                                                            <  2:3009812/96‑1 (MQ=255)
                                          aCGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTg                                                           <  1:1783192/100‑1 (MQ=255)
                                          aCGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTg                                                           >  2:3447352/1‑100 (MQ=255)
                                           cGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGa                                                          <  1:3482677/100‑1 (MQ=255)
                                           cGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGa                                                          <  1:2429667/100‑1 (MQ=255)
                                            gAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGAt                                                         >  1:1064032/1‑100 (MQ=255)
                                                 ttACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCt                                                    <  1:2840127/100‑1 (MQ=255)
                                                   aCCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTAc                                                  >  2:1347180/1‑100 (MQ=255)
                                                   aCCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTAc                                                  >  2:736987/1‑100 (MQ=255)
                                                   aCCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTAc                                                  <  1:830778/100‑1 (MQ=255)
                                                     cTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGc                                                <  2:578710/100‑1 (MQ=255)
                                                     cTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGc                                                <  1:1230747/100‑1 (MQ=255)
                                                           agtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCa                                          <  2:1302719/100‑1 (MQ=255)
                                                            gtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAg                                         >  1:952581/1‑100 (MQ=255)
                                                             taAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAgg                                        >  2:3272211/1‑100 (MQ=255)
                                                               aaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGcc                                      >  2:1055655/1‑100 (MQ=255)
                                                               aaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGcc                                      >  2:1051071/1‑100 (MQ=255)
                                                                 gggATCGAGTAAGTGATTATCCGTCAACCCGATACTGGCGGGGTGACGGCGGCTGATGAGAAATGCAGAATAGCCTGATGCGCTACGCTTa                                             <  2:3314513/91‑1 (MQ=38)
                                                                   gATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACg                                  >  1:1191489/1‑100 (MQ=255)
                                                                    aTCGAGTAAGTGATTAGCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCTCTACGCTTATCAGGCCTACGc                                 <  2:183427/100‑1 (MQ=255)
                                                                      cGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATCGCCTGATGCGCTACGCTTATCAGGCCTAc                                   >  2:2247496/1‑96 (MQ=255)
                                                                       gAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGt                              >  1:425385/1‑100 (MQ=255)
                                                                       gAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGt                              >  2:3722020/1‑100 (MQ=255)
                                                                       gAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGt                              <  2:3502779/100‑1 (MQ=255)
                                                                        aGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTa                             <  1:3574812/100‑1 (MQ=255)
                                                                          tAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAcc                           <  1:3855896/100‑1 (MQ=255)
                                                                          tAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAcc                           <  1:3475781/100‑1 (MQ=255)
                                                                           aaGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCt                          >  2:1015591/1‑100 (MQ=255)
                                                                            catcagccgccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTg                         >  1:3139644/10‑100 (MQ=25)
                                                                             gTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGc                        >  1:572854/1‑100 (MQ=255)
                                                                              tGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCa                       >  1:3159099/1‑100 (MQ=255)
                                                                               gATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACTCAGTACCTGCaa                      <  1:2600123/100‑1 (MQ=255)
                                                                                aTTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAat                     >  1:307336/1‑100 (MQ=255)
                                                                                aTTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAat                     <  2:409510/100‑1 (MQ=255)
                                                                                aTTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAat                     <  1:3926870/100‑1 (MQ=255)
                                                                                aTTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAAAAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAat                     <  2:22216/100‑1 (MQ=255)
                                                                                   aTCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatata                  <  2:2677083/100‑1 (MQ=255)
                                                                                    gccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat                 >  1:1952118/2‑100 (MQ=25)
                                                                                    gccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat                 >  2:2372591/2‑100 (MQ=25)
                                                                                    gccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat                 >  2:2374044/2‑100 (MQ=25)
                                                                                    gccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat                 >  2:247998/2‑100 (MQ=25)
                                                                                    gccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat                 >  1:3404587/2‑100 (MQ=25)
                                                                                    gccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat                 >  1:130763/2‑100 (MQ=25)
                                                                                     ccGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt                >  2:81868/1‑100 (MQ=255)
                                                                                     ccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATAACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt                >  2:2103186/1‑100 (MQ=25)
                                                                                      cGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTg               >  1:1608506/1‑100 (MQ=255)
                                                                                      cGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTg               <  2:89568/100‑1 (MQ=255)
                                                                                      cGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACACATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTg               >  2:1670352/1‑100 (MQ=25)
                                                                                      cGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTg               >  2:3901065/1‑100 (MQ=25)
                                                                                       gTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTg               <  1:191554/99‑1 (MQ=255)
                                                                                       gTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGAGCTACGCATATCAGGCCTACGCAGTACCTGCAATATAtt                <  1:1827685/98‑1 (MQ=255)
                                                                                       gTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGa              >  1:3123009/1‑100 (MQ=255)
                                                                                       gTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGa              >  1:2287284/1‑100 (MQ=25)
                                                                                        tCACCCCGATACTGTCGGGGTGACGGCTGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGaa             <  1:3496364/100‑1 (MQ=255)
                                                                                        tCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGaa             >  2:2942680/1‑100 (MQ=25)
                                                                                         cACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAAt            >  1:533787/1‑100 (MQ=25)
                                                                                         cACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAAt            >  2:1076705/1‑100 (MQ=21)
                                                                                         cACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAAt            >  1:458532/1‑100 (MQ=25)
                                                                                           ccccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCt          >  1:2762193/1‑100 (MQ=255)
                                                                                           ccccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCt          >  1:2652041/1‑100 (MQ=255)
                                                                                           ccccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCt          >  2:1990060/1‑100 (MQ=255)
                                                                                           ccccGATACTGCCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCt          >  2:3057920/1‑100 (MQ=255)
                                                                                             ccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGc        <  1:1841966/100‑1 (MQ=255)
                                                                                             ccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGc        <  2:1250610/100‑1 (MQ=255)
                                                                                             ccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGc        >  2:3647072/1‑100 (MQ=255)
                                                                                             ccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGc        <  2:3605870/100‑1 (MQ=255)
                                                                                               gATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAg      <  2:2707421/100‑1 (MQ=255)
                                                                                                aTACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAgg     <  2:1723125/100‑1 (MQ=255)
                                                                                                aTACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAgg     <  1:1428355/100‑1 (MQ=255)
                                                                                                   cTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACt  >  2:876925/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACT  >  NC_000913/857479‑857677

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: