Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,579,563 G→T 16.0% intergenic (‑221/+70) yddA ← / ← ydeM putative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,579,5630GT16.0% 416.5 / 59.7 159intergenic (‑221/+70)yddA/ydeMputative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base G (58/73);  new base T (14/11);  total (73/86)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

Only 100 of 165 total aligned reads displayed.

TTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATC  >  NC_000913/1579464‑1579662
                                                                                                   |                                                                                                   
ttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAg                                                                                                     <  1:831029/100‑1 (MQ=255)
 tatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAgg                                                                                                    >  2:4053300/1‑100 (MQ=255)
     cTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAtt                                                                                                <  2:3505869/100‑1 (MQ=255)
     cTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAtt                                                                                                <  1:339794/100‑1 (MQ=255)
      tATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAttt                                                                                               >  1:1019793/1‑100 (MQ=255)
      tATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAttt                                                                                               >  2:3008030/1‑100 (MQ=255)
         tGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATg                                                                                            <  1:1802103/100‑1 (MQ=255)
          gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGGCGAAACACCCGACAGGAATTTATGt                                                                                           >  1:1592649/1‑100 (MQ=255)
           aTAACCAAGGATAGTGCTAATGCGCTGATTCATACCACATAAGCCAGGATTTTAGCAAGGGTAGGGGATGATGACGAAACACCCGACAGGAATTTAt                                                                                             <  1:3167691/97‑1 (MQ=255)
              aTCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCggg                                                                                       <  1:4335713/100‑1 (MQ=255)
               tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGTATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                      <  1:3344257/100‑1 (MQ=255)
               tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                      <  2:947489/100‑1 (MQ=255)
               tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                      <  1:1754840/100‑1 (MQ=255)
               tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                      <  1:1224811/100‑1 (MQ=255)
                cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTg                                                                                     >  2:3501661/1‑100 (MQ=255)
                  aGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGcc                                                                                   <  2:491393/100‑1 (MQ=255)
                     atatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTgtttcg                                                                                <  2:2311406/100‑6 (MQ=37)
                     atatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTACCTTg                                                                                <  1:3587024/100‑1 (MQ=255)
                     atatTGCTAATGCCCTGACGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCTGACAGGAATTTATGTCGGGTGCCTTg                                                                                <  2:2677675/100‑1 (MQ=255)
                      tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGt                                                                               <  2:4029375/100‑1 (MQ=255)
                      tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGt                                                                               >  2:356637/1‑100 (MQ=255)
                           cTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                          <  1:1741560/100‑1 (MQ=255)
                           cTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                          <  1:1523218/100‑1 (MQ=255)
                            tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                         <  1:3293018/100‑1 (MQ=255)
                            tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                         >  2:3251694/1‑100 (MQ=255)
                             aaTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTc                                                                        >  2:2903923/1‑100 (MQ=255)
                              aTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTc                                                                        <  2:1391972/99‑1 (MQ=255)
                                    tGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGCAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaagaa                                                                 >  1:3711361/1‑100 (MQ=255)
                                     gATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaagaag                                                                <  2:270405/100‑1 (MQ=255)
                                      aTGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAagg                                                               <  2:3605552/100‑1 (MQ=255)
                                      aTGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAagg                                                               <  2:4131649/100‑1 (MQ=255)
                                      aTGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAagg                                                               <  2:426098/100‑1 (MQ=255)
                                         cATACCACGTAAGCCAGGATTTTCGCAAGGGAGGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggagg                                                            <  1:1172696/100‑1 (MQ=255)
                                         cATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggagg                                                            <  2:4369770/100‑1 (MQ=255)
                                          aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                           <  2:1846505/100‑1 (MQ=255)
                                           tACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCt                                                          >  1:1932547/1‑100 (MQ=255)
                                           tACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCt                                                          >  1:925923/1‑100 (MQ=255)
                                            aCCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTa                                                         >  2:1700002/1‑100 (MQ=255)
                                                   aaGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATgg                                                  >  2:2756333/1‑100 (MQ=255)
                                                    aGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGa                                                 <  2:3741955/100‑1 (MQ=255)
                                                     gCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCGTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                <  2:1081998/100‑1 (MQ=255)
                                                      ccAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTCAGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGCAGt                                               >  1:4171020/1‑100 (MQ=255)
                                                      ccAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGt                                               >  1:4019508/1‑100 (MQ=255)
                                                      ccAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGt                                               >  2:1535269/1‑100 (MQ=255)
                                                      ccAGGATCTTCGCAAGGAAAGATGATGATGACGAAACACCCGACCGGAATTTATGTCGGGTGCCTTGTCAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                <  1:1930972/99‑1 (MQ=255)
                                                        aGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTa                                             <  1:707821/100‑1 (MQ=255)
                                                        aGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTa                                             >  2:3581585/1‑100 (MQ=255)
                                                        aGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTa                                             <  1:2903907/100‑1 (MQ=255)
                                                          gATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaa                                           <  2:3291529/100‑1 (MQ=255)
                                                             tttCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAaga                                        <  2:2242376/100‑1 (MQ=255)
                                                              ttCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagag                                       >  2:1604933/1‑100 (MQ=255)
                                                               tCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagaga                                      <  2:487224/100‑1 (MQ=255)
                                                               tCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagag                                       <  1:1439910/99‑1 (MQ=255)
                                                                cGCAAGGGAAGATTGTGGTGACGAAACACCCTACAGGAATTTATGTCGGGTGGCTTGTTAAGCTCATAAGAAGTAGGTTAAGAAGTGAGGAAAAagc                                        <  1:2160870/97‑2 (MQ=255)
                                                                cGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagagag                                     >  2:4320062/1‑100 (MQ=255)
                                                                cGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagagag                                     >  1:4258217/1‑100 (MQ=255)
                                                                cGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagagag                                     <  2:2697586/100‑1 (MQ=255)
                                                                 gCAAGGGGAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGc                                    >  1:4124586/1‑100 (MQ=255)
                                                                 gCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagagaac                                    >  1:707987/1‑98 (MQ=255)
                                                                 gCAAGGGAAGATGATGATGACGAAAAACCCGACAGGAATTTATGTCGGGTGCCTAGTTAAGGTCCTAAGAAGTAGGCTAAGAATGGAGTTAAAAGAGAGc                                    >  1:1015915/1‑100 (MQ=255)
                                                                    aGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGtt                                 >  1:3549143/1‑100 (MQ=255)
                                                                     gggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGtta                                <  1:3573417/100‑1 (MQ=255)
                                                                      ggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttat                               >  2:11821/1‑100 (MQ=255)
                                                                       gAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttatt                              <  2:2000596/100‑1 (MQ=255)
                                                                       gAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttatt                              <  2:634120/100‑1 (MQ=255)
                                                                       gAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttatt                              <  2:2940097/100‑1 (MQ=255)
                                                                          gatgatgatGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaat                           <  1:4055092/100‑1 (MQ=255)
                                                                          gatgatgatGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaat                           <  1:845101/100‑1 (MQ=255)
                                                                          gatgatgatGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaan                           <  2:3922355/100‑2 (MQ=255)
                                                                           atgatgatgACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatt                          <  1:4031858/100‑1 (MQ=255)
                                                                             gatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAattat                        >  2:2837990/1‑100 (MQ=255)
                                                                               tgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtc                      <  2:3770741/100‑1 (MQ=255)
                                                                                   gaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCa                  >  2:2754285/1‑100 (MQ=255)
                                                                                      gAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGTAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAttt               <  1:3244475/100‑1 (MQ=255)
                                                                                      gAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAttt               >  2:3356905/1‑100 (MQ=255)
                                                                                      aaggcacCCGACATAAATTCCTGTCGGGTGTCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAttt               >  2:4148844/5‑100 (MQ=37)
                                                                                      aaggcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGACTGGAGTTAAAAGAGAGCGTTATTAATTATTCTACAttt               >  1:351365/5‑100 (MQ=37)
                                                                                       aaaCACCCTACAGGAATTTATGTCGGGTGCCGTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAtt                <  1:1422281/98‑1 (MQ=255)
                                                                                       aaaCACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTg              >  2:1285213/1‑100 (MQ=255)
                                                                                       aaaCACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTg              <  2:3151625/100‑1 (MQ=255)
                                                                                         gcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtt            >  2:112497/2‑100 (MQ=38)
                                                                                         gcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtt            >  1:2647824/2‑100 (MQ=38)
                                                                                         gcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtt            >  1:966751/2‑100 (MQ=38)
                                                                                         acacCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtt            <  2:3119241/100‑1 (MQ=255)
                                                                                          cacCCGACATAAATTCCTGTCGGGTGCCTTGTTACGGCCATAAGAAGGAGGATAAGACTGGAGTTAAAAGCGAGCGTTATTAATTATTCTCCATTTGt             >  1:2409201/1‑98 (MQ=25)
                                                                                          cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt           >  1:952075/1‑100 (MQ=38)
                                                                                          cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt           >  2:3968126/1‑100 (MQ=38)
                                                                                          cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt           >  1:1134056/1‑100 (MQ=38)
                                                                                          cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt           >  1:2392758/1‑100 (MQ=38)
                                                                                          cacCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt           >  2:3215008/1‑100 (MQ=255)
                                                                                           acCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtttt          <  2:509540/100‑1 (MQ=255)
                                                                                           acCCGACACGAATTTATGTCGGGCGCCTTGTTAAGGCCCTAAGAAGGAGGGTAAGAATGGATTTAAAAGAGAGCGTTATTAATTCTTCTCCATTTGtttt          >  1:4118852/1‑100 (MQ=255)
                                                                                             ccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGc        <  2:1013976/100‑1 (MQ=255)
                                                                                             ccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGc        <  1:1627682/100‑1 (MQ=255)
                                                                                             ccGACAGGAATTTATGTCGGGTGACTTGTTACGGTCATAAAAAGGACTCTAAGAATACAGTTAAAAGAGAGCGTCATTAATTATTCTCGATTTGtttt          >  1:4226133/1‑98 (MQ=255)
                                                                                              cGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCa       <  2:1840604/100‑1 (MQ=255)
                                                                                                 cAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACa    <  1:1646810/100‑1 (MQ=255)
                                                                                                 cAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACa    >  2:918218/1‑100 (MQ=255)
                                                                                                   ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc  >  2:1115760/1‑100 (MQ=255)
                                                                                                   ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc  >  2:837152/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATC  >  NC_000913/1579464‑1579662

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: