Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,059,704 A→T 16.4% intergenic (‑3/‑49) 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,7040AT16.4% 431.1 / 71.4 177intergenic (‑3/‑49)yqfE/argPpseudogene, LysR family/transcriptional regulator for arginine transport and DNA replication genes; replication initiation inhibitor
Reads supporting (aligned to +/- strand):  ref base A (71/76);  new base T (12/17);  total (84/93)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 187 total aligned reads displayed.

AATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGT  >  NC_000913/3059605‑3059803
                                                                                                   |                                                                                                   
aaTGCATTAGCGCAAATCCCTGATATTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACaaaa                                                                                                     <  1:1779675/100‑1 (MQ=255)
 acgCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAAt                                                                                                    <  2:2318059/98‑1 (MQ=255)
 aTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAAt                                                                                                    <  1:1434663/100‑1 (MQ=255)
     ataaccgcAAATCCCTGATCTTTTAGGCTGCGCTGTTAACCAGGGCTTCAAATGTTACTCGCCGACGGCTTCGATATATACAACCTGACACAAAATtgtg                                                                                                >  2:2932127/6‑100 (MQ=255)
     aTTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                <  1:605994/100‑1 (MQ=255)
     aTTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                <  1:2516393/100‑1 (MQ=255)
       tAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTc                                                                                              >  1:3675156/1‑100 (MQ=255)
         gcgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCac                                                                                            <  1:3501318/100‑2 (MQ=255)
          cgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCagg                                                                                           <  2:1057221/100‑3 (MQ=255)
          cgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCagg                                                                                           <  2:1610031/100‑3 (MQ=255)
           gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACTACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGCGTCaggt                                                                                          <  2:2000333/100‑4 (MQ=255)
           gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggt                                                                                          <  2:2464130/100‑4 (MQ=255)
           gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggt                                                                                          <  1:1075866/100‑4 (MQ=255)
           gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggt                                                                                          <  1:2637377/100‑4 (MQ=255)
           gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAg                                                                                          >  1:3098151/1‑100 (MQ=255)
           gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAg                                                                                          >  1:1614955/1‑100 (MQ=255)
            cAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggtt                                                                                         <  1:934424/100‑5 (MQ=255)
            cAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggta                                                                                         >  2:1255861/1‑96 (MQ=255)
            cAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggt                                                                                          <  2:2609715/99‑4 (MQ=255)
            cAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAg                                                                                          <  2:2128531/99‑1 (MQ=255)
            cAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGt                                                                                         >  2:2075299/1‑100 (MQ=255)
              aaTCCCTGATCTGTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCCTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGc                                                                                       <  1:567771/100‑1 (MQ=255)
               aTCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCa                                                                                      <  1:2238999/100‑1 (MQ=255)
                tCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggttgcat                                                                                     <  1:3392022/100‑9 (MQ=38)
                tCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAg                                                                                     <  1:2876785/100‑1 (MQ=255)
                tCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGACACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTg                                                                                        <  2:1704861/97‑1 (MQ=255)
                 cccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAgg                                                                                    >  2:3153903/1‑100 (MQ=255)
                  ccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGa                                                                                   >  2:140400/1‑100 (MQ=255)
                  ccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGa                                                                                   <  1:1320324/100‑1 (MQ=255)
                  ccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGa                                                                                   >  2:2907478/1‑100 (MQ=255)
                   cTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaa                                                                                  >  1:3727438/1‑100 (MQ=255)
                     gATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaaaa                                                                                <  2:2240611/100‑1 (MQ=255)
                        cTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGc                                                                             <  1:576615/100‑1 (MQ=255)
                         ttttAAGCTGCGCTTTTAACCACGGATACAAATGTTATCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGc                                                                             <  2:1090092/99‑1 (MQ=255)
                         ttttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCa                                                                            >  1:1981361/1‑100 (MQ=255)
                          tttAAGCTGCGCTTGTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAt                                                                           <  2:3181301/100‑1 (MQ=255)
                             aaGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTa                                                                        >  2:1476366/1‑100 (MQ=255)
                             aaGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTa                                                                        <  2:3794869/100‑1 (MQ=255)
                              aGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAc                                                                       >  2:795438/1‑100 (MQ=255)
                               gCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAcc                                                                      >  2:3364622/1‑100 (MQ=255)
                                cTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                     >  1:1929956/1‑100 (MQ=255)
                                cTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                     >  1:2544400/1‑100 (MQ=255)
                                cTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                     <  1:1266499/100‑1 (MQ=255)
                                cTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                     >  1:889547/1‑100 (MQ=255)
                                 tGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAcc                                                                      <  2:487199/98‑1 (MQ=255)
                                   cgcTTTTAACCACGGATACAAATGTTACCCGCGGACGGCTTCGGTATATGCAACCTGAAACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                     <  2:2307470/97‑1 (MQ=255)
                                      ttttAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCa                                                               <  1:1491223/100‑1 (MQ=255)
                                       tttAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAg                                                              <  2:1163226/100‑1 (MQ=255)
                                       tttAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAg                                                              <  1:466784/100‑1 (MQ=255)
                                         tAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGac                                                            >  1:3718886/1‑100 (MQ=255)
                                           aCCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACATGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaa                                                          >  2:1672073/1‑100 (MQ=255)
                                            ccACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaac                                                         <  1:833805/100‑1 (MQ=255)
                                             cACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaaca                                                        <  2:3667025/100‑1 (MQ=255)
                                               cGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGAAGAAAAAAGCATTTACCAGGAGCAGACAACAGc                                                      <  2:1211182/100‑1 (MQ=255)
                                                ggATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCa                                                     <  2:556749/100‑1 (MQ=255)
                                                  aTACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAAt                                                   <  1:166114/100‑1 (MQ=255)
                                                   tACCAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGACCAGACAACAGCAATg                                                  <  2:1048198/100‑1 (MQ=255)
                                                     cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                                >  1:1112527/1‑100 (MQ=255)
                                                     cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                                >  1:3290734/1‑100 (MQ=255)
                                                     cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                                >  2:3345949/1‑100 (MQ=255)
                                                        aTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACg                                             >  1:2962747/1‑100 (MQ=255)
                                                           ttACGCGCCGACGGCTGCGGTATATGCAACCGGACACAAAATGGTGTGATAGTGCAGGAAAAACCATTTACCAGGAGCAGACAACAGCAATGAAACGccc                                          <  2:3746817/100‑1 (MQ=255)
                                                           ttACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGccc                                          <  2:3135359/100‑1 (MQ=255)
                                                                cgccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACt                                     >  1:3035128/1‑100 (MQ=255)
                                                                   cgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTAc                                   <  2:1567257/99‑1 (MQ=255)
                                                                    gACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAg                                 <  1:3070890/100‑1 (MQ=255)
                                                                       ggCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAAc                              <  1:2221032/100‑1 (MQ=255)
                                                                        gCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACa                             >  1:1867608/1‑100 (MQ=255)
                                                                        ccctCGGGATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAATAACTCACCAGGAGCAGAAAACAGCAATCAAACGCCCGGCCAACAGaaaa                             >  2:1500420/4‑98 (MQ=255)
                                                                         cTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACAt                            >  2:2891829/1‑100 (MQ=255)
                                                                         cTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACAt                            >  1:104763/1‑100 (MQ=255)
                                                                         cTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACAt                            >  1:2479424/1‑100 (MQ=255)
                                                                         cTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACAt                            >  2:1332910/1‑100 (MQ=255)
                                                                                atatGCAAACTGACACAAAATTGTGTCATAGTGCAGTAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAg                       <  2:2209132/98‑1 (MQ=255)
                                                                                 tatGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCa                    >  1:1545550/1‑100 (MQ=255)
                                                                                  ctgcactatGACACAATTTTGTGTCAGGTTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCAc                   >  1:2610341/9‑100 (MQ=25)
                                                                                   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)
                                                                                     cAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTgg                <  1:717270/100‑1 (MQ=255)
                                                                                       ctatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAAGCACGGGAt              >  2:668374/4‑100 (MQ=255)
                                                                                       aCCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGAt              >  1:762621/1‑100 (MQ=255)
                                                                                       aCCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGAt              <  1:3460788/100‑1 (MQ=255)
                                                                                        ccTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATg             <  2:729895/100‑1 (MQ=255)
                                                                                         cTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGc            >  2:2764533/1‑100 (MQ=255)
                                                                                         atGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGc            >  1:3315286/2‑100 (MQ=255)
                                                                                         atGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGc            >  1:2066097/2‑100 (MQ=255)
                                                                                            acacAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGt         <  1:3473209/100‑1 (MQ=255)
                                                                                            acacAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGt         <  1:2829885/100‑1 (MQ=255)
                                                                                            acacAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGt         <  1:1559450/100‑1 (MQ=255)
                                                                                             cacaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg        >  2:3652688/1‑100 (MQ=255)
                                                                                             cacaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg        <  1:757533/100‑1 (MQ=255)
                                                                                             cacaAAATTGTGTCATAGTGCAGGAAAAAGCACTTACCAGGAGCAGACAACAGCAACGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCg           >  2:3238050/1‑97 (MQ=255)
                                                                                              acaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGa       >  1:781386/1‑100 (MQ=255)
                                                                                                 aaaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATAc    <  1:2204499/100‑1 (MQ=255)
                                                                                                  aaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACg   <  2:1907751/100‑1 (MQ=255)
                                                                                                  aaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACg   <  1:2248015/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGT  >  NC_000913/3059605‑3059803

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: