Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,059,703 A→T 20.6% 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,7030AT20.6% 395.4 / 106.8 191intergenic (‑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/75);  new base T (19/20);  total (96/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-01

Only 100 of 198 total aligned reads displayed.

AAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACG  >  NC_000913/3059604‑3059802
                                                                                                   |                                                                                                   
aaaTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACaaa                                                                                                     <  1:2673827/100‑1 (MQ=255)
aaaTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACaaa                                                                                                     <  1:493286/100‑1 (MQ=255)
 aaTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACaaaa                                                                                                    <  1:100144/100‑1 (MQ=255)
  aTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAAt                                                                                                   <  2:1405345/100‑1 (MQ=255)
    gCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtg                                                                                                 <  1:298903/100‑1 (MQ=255)
     cATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgt                                                                                                >  2:3256310/1‑100 (MQ=255)
     cATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgt                                                                                                >  2:3737520/1‑100 (MQ=255)
      aTTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                               <  1:2682861/100‑1 (MQ=255)
      aTTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                               <  1:4133173/100‑1 (MQ=255)
       ttAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCTCGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTgtgt                                                                                              <  1:3987269/100‑1 (MQ=255)
       ttAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTgtgt                                                                                              <  2:3100107/100‑1 (MQ=255)
        tAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTc                                                                                             <  1:262933/100‑1 (MQ=255)
          gcgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCag                                                                                           <  1:1992973/100‑2 (MQ=255)
          gcgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCag                                                                                           <  2:1936986/100‑2 (MQ=255)
          gcgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCag                                                                                           <  1:1225215/100‑2 (MQ=255)
          gcgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCAt                                                                                           <  2:2389565/100‑1 (MQ=255)
          gcgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCAt                                                                                           <  2:817705/100‑1 (MQ=255)
           ctcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCagg                                                                                          <  2:3226309/98‑3 (MQ=255)
           cgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGTATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCagg                                                                                          <  2:957875/100‑3 (MQ=255)
           cgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTGCCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCagg                                                                                          <  2:2262608/100‑3 (MQ=255)
           cgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGTTTCGGTATATGCAACCTGACACAATTTTGTGTCagg                                                                                          <  1:2962812/100‑3 (MQ=255)
           cgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATa                                                                                          <  2:1914091/100‑1 (MQ=255)
           cgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATa                                                                                          >  1:1671313/1‑100 (MQ=255)
            gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGTATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggt                                                                                         <  2:1062226/100‑4 (MQ=255)
            gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggt                                                                                         <  2:261844/100‑4 (MQ=255)
             cAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTgtgtaaggag                                                                                        >  1:1029713/1‑94 (MQ=255)
             cAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggag                                                                                        >  1:3226323/1‑96 (MQ=255)
                  cccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAgg                                                                                   >  2:646705/1‑100 (MQ=255)
                  cccTGATCTTTTAAGCTGCGCTTTTAACCACGGAGAGAAATGTTACGCGCCGACGGCTTCGGTATATGCAAACTGACACAAAATTGTGTCATAGTGCAgg                                                                                   >  2:1055838/1‑100 (MQ=255)
                     tGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaaa                                                                                >  2:2404503/1‑100 (MQ=255)
                      gATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaaaa                                                                               >  1:2189541/1‑100 (MQ=255)
                      gATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaaaa                                                                               >  2:3819764/1‑100 (MQ=255)
                        tCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAg                                                                             <  1:2277836/100‑1 (MQ=255)
                        tCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAg                                                                             <  1:472810/100‑1 (MQ=255)
                         cTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGc                                                                            >  1:591535/1‑100 (MQ=255)
                            ttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAtt                                                                         <  2:917366/100‑1 (MQ=255)
                             tAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAttt                                                                        >  1:3043832/1‑100 (MQ=255)
                             tAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAttt                                                                        >  1:512227/1‑100 (MQ=255)
                              aaGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTa                                                                       <  1:3167111/100‑1 (MQ=255)
                                 cTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                    <  1:3738436/100‑1 (MQ=255)
                                  tGCGCTTTTAACCACGGATACAAATGTTACGCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAg                                                                   <  2:2567848/100‑1 (MQ=255)
                                  tGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAg                                                                   <  2:1640891/100‑1 (MQ=255)
                                  tGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGAAAAAAAAAAAAAAAAAGCATTTACCAg                                                                   <  2:4132712/100‑1 (MQ=255)
                                   gcgcTTTTACCCACGGATACAAATGTTCCCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGTAAAAAGCATTTACCAgg                                                                  <  2:1357128/100‑1 (MQ=255)
                                   gcgcTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAgg                                                                  <  2:2213899/100‑1 (MQ=255)
                                    cgcTTTTAACCACGGATACAAATGTTACCCGCCGACGACTTCGCCACATGCAACCAGACACAAAATGCTGCCATTGTGCAGGAAGAAGCATTTACCAGGa                                                                 >  2:2106793/1‑100 (MQ=255)
                                      cTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGc                                                               >  1:1745868/1‑100 (MQ=255)
                                      cTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGc                                                               <  1:2590913/100‑1 (MQ=255)
                                      cTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGc                                                               >  2:223601/1‑100 (MQ=255)
                                      cTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGc                                                               >  1:2885276/1‑100 (MQ=255)
                                       ttttAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCa                                                              <  2:1198251/100‑1 (MQ=255)
                                        tttAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGc                                                               <  1:1271728/98‑1 (MQ=255)
                                            aCCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaa                                                         >  1:959782/1‑100 (MQ=255)
                                             ccACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaac                                                        >  1:2620184/1‑100 (MQ=255)
                                                cGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGc                                                     >  1:1493617/1‑100 (MQ=255)
                                                  gATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                   <  2:2712494/100‑1 (MQ=255)
                                                     aCAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGa                                                <  1:4185278/100‑1 (MQ=255)
                                                      cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                               <  2:3145329/100‑1 (MQ=255)
                                                      cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                               >  1:2532603/1‑100 (MQ=255)
                                                      cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAGTTACCAGGAGGCGACAACAGCAATGaa                                               >  2:1232119/1‑100 (MQ=255)
                                                       aaaTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaaa                                              <  1:1144858/100‑1 (MQ=255)
                                                          tGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACg                                            <  1:2732494/99‑1 (MQ=255)
                                                           gTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGcc                                          >  2:3473802/1‑100 (MQ=255)
                                                              aCCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCgg                                       <  2:1555365/100‑1 (MQ=255)
                                                               cccGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGa                                      >  1:3352693/1‑100 (MQ=255)
                                                                ccgccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGAc                                     >  1:659298/1‑100 (MQ=255)
                                                                 cgccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACt                                    <  2:1900562/100‑1 (MQ=255)
                                                                         gCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACa                            >  1:162455/1‑100 (MQ=255)
                                                                            tCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTa                         >  1:3689049/1‑100 (MQ=255)
                                                                              ggTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACa                       >  1:2238731/1‑100 (MQ=255)
                                                                              ggTATATGCAACCTGACACAAAATTGTGTCATAGTGCAAGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACa                       <  2:1694794/100‑1 (MQ=255)
                                                                              ggTATATGCAACCCGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACa                       >  1:4097083/1‑100 (MQ=255)
                                                                               gTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAg                      <  2:39647/100‑1 (MQ=255)
                                                                               gTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAg                      <  2:696508/100‑1 (MQ=255)
                                                                                 tcctgcactatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGc                    >  1:3577434/11‑100 (MQ=38)
                                                                                  tatGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCa                   >  1:3695004/1‑100 (MQ=255)
                                                                                  tatGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCa                   >  2:4113992/1‑100 (MQ=255)
                                                                                  cctgcactatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCa                   >  2:603139/10‑100 (MQ=38)
                                                                                    tgcactatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACt                 >  1:2424871/8‑100 (MQ=255)
                                                                                    tGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACACTACAGGCACt                 >  2:2743961/1‑100 (MQ=255)
                                                                                      cactatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACTAGGAGCAGACAACAGCAATGAAACGACCGAACTAAAGAACATTACAGGCACggg               >  2:4196597/6‑100 (MQ=255)
                                                                                      cAACCTGACACAAAATTGTGTCATAGTGCAGGTAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTgg               <  2:2939826/100‑1 (MQ=255)
                                                                                      cAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATGTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTgg               >  2:1798569/1‑100 (MQ=255)
                                                                                       actatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGAACTGGa              >  2:1030275/5‑100 (MQ=255)
                                                                                       aaCCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGa              <  2:1465114/100‑1 (MQ=255)
                                                                                          cTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGc           <  2:1493614/100‑1 (MQ=255)
                                                                                          atGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGc           >  1:3119733/2‑100 (MQ=255)
                                                                                            gACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCgg         >  1:4068476/1‑100 (MQ=255)
                                                                                            gACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCgg         >  2:3442329/1‑100 (MQ=255)
                                                                                             acacAAACTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGt        <  2:2000507/100‑1 (MQ=255)
                                                                                             acacAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGt        >  2:2356214/1‑100 (MQ=255)
                                                                                              cacaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg       >  1:3741816/1‑100 (MQ=255)
                                                                                              cacaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg       >  1:3469766/1‑100 (MQ=255)
                                                                                              cacaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACCGGATGCGGTg       >  2:2764897/1‑100 (MQ=255)
                                                                                               acaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGCTGa      >  2:1261110/1‑100 (MQ=255)
                                                                                                caAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGtaga     >  1:2403023/1‑97 (MQ=255)
                                                                                                caAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGAt     >  1:1378475/1‑100 (MQ=255)
                                                                                                 aaaaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATa    <  1:801391/100‑1 (MQ=255)
                                                                                                  aaaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATa    <  2:3658576/99‑1 (MQ=255)
                                                                                                   aaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACg  >  2:3182715/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: