Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,579,570 A→C 17.9% intergenic (‑228/+63) 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,5700AC17.9% 246.7 / 46.3 104intergenic (‑228/+63)yddA/ydeMputative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base A (35/50);  new base C (8/11);  total (43/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

ATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCA  >  NC_000913/1579471‑1579664
                                                                                                   |                                                                                              
aTTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACGCCCGACAGGAATTTa                                                                                                >  1:2520104/1‑100 (MQ=255)
 ttGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTAt                                                                                               <  1:876083/100‑1 (MQ=255)
  tGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATg                                                                                              >  2:1253369/1‑100 (MQ=255)
  tGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATg                                                                                              <  2:2214751/100‑1 (MQ=255)
   gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGt                                                                                             >  2:22322/1‑100 (MQ=255)
   gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGt                                                                                             >  2:1675281/1‑100 (MQ=255)
     tAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCg                                                                                           <  1:1060245/100‑1 (MQ=255)
      aaTCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCgg                                                                                          >  1:1617816/1‑100 (MQ=255)
       aTCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACCCCCGACAGGAATTTATGTCggg                                                                                         <  2:1914097/100‑1 (MQ=255)
        tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                        <  1:2198689/100‑1 (MQ=255)
        tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                        <  2:528287/100‑1 (MQ=255)
        tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                        <  1:2141916/100‑1 (MQ=255)
        tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                        <  2:1924648/100‑1 (MQ=255)
        tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                        <  2:490565/100‑1 (MQ=255)
        tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                        <  1:2287845/100‑1 (MQ=255)
        tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGt                                                                                        <  1:2218127/100‑1 (MQ=255)
        tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGt                                                                                        <  1:745820/100‑1 (MQ=255)
         cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATTTCGGGTg                                                                                       <  1:1350305/100‑1 (MQ=255)
         cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTg                                                                                       <  1:1843814/100‑1 (MQ=255)
         cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTg                                                                                       <  1:552619/100‑1 (MQ=255)
          aaGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTgt                                                                                      <  1:2345974/100‑2 (MQ=38)
          aaGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTgt                                                                                      <  1:2345961/100‑2 (MQ=38)
          aaGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGc                                                                                      >  2:659872/1‑100 (MQ=255)
             gATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCtt                                                                                   <  1:2036645/100‑1 (MQ=255)
               tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGt                                                                                 <  1:765410/100‑1 (MQ=255)
               tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGt                                                                                 >  1:202952/1‑100 (MQ=255)
               tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGt                                                                                 <  2:804566/100‑1 (MQ=255)
                 ttGCTAATGCCCTGATGCATACCATGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTgtttcgt                                                                                 <  2:1115471/98‑7 (MQ=25)
                 ttGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGtt                                                                                >  2:2216678/1‑99 (MQ=255)
                  tGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTaa                                                                              >  1:1968000/1‑100 (MQ=255)
                    cTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                            >  2:1640492/1‑100 (MQ=255)
                    cTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                            <  1:1611693/100‑1 (MQ=255)
                    cTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                            <  2:378384/100‑1 (MQ=255)
                     tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                           <  1:898051/100‑1 (MQ=255)
                     tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                           <  1:39030/100‑1 (MQ=255)
                        tGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCAt                                                                        >  2:1512682/1‑100 (MQ=255)
                          cccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaa                                                                      <  2:133761/100‑1 (MQ=255)
                          cccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaa                                                                      >  1:1871889/1‑100 (MQ=255)
                           ccTGATGCATACCACGTAAGCCAGGATTTTCGCCAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaag                                                                     >  2:1595778/1‑100 (MQ=255)
                           ccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaa                                                                      >  1:841815/1‑99 (MQ=255)
                              gATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaagaag                                                                  >  2:1278213/1‑100 (MQ=255)
                              gATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaagaag                                                                  <  1:1540485/100‑1 (MQ=255)
                              gATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaagaa                                                                   >  2:67114/1‑99 (MQ=255)
                              gATGCATACCACGTAAGCCAGGATTTTCACAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaaga                                                                    <  2:510274/98‑1 (MQ=255)
                                tGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAagga                                                                >  1:1934574/1‑100 (MQ=255)
                                 gCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggag                                                               <  1:969098/100‑1 (MQ=255)
                                  cATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggagg                                                              >  2:599054/1‑100 (MQ=255)
                                   aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                             <  2:1301681/100‑1 (MQ=255)
                                   aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                             <  1:920167/100‑1 (MQ=255)
                                   aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                             <  2:649327/100‑1 (MQ=255)
                                   aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                             <  1:22323/100‑1 (MQ=255)
                                     aCCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                             <  1:2179193/98‑1 (MQ=255)
                                      ccACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTaa                                                          >  2:497116/1‑100 (MQ=255)
                                       cACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAg                                                         <  1:2085922/100‑1 (MQ=255)
                                        aCGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGa                                                        >  1:1140854/1‑100 (MQ=255)
                                            aaGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATgg                                                    <  1:2214384/100‑1 (MQ=255)
                                              gCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                  <  2:1779916/100‑1 (MQ=255)
                                              gCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                  <  1:1085031/100‑1 (MQ=255)
                                              gCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                  <  1:1912399/100‑1 (MQ=255)
                                               ccAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGt                                                 >  1:2157531/1‑100 (MQ=255)
                                                   gATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATCAGAAGGAGGCTAAGAATGGAGTTaaa                                             <  1:1561415/100‑1 (MQ=255)
                                                   gATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaa                                             >  1:1579953/1‑100 (MQ=255)
                                                    aTTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCGTGTTAAGGTCACAAGAAGGAGGCTAAGAATGGAGTTaaaa                                            <  2:515883/100‑1 (MQ=255)
                                                        tCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagaga                                        >  2:320852/1‑100 (MQ=255)
                                                         cGCAAGGGAAGATGATGATGACGAAAAACCAGACAGGAATTTATGTAGGGTGCCTAGTTAAGGTCATGAGAAGGTGGATAAGAATGGAATTAAAagagag                                       >  2:226237/1‑100 (MQ=255)
                                                           cAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCg                                     >  1:1753031/1‑100 (MQ=255)
                                                             atggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGTTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGtt                                   <  2:2507624/98‑1 (MQ=255)
                                                              gggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGtta                                  <  2:2097408/100‑1 (MQ=255)
                                                                gAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttatt                                <  1:196701/100‑1 (MQ=255)
                                                                     tgatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatta                           >  1:519389/1‑100 (MQ=255)
                                                                          atgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtctc                      >  1:24886/1‑100 (MQ=255)
                                                                           tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtctcn                     <  1:1275631/100‑2 (MQ=255)
                                                                           tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTcc                     >  1:1877749/1‑100 (MQ=255)
                                                                           tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTcc                     >  2:1245328/1‑100 (MQ=255)
                                                                             acaaggcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAt                   >  2:453148/7‑100 (MQ=37)
                                                                             acaaggcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAt                   >  2:2305468/7‑100 (MQ=37)
                                                                             aCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAt                   <  2:1033793/100‑1 (MQ=255)
                                                                             aCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAt                   <  2:2157517/100‑1 (MQ=255)
                                                                              cGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAtt                  >  2:322701/1‑100 (MQ=255)
                                                                                aaaCACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTg                <  1:1245336/100‑1 (MQ=255)
                                                                                aaaCACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTg                >  1:1806636/1‑100 (MQ=255)
                                                                                  gcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtt              >  2:2503660/2‑100 (MQ=38)
                                                                                  gcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtt              >  2:2549832/2‑100 (MQ=38)
                                                                                  gcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtt              >  2:2275948/2‑100 (MQ=38)
                                                                                   nacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGCCAAAAGAAGGCGGCTAAGAATGGAGTTAAAAGAGAGCTTTATTACTTATTCTACATTTGttt             >  2:790471/2‑100 (MQ=255)
                                                                                   cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt             >  1:1987630/1‑100 (MQ=38)
                                                                                   cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt             >  2:1442208/1‑100 (MQ=38)
                                                                                    acCCTACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGTATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtttt            <  2:883278/100‑1 (MQ=255)
                                                                                     cccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATTGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtttt            >  2:16082/1‑99 (MQ=255)
                                                                                      ccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGc          <  2:800481/100‑1 (MQ=255)
                                                                                       cGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCa         >  2:2333834/1‑100 (MQ=255)
                                                                                        gACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCaa        <  1:1150400/100‑1 (MQ=255)
                                                                                        gACAGGAATTTATGTCAGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAAGTATTCTCCATTTGTTTTGCaa        <  2:1317107/100‑1 (MQ=255)
                                                                                         aCAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCaa        <  2:2463681/99‑1 (MQ=255)
                                                                                           aGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACAt     <  2:545770/100‑1 (MQ=255)
                                                                                           aGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACAt     <  1:1978670/100‑1 (MQ=255)
                                                                                            ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACAt     >  2:257271/1‑99 (MQ=255)
                                                                                            ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc    <  1:1601794/100‑1 (MQ=255)
                                                                                             gAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATcc   <  1:497122/100‑1 (MQ=255)
                                                                                              aaTTTATGTCGGGTGCCTTGTTAATGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCa  <  2:1885731/100‑1 (MQ=255)
                                                                                                   |                                                                                              
ATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCA  >  NC_000913/1579471‑1579664

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: