Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,059,703 A→T 16.4% 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,7030AT16.4% 434.2 / 72.4 177intergenic (‑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 (72/76);  new base T (12/17);  total (84/93)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

Only 100 of 185 total aligned reads displayed.

AATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACG  >  NC_000913/3059605‑3059802
                                                                                                  |                                                                                                   
aaTGCATTAGCGCAAATCCCTGATATTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACaaaa                                                                                                    <  1:1779675/100‑1 (MQ=255)
 aTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAAt                                                                                                   <  1:1434663/100‑1 (MQ=255)
  tGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAAtt                                                                                                  >  2:3756288/1‑100 (MQ=255)
     ataaccgcAAATCCCTGATCTTTTAGGCTGCGCTGTTAACCAGGGCTTCAAATGTTACTCGCCGACGGCTTCGATATATACAACCTGACACAAAATtgtg                                                                                               >  2:2932127/6‑100 (MQ=255)
      ttAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTgtgt                                                                                              <  1:257218/100‑1 (MQ=255)
        aGCGCAAATCCCTGATCTTTTAAGCGGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCa                                                                                            <  1:747006/100‑1 (MQ=255)
         gcgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCag                                                                                           <  1:627523/100‑2 (MQ=255)
          cgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCagg                                                                                          <  2:3154496/100‑3 (MQ=255)
          cgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCagg                                                                                          <  2:1057221/100‑3 (MQ=255)
          cgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCagg                                                                                          <  2:1182268/100‑3 (MQ=255)
          cgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCagg                                                                                          <  2:1610031/100‑3 (MQ=255)
           gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggt                                                                                         <  1:1075866/100‑4 (MQ=255)
           gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggt                                                                                         <  1:2637377/100‑4 (MQ=255)
           gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAg                                                                                         >  1:1614955/1‑100 (MQ=255)
            cAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAg                                                                                         <  2:2128531/99‑1 (MQ=255)
             aaaTCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGt                                                                                        >  2:1387182/1‑99 (MQ=255)
               aTCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCa                                                                                     <  1:2238999/100‑1 (MQ=255)
                tCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAg                                                                                    <  2:2858541/100‑1 (MQ=255)
                tCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAg                                                                                    <  1:2876785/100‑1 (MQ=255)
                tCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGACACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTg                                                                                       <  2:1704861/97‑1 (MQ=255)
                 cccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAgg                                                                                   >  1:587816/1‑100 (MQ=255)
                  ccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGa                                                                                  >  2:2907478/1‑100 (MQ=255)
                  ccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGa                                                                                  <  1:1320324/100‑1 (MQ=255)
                   cTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaa                                                                                 >  1:3727438/1‑100 (MQ=255)
                   cTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaa                                                                                 >  1:304568/1‑100 (MQ=255)
                   cTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaa                                                                                 >  1:876876/1‑100 (MQ=255)
                     gATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaaaa                                                                               <  2:2240611/100‑1 (MQ=255)
                       tCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAg                                                                             <  1:2017365/100‑1 (MQ=255)
                       tCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAg                                                                             <  2:1145023/100‑1 (MQ=255)
                        cTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGc                                                                            >  2:1404791/1‑100 (MQ=255)
                        cTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGc                                                                            <  1:1976666/100‑1 (MQ=255)
                             aaGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTa                                                                       <  2:3794869/100‑1 (MQ=255)
                             aaGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTa                                                                       >  2:1476366/1‑100 (MQ=255)
                              aGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAc                                                                      >  2:795438/1‑100 (MQ=255)
                               gCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAcc                                                                     <  1:3015378/100‑1 (MQ=255)
                               gCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAcc                                                                     >  2:3364622/1‑100 (MQ=255)
                                cTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                    >  1:1929956/1‑100 (MQ=255)
                                cTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                    <  1:1266499/100‑1 (MQ=255)
                                 tGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCGACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAg                                                                   >  1:2505328/1‑100 (MQ=255)
                                 tGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAcc                                                                     <  2:487199/98‑1 (MQ=255)
                                 tGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAg                                                                   >  1:3406508/1‑100 (MQ=255)
                                   cgcTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGa                                                                 <  1:89023/100‑1 (MQ=255)
                                      ttttAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCa                                                              <  1:1491223/100‑1 (MQ=255)
                                       tttAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAg                                                             <  1:567828/100‑1 (MQ=255)
                                         tAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGac                                                           >  1:3718886/1‑100 (MQ=255)
                                           aCCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACATGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaa                                                         >  2:1672073/1‑100 (MQ=255)
                                             cACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaaca                                                       <  2:3667025/100‑1 (MQ=255)
                                              aCGGCTACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAg                                                      <  2:1970940/100‑1 (MQ=255)
                                              aCGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAg                                                      <  1:2532651/100‑1 (MQ=255)
                                               cGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGc                                                     >  2:2530778/1‑100 (MQ=255)
                                                ggATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAg                                                      <  2:2314400/98‑1 (MQ=255)
                                                ggATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCa                                                    <  2:556749/100‑1 (MQ=255)
                                                 gATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                   >  2:2263985/1‑100 (MQ=255)
                                                  aTACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAAt                                                  <  1:166114/100‑1 (MQ=255)
                                                     cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                               >  1:2947074/1‑100 (MQ=255)
                                                     cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                               >  1:3290734/1‑100 (MQ=255)
                                                     cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                               >  2:3345949/1‑100 (MQ=255)
                                                     cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                               >  2:1906850/1‑100 (MQ=255)
                                                     cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                               >  1:3086131/1‑100 (MQ=255)
                                                     cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                               >  1:1112527/1‑100 (MQ=255)
                                                           ttACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGccc                                         <  2:3135359/100‑1 (MQ=255)
                                                              cccGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGa                                      >  1:1534741/1‑100 (MQ=255)
                                                              cccGCCGACGGCTTCGGGATATGCAACCTGACACAAAATTGCGTCATAGTGCAGGAAAAAGGATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGa                                      >  2:2072658/1‑100 (MQ=255)
                                                               ccgccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGAc                                     <  1:3778131/100‑1 (MQ=255)
                                                                 gccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTa                                   >  1:829490/1‑100 (MQ=255)
                                                                       ggCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAAc                             <  1:2221032/100‑1 (MQ=255)
                                                                        gCTTCGGTATATGCAACCTGACACAAAGTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACa                            >  2:2893098/1‑100 (MQ=255)
                                                                        gCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACa                            >  1:1867608/1‑100 (MQ=255)
                                                                        gCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACa                            >  2:1056853/1‑100 (MQ=255)
                                                                        ccctCGGGATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAATAACTCACCAGGAGCAGAAAACAGCAATCAAACGCCCGGCCAACAGaaaa                            >  2:1500420/4‑98 (MQ=255)
                                                                         cTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACAt                           >  1:104763/1‑100 (MQ=255)
                                                                         cTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACAt                           >  1:2479424/1‑100 (MQ=255)
                                                                                atatGCAAACTGACACAAAATTGTGTCATAGTGCAGTAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAg                      <  2:2209132/98‑1 (MQ=255)
                                                                                   tgcactatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACt                 >  2:3291810/8‑100 (MQ=255)
                                                                                   tGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACt                 >  1:3127983/1‑100 (MQ=255)
                                                                                   tGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAACAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATGACAGGCACt                 >  2:2434380/1‑100 (MQ=255)
                                                                                    gCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTg                <  2:1877985/100‑1 (MQ=255)
                                                                                     cactatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTAACAGGAGCAGACAACAGCAATGAAACGCCTGGACTACAGAAAATTACAGGCACTgg               >  2:3025614/6‑100 (MQ=255)
                                                                                     cAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCTGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTgg               >  2:2117011/1‑100 (MQ=255)
                                                                                       ctatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGAt             >  2:320505/4‑100 (MQ=255)
                                                                                       ctatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGAt             >  1:3567289/4‑100 (MQ=255)
                                                                                       ctatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGAt             >  1:954701/4‑100 (MQ=255)
                                                                                       ctatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAAGCACGGGAt             >  2:668374/4‑100 (MQ=255)
                                                                                        ccTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATg            <  2:729895/100‑1 (MQ=255)
                                                                                         cTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGc           >  2:2763539/1‑100 (MQ=255)
                                                                                         atGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGc           >  1:3315286/2‑100 (MQ=255)
                                                                                          tGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACGAGGAGCATACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCg          >  2:1270523/1‑100 (MQ=255)
                                                                                           gACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCgg         <  1:3345974/100‑1 (MQ=255)
                                                                                            acacAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGt        <  1:2829885/100‑1 (MQ=255)
                                                                                            acacAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGt        <  1:3462550/100‑1 (MQ=255)
                                                                                            acacAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGt        <  1:1559450/100‑1 (MQ=255)
                                                                                             cacaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg       <  2:1244931/100‑1 (MQ=255)
                                                                                             cacaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg       <  1:757533/100‑1 (MQ=255)
                                                                                             cacaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg       >  2:3652688/1‑100 (MQ=255)
                                                                                             cacaAAATTGTGTCATAGTGCAGGAAAAAGCACTTACCAGGAGCAGACAACAGCAACGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCg          >  2:3238050/1‑97 (MQ=255)
                                                                                              acaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTgg      >  1:780230/1‑99 (MQ=255)
                                                                                                aaaaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATa    >  1:3485783/1‑100 (MQ=255)
                                                                                                 aaaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATAc   <  1:2204499/100‑1 (MQ=255)
                                                                                                  aaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACg  <  1:2248015/100‑1 (MQ=255)
                                                                                                  aaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACg  <  1:506967/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
AATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACG  >  NC_000913/3059605‑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: