Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,059,703 A→T 15.0% intergenic (‑2/‑50) yqfE ← / → argP pseudogene, LysR family/transcriptional regulator for arginine transport and DNA replication genes; replication initiation inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,059,7030AT15.0% 501.4 / 74.2 199intergenic (‑2/‑50)yqfE/argPpseudogene, LysR family/transcriptional regulator for arginine transport and DNA replication genes; replication initiation inhibitor
Reads supporting (aligned to +/- strand):  ref base A (76/93);  new base T (16/14);  total (92/107)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

Only 100 of 205 total aligned reads displayed.

AAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACG  >  NC_000913/3059604‑3059802
                                                                                                   |                                                                                                   
aaaTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACaaa                                                                                                     <  1:904483/100‑1 (MQ=255)
   tGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAAtt                                                                                                  >  1:1637978/1‑100 (MQ=255)
    gCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtg                                                                                                 <  2:3982356/100‑1 (MQ=255)
     cATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgt                                                                                                >  2:1632551/1‑100 (MQ=255)
      aTTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                               >  1:1297058/1‑100 (MQ=255)
        tAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTc                                                                                             <  2:1106625/100‑1 (MQ=255)
         aGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCa                                                                                            >  1:3386376/1‑100 (MQ=255)
          gcgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCan                                                                                           <  2:3352856/100‑2 (MQ=255)
          gcgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCAt                                                                                           >  2:3235149/1‑100 (MQ=255)
            gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGATGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAg                                                                                         >  1:3975726/1‑100 (MQ=255)
            gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggt                                                                                         <  2:3169443/100‑4 (MQ=255)
            gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggt                                                                                         <  2:2734577/100‑4 (MQ=255)
            gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAg                                                                                         <  1:2776532/100‑1 (MQ=255)
             cAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggtt                                                                                        <  1:2799903/100‑5 (MQ=255)
              aaaTCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggttg                                                                                       <  1:3007968/100‑6 (MQ=255)
               aaTCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggttgc                                                                                      <  2:877721/100‑7 (MQ=255)
               aaTCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGc                                                                                      >  1:2803913/1‑100 (MQ=255)
                 tCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAg                                                                                    <  1:3431906/100‑1 (MQ=255)
                  cccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAgg                                                                                   <  2:3698842/100‑1 (MQ=255)
                  cccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAgg                                                                                   >  2:2102481/1‑100 (MQ=255)
                   ccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGa                                                                                  >  2:1015639/1‑100 (MQ=255)
                      gATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaaaa                                                                               >  2:3856548/1‑100 (MQ=255)
                       aTCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaaaa                                                                               >  2:670907/1‑99 (MQ=255)
                        tCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAg                                                                             <  2:3386394/100‑1 (MQ=255)
                        tCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAg                                                                             <  2:229232/100‑1 (MQ=255)
                         cTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGc                                                                            >  2:3242674/1‑100 (MQ=255)
                          ttttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCa                                                                           <  1:2035550/100‑1 (MQ=255)
                           tttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAt                                                                          <  1:2291523/100‑1 (MQ=255)
                                gCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAcc                                                                     >  2:2926106/1‑100 (MQ=255)
                                gCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAcc                                                                     >  2:3500245/1‑100 (MQ=255)
                                gCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACCGCTTCGGTAAAAGCAACCTGACACAAAATTGTGTCAATGTGCAAGAAAAAGCATTTAcc                                                                     <  1:3216168/100‑1 (MQ=255)
                                   gcgcTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAgg                                                                  <  1:2247734/100‑1 (MQ=255)
                                      cTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGc                                                               >  2:3144908/1‑100 (MQ=255)
                                       ttttAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCa                                                              <  1:690977/100‑1 (MQ=255)
                                       ttttAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCa                                                              <  2:785625/100‑1 (MQ=255)
                                          tAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGac                                                           <  2:2743973/100‑1 (MQ=255)
                                          tAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGac                                                           >  1:1624701/1‑100 (MQ=255)
                                          tAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGaa                                                           <  2:1974116/100‑2 (MQ=255)
                                           aacaacGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCAAAAAAAAAAAAAAAGCATTTACCAGGAGCAGaca                                                          <  1:1029896/100‑1 (MQ=255)
                                           aaCCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGaca                                                          <  1:2572467/100‑1 (MQ=255)
                                              cACGGATACAAATCTTACCCACCGACGGCTTCGTTATATGCAACCTGACACAAAATTGTGTCAAAGAGCAAGAAAAAGCATTTACCAGGAGCAGacaac                                                        <  1:819382/99‑1 (MQ=255)
                                                  gATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                   <  1:1996652/100‑1 (MQ=255)
                                                   aTACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAATGAAAGAAAAAGCATTTACCAGGAGCAGACAACAGCAAt                                                  <  1:1547463/100‑1 (MQ=255)
                                                    tACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATg                                                 >  2:2601194/1‑100 (MQ=255)
                                                     aCAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGa                                                >  1:526532/1‑100 (MQ=255)
                                                      cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                               <  2:2639189/100‑1 (MQ=255)
                                                        aaTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTTTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAAc                                             <  2:3267787/100‑1 (MQ=255)
                                                             tACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCg                                        <  1:2983970/100‑1 (MQ=255)
                                                               cccGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGa                                      >  1:1887250/1‑100 (MQ=255)
                                                                ccgccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGAc                                     <  2:2359817/100‑1 (MQ=255)
                                                                ccgccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGAc                                     <  2:1820007/100‑1 (MQ=255)
                                                                 cgccgACGGCATCGGTATATGCAACCTGACACAAAATTGTGGCAGAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACt                                    <  2:3076526/100‑1 (MQ=255)
                                                                 agccgccGGCTTCGGTATATGCAAACTAACAAAAAAATGTGTCAAAGTGCAGAAAAAAGAATTTACAAGGAGCAAACAACAGCAATGAAACGCACgaa                                      >  1:218641/2‑96 (MQ=255)
                                                                  accgACGGCTTCGGGATATGCAAACTGACACAAAATTGTGAAAAAAAGAAAAAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGAc                                     <  1:3069623/97‑1 (MQ=255)
                                                                   ccgACGGCTTCGGTATATGCAACCTGACACAAAATTGCGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTAc                                  <  1:2147296/100‑1 (MQ=255)
                                                                    cgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACa                                 <  2:1119808/100‑1 (MQ=255)
                                                                    cgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACa                                 >  1:3628327/1‑100 (MQ=255)
                                                                     gACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAg                                >  2:4140394/1‑100 (MQ=255)
                                                                      aCGGCTTCGGTATATGCAACCTGACTCAAAATTGGGTCATAGTGCATGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGa                               <  1:1934770/100‑1 (MQ=255)
                                                                      aCGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGa                               <  1:4154122/100‑1 (MQ=255)
                                                                        ggCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAAc                             <  2:3975394/100‑1 (MQ=255)
                                                                        ggCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAAc                             <  1:3882047/100‑1 (MQ=255)
                                                                         gCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGAAAAAAGCATTTACAAAGAGCAGACAACAGCAATGAAACACACGGAATACAAAACa                            >  1:1380283/1‑100 (MQ=255)
                                                                          cTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAAc                             >  1:562984/1‑98 (MQ=255)
                                                                           ttCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACAtt                          <  2:489962/100‑1 (MQ=255)
                                                                           ttCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACAtt                          <  2:2536868/100‑1 (MQ=255)
                                                                             cGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTAc                        <  2:588830/100‑1 (MQ=255)
                                                                             cGGTATATGCAACCTGACACAAAATTGTGTCATAATACAAGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTa                         <  1:3460521/99‑1 (MQ=255)
                                                                              ggTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACa                       >  1:1202315/1‑100 (MQ=255)
                                                                              ggTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACa                       <  2:1419487/100‑1 (MQ=255)
                                                                               gTATATGCAACCCGACACAAAAGTGTGTCATAGTGCAGGAAAAAGCATTAACCAGGAGCAGACAACAGCAATGAAACGCCAGGACTACAGAACATTACAg                      <  1:4039623/100‑1 (MQ=255)
                                                                                tataTGCAACCTGACACAAAATTGTGTCATAGTGAAGAAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAgg                     <  1:3140643/100‑1 (MQ=255)
                                                                                   atGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCAc                  >  2:59506/1‑100 (MQ=255)
                                                                                    tGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAgagatc                 >  1:1735112/1‑95 (MQ=255)
                                                                                    tGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACt                 >  1:1568128/1‑100 (MQ=255)
                                                                                    tGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACt                 >  2:1368363/1‑100 (MQ=255)
                                                                                    tGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACt                 >  1:871500/1‑100 (MQ=255)
                                                                                    tGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGAAATGAAAAGCCCGCACTACAGAACATTAAAGGCAc                  >  1:4382560/1‑99 (MQ=255)
                                                                                     gCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTg                >  1:918624/1‑100 (MQ=255)
                                                                                      cactatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTgg               >  1:180417/6‑100 (MQ=255)
                                                                                      cAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTgg               <  1:2102473/100‑1 (MQ=255)
                                                                                       actatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGa              >  2:3624874/5‑100 (MQ=255)
                                                                                       actatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGa              >  2:4058117/5‑100 (MQ=255)
                                                                                       aaCCTGACACAAAATTGTGACATAGTGCAGGAAAAAGCATTCACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTgg               >  1:270188/1‑99 (MQ=255)
                                                                                        ctatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGAt             >  2:710586/4‑100 (MQ=255)
                                                                                        ctatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGAt             >  1:2254209/4‑100 (MQ=255)
                                                                                        aCCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGAt             >  2:121106/1‑100 (MQ=255)
                                                                                         ccTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATg            >  1:4296524/1‑100 (MQ=255)
                                                                                         ccTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATg            <  1:77952/100‑1 (MQ=255)
                                                                                          atGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATg            >  1:4126609/2‑99 (MQ=255)
                                                                                          atGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGc           >  2:1891296/2‑100 (MQ=255)
                                                                                           tGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCg          >  2:2235280/1‑100 (MQ=255)
                                                                                           tGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCg          <  1:3150111/100‑1 (MQ=255)
                                                                                            gACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCgg         <  2:4203854/100‑1 (MQ=255)
                                                                                              cacaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg       <  2:3024839/100‑1 (MQ=255)
                                                                                               acaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGtca      >  1:2675542/1‑98 (MQ=255)
                                                                                               acaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGa      >  2:4263259/1‑100 (MQ=255)
                                                                                               acaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGa      <  1:385118/100‑1 (MQ=255)
                                                                                               acaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGa      <  2:2817627/100‑1 (MQ=255)
                                                                                                caAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGAt     >  1:3919476/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACG  >  NC_000913/3059604‑3059802

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: