Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,260,901 C→T 20.9% intergenic (‑98/+27) dauA ← / ← prs C4‑dicarboxylic acid transporter/phosphoribosylpyrophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,260,9010CT20.9% 235.4 / 54.8 119intergenic (‑98/+27)dauA/prsC4‑dicarboxylic acid transporter/phosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (46/48);  new base T (16/9);  total (62/57)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

Only 100 of 132 total aligned reads displayed.

CACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGAC  >  NC_000913/1260801‑1261000
                                                                                                    |                                                                                                   
cacaATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAACccgcc                                                                                                     >  2:3281166/1‑100 (MQ=255)
 acaATACTGTTTCCTCCGGCAGCACATCATCTGTCAACCAAATGGCAGGTGAAGGAGGCATAGCTCATACAAATGGATATTAAAGACAAGAGAACccgcc                                                                                                     <  2:3073845/100‑1 (MQ=255)
 acaATAATGTTTCATCCGTGAGCGCATCATCTATCAACTAAATGGCAGTTGAAGGAGGCATAGGTCATACAAATTGATATTACAGAAAAAAAAACccgcc                                                                                                     <  2:1046712/100‑1 (MQ=255)
         gTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCggg                                                                                             >  1:1463995/1‑100 (MQ=255)
         gTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCggg                                                                                             >  1:749860/1‑100 (MQ=255)
          tttCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGt                                                                                            <  2:2428561/100‑1 (MQ=255)
               tCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCTGCGGCGGGTttttt                                                                                       <  1:3766436/100‑2 (MQ=255)
               tCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCTGCGGCGGGTttttt                                                                                       <  1:23727/100‑2 (MQ=255)
               tCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTg                                                                                       <  1:118551/100‑1 (MQ=255)
               tCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTg                                                                                       <  2:2708907/100‑1 (MQ=255)
                ccGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCTGCGGCGGGTtttttt                                                                                      <  1:3455015/100‑3 (MQ=255)
                ccGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCTGCGGCGGGTtttttt                                                                                      <  1:1437426/100‑3 (MQ=255)
                 cGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCtttttg                                                                                     <  2:3693820/100‑4 (MQ=255)
                 cGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAg                                                                                     >  1:1049760/1‑100 (MQ=255)
                 cGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAg                                                                                     <  2:2558258/100‑1 (MQ=255)
                  gTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCTGCGGCGGGTttttttg                                                                                     <  2:1349527/99‑4 (MQ=255)
                  gTGAGCGCATCATCTGTCAAACAACTGGCAGGTGAAGGAGGCATAGGTCATACAAATGGCTATTACAGACAAATATATCCGCTGCTGCGGGTttttttg                                                                                     <  2:3184890/99‑4 (MQ=255)
                      gcgcATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCggg                                                                                >  1:2088889/1‑100 (MQ=255)
                      gcgcATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCggg                                                                                >  2:199783/1‑100 (MQ=255)
                       cgcATCATCTGTCAACTAAATGGCAGGTGAAGGAGGTATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGt                                                                               <  1:2679573/100‑1 (MQ=255)
                       cgcATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGt                                                                               >  2:1123501/1‑100 (MQ=255)
                        gcAACAACTGTCAACTAAGTGGCAGGTGAAGGAGGCATAGCTCATACAAATGTATATTACAGACAAAAAGACCCGCTGCGGCGGGTttttttgtctgta                                                                               <  2:2545005/99‑10 (MQ=21)
                         catcatCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTc                                                                             >  1:1563825/1‑100 (MQ=255)
                         catcatCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTc                                                                             >  1:2831081/1‑100 (MQ=255)
                          atcatcTGCCAAGTAAATGGCAGGTGAAGGAGTCATAGGTCATACAAATGGATATTACAGACAAAAATACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCg                                                                            <  1:3628342/100‑1 (MQ=255)
                            catcTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGAt                                                                          <  2:409631/100‑1 (MQ=255)
                              tcTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTa                                                                        <  2:519072/100‑1 (MQ=255)
                                tGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCGGGGTTCGATTAgt                                                                      <  2:2363861/100‑1 (MQ=255)
                                tGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgt                                                                      >  2:1835015/1‑100 (MQ=255)
                                tGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgt                                                                      >  2:677766/1‑100 (MQ=255)
                                 gTCAACTAAATGGCAGGTGAAGGAGGCATAGTTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgtg                                                                     <  1:865529/100‑1 (MQ=255)
                                 gTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgtg                                                                     <  1:1042980/100‑1 (MQ=255)
                                   cAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGtt                                                                   >  1:978463/1‑100 (MQ=255)
                                   cAACTAAATGGCAGGTGAAGGAGGCATAGGGCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGtt                                                                   <  2:2011096/100‑1 (MQ=255)
                                      cTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGa                                                                <  1:3292591/100‑1 (MQ=255)
                                      cTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGa                                                                >  2:2531426/1‑100 (MQ=255)
                                      cTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAACGGATATTACAGACAACACAAACCGCCGCAACGGGTCTTGGAACCGGGTTCGATTAGTGTTCGa                                                                >  2:3117748/1‑100 (MQ=255)
                                        aaaTGGCAGGTGAAGGAGGCATATGTCAGACAATTGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGa                                                                <  2:3574055/98‑1 (MQ=255)
                                           tGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATg                                                           >  1:3082298/1‑100 (MQ=255)
                                               aGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCag                                                       >  1:213211/1‑100 (MQ=255)
                                               aGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCa                                                        >  2:2168831/1‑99 (MQ=255)
                                                ggTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCaga                                                      >  1:2131048/1‑100 (MQ=255)
                                                 gTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCagag                                                     <  1:375406/100‑1 (MQ=255)
                                                     aggaggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCg                                                 >  1:95094/1‑100 (MQ=255)
                                                       gaggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAt                                               <  1:1912384/100‑1 (MQ=255)
                                                           cATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttctt                                           >  1:2197968/1‑100 (MQ=255)
                                                           cATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttctt                                           >  2:2791212/1‑100 (MQ=255)
                                                             tAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCg                                         <  1:3080575/100‑1 (MQ=255)
                                                              aGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGt                                        <  1:199783/100‑1 (MQ=255)
                                                               ggTCATACAAATGGATATTACAGACAAAAACACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGCGTTCGAACATGGCAGAGATCGATTCTTCg                                         >  2:1864341/1‑98 (MQ=255)
                                                               ggTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGtt                                       <  1:2652158/100‑1 (MQ=255)
                                                               ggTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGtt                                       >  2:1538615/1‑100 (MQ=255)
                                                                gccATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTg                                      <  2:3721352/98‑1 (MQ=255)
                                                                gTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTg                                      <  2:2378007/100‑1 (MQ=255)
                                                                 tCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGc                                     <  1:1752086/100‑1 (MQ=255)
                                                                 tCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGc                                     <  1:3417484/100‑1 (MQ=255)
                                                                  cATACAAATGGATATTACATACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGc                                     <  2:486542/99‑1 (MQ=255)
                                                                  cATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCt                                    <  1:564597/100‑1 (MQ=255)
                                                                   aTACAAATGGTTATTACAGACAAAAAAGCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGATATCGATTCTTCGTTGCTg                                   <  2:109856/100‑1 (MQ=255)
                                                                      cAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATa                                >  1:1602925/1‑100 (MQ=255)
                                                                       aaaTGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATac                               <  1:2492762/100‑1 (MQ=255)
                                                                        aaTGGATATTACAGACAAAAAACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacg                              <  1:2179081/100‑1 (MQ=255)
                                                                        aaTGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGATAGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacg                              >  2:3397763/1‑100 (MQ=255)
                                                                          tGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacgac                            >  1:2127251/1‑100 (MQ=255)
                                                                          tGGATATAACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacgac                            <  2:2980364/100‑1 (MQ=255)
                                                                               atTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATc                       <  2:1380104/100‑1 (MQ=255)
                                                                                 tACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGc                     >  1:1109414/1‑100 (MQ=255)
                                                                                   cAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCtt                   >  2:449500/1‑100 (MQ=255)
                                                                                   cAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCtt                   >  1:2749904/1‑100 (MQ=255)
                                                                                    aGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTc                  >  1:3543294/1‑100 (MQ=255)
                                                                                     gACAAAAAAGCCCGCCGCAGCGGGTATTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTc                  <  2:143547/99‑1 (MQ=255)
                                                                                     gACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCg                 <  1:3384316/100‑1 (MQ=255)
                                                                                     gACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTCCGTTGCTGATACGACGAATCGCTTCg                 <  1:2315730/100‑1 (MQ=255)
                                                                                      gctcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGTACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgg                >  2:1289804/9‑100 (MQ=38)
                                                                                      gctcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGATATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgg                >  2:3687738/9‑100 (MQ=37)
                                                                                      gctcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgg                >  1:2975324/9‑100 (MQ=38)
                                                                                       ctcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGc               >  1:2239468/8‑100 (MQ=38)
                                                                                       ctcaaagaCCCGCTGCGGCGGGTCTTTGAGCCGGGTTCGATCAGCCTTCGAACATGCCAGCATTCGATTCTTCGTTGCTGCTACGACGAATCGCTTCgg                >  2:1733787/4‑99 (MQ=255)
                                                                                       cAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGc               >  1:1754010/1‑100 (MQ=255)
                                                                                       cAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGc               <  1:3715689/100‑1 (MQ=255)
                                                                                        tcaaagaCCCGCTGCGGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGcc              >  1:1961836/7‑100 (MQ=255)
                                                                                        aaaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGcc              <  1:2925527/100‑1 (MQ=255)
                                                                                        aaaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGcc              >  1:935827/1‑100 (MQ=255)
                                                                                         naaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGcc              >  2:1545068/2‑99 (MQ=255)
                                                                                         caaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGcc              >  2:81520/2‑99 (MQ=39)
                                                                                         caaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCa             >  2:1541849/2‑100 (MQ=255)
                                                                                         caaagaCCCGCTGCGGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCa             >  1:3337507/2‑100 (MQ=255)
                                                                                         aaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCa             >  2:3794472/1‑100 (MQ=255)
                                                                                          aaagaCCCGCTTCGGCGGGTCTTTGAGCCGGGTTCGGTTAGTGTTCGAACATGGCATGGATCGATTGTTCGTTGCTGATACGACGAATCGCTTCGGcc              >  2:766276/1‑98 (MQ=255)
                                                                                          aaagaCCCGCTGCGGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAg            >  1:3203610/1‑100 (MQ=255)
                                                                                          aaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGACTCGCTTCGGCCAg            >  2:1487334/1‑100 (MQ=255)
                                                                                            aaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCa          <  2:260398/100‑1 (MQ=255)
                                                                                               cccGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATAc       >  1:754342/1‑100 (MQ=255)
                                                                                                ccgcggCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCTCTTCGGCCAGCAt         <  2:1303032/97‑1 (MQ=255)
                                                                                                ccgccgCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATAcc      <  1:1620441/100‑1 (MQ=255)
                                                                                                ccgccgCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATAcc      <  1:2492376/100‑1 (MQ=255)
                                                                                                  gccgccGCGCGTCTTTGAGCCGGGTTCTAGTAGTGTTCGAACATGGCAGAGTTCTATTCTTCGTTGCTCATACGACGAATCGCTTCGGCCAGCATACCTg    <  2:864170/100‑1 (MQ=255)
                                                                                                  gccgcAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTg    >  1:146265/1‑100 (MQ=255)
                                                                                                    cgcAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGAc  >  1:3314753/1‑100 (MQ=255)
                                                                                                    cgcAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTATTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGAc  >  1:2321425/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGAC  >  NC_000913/1260801‑1261000

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: