Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,152,760 G→T 100% N113K (AAC→AAA ppc ← phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,152,7600GT97.2% 337.5 / 1.4 106N113K (AAC→AAAppcphosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (53/50);  total (55/51)
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.82e-01

Only 100 of 107 total aligned reads displayed.

CGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCC  >  NC_000913/4152665‑4152856
                                                                                               |                                                                                                
cgacgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTtttt                                                                                              >  1:3032445/1‑100 (MQ=255)
cgacgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTtttt                                                                                              >  2:2214957/1‑100 (MQ=255)
cgacgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCAACTGCTTTTTTGCTGGTGTCTTCGCTCAGTTCCGGCTGttt                                                                                                >  2:2321123/1‑98 (MQ=255)
  acgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTc                                                                                            <  1:2923807/100‑1 (MQ=255)
  acgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTc                                                                                            <  2:189469/100‑1 (MQ=255)
   cgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCa                                                                                           >  1:2797707/1‑100 (MQ=255)
   cgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCa                                                                                           >  2:1372502/1‑100 (MQ=255)
    gggTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAg                                                                                          <  1:2833627/100‑1 (MQ=255)
     ggTAATTTCGGTTGGGTGAGCCGTGATGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGt                                                                                         <  1:393653/100‑1 (MQ=255)
       tAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGttt                                                                                       <  1:810/100‑1 (MQ=255)
       tAATGTGGGTTGGCTGAGCCGTGAGGTTTAGTTCCAGCGACAGCGATTCCTTGGCTTTTTTGCTGTTGTCTTCGTTTAGTTCCGACCTTTTTTTCAGtat                                                                                       <  2:500853/100‑3 (MQ=255)
         aTTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTAc                                                                                     >  1:973723/1‑100 (MQ=255)
          tttCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACg                                                                                    >  1:1693842/1‑100 (MQ=255)
          tttCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACg                                                                                    >  1:2263618/1‑100 (MQ=255)
           ttCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGc                                                                                   <  2:196088/100‑1 (MQ=255)
            tctgTTGGGTGAGCCGTGAGGACCAGTTCCAGTGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCa                                                                                  <  1:1029547/97‑1 (MQ=255)
             cGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCa                                                                                  <  2:2900514/99‑1 (MQ=255)
              ggTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAgg                                                                                >  1:1834249/1‑100 (MQ=255)
              ggTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTAATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCATTTTACGCAg                                                                                 >  1:781251/1‑99 (MQ=255)
               gTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAggg                                                                               <  2:2038270/100‑1 (MQ=255)
                   ggTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGc                                                                            >  2:1547032/1‑99 (MQ=255)
                    gTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCgg                                                                          >  1:575834/1‑100 (MQ=255)
                    gTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCgg                                                                          <  1:2360668/100‑1 (MQ=255)
                     tGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCggg                                                                         <  1:290261/100‑1 (MQ=255)
                      gAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCggtc                                                                        <  1:1220061/100‑3 (MQ=255)
                       aGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCg                                                                       >  1:52536/1‑100 (MQ=255)
                          cGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATc                                                                    <  2:499509/100‑1 (MQ=255)
                           gTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATc                                                                    <  1:260624/99‑1 (MQ=255)
                           gTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCa                                                                   <  1:751421/100‑1 (MQ=255)
                               ggACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTc                                                               >  1:1706073/1‑100 (MQ=255)
                                gACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGCTGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGAGGGCGATCACTTcc                                                              >  1:359046/1‑100 (MQ=255)
                                gACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTcc                                                              >  2:2981718/1‑100 (MQ=255)
                                  ccAGTTCCAGCGACAGCGATTCCACTGCTTTTTTTATTTTTTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCgg                                                            >  2:800087/1‑100 (MQ=255)
                                  ccAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCgg                                                            >  1:2209673/1‑100 (MQ=255)
                                  ccAGTTCCAGCGACAGCAATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCgg                                                            <  1:543347/100‑1 (MQ=255)
                                   cAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCggg                                                           <  1:1192161/100‑1 (MQ=255)
                                   cAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCggg                                                           <  2:1391187/100‑1 (MQ=255)
                                   cAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCggg                                                           >  2:606331/1‑100 (MQ=255)
                                   cAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCggg                                                           <  1:1386732/100‑1 (MQ=255)
                                    acttCCATCGCCCGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTTCGGCTGTTTTTGCAGTGTAAGCAGGGTTTGGGCGACCACTTc                                                               >  2:1399982/3‑95 (MQ=255)
                                     gTTCCAGCGACAGCGATTCCACTGCTTTTTTTATTGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGtt                                                         >  2:1621936/1‑100 (MQ=255)
                                     gTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGAGTGCGGGCGATCACTTCCGGGtt                                                         >  1:2542129/1‑100 (MQ=255)
                                     gTTCCAGCGACAGCAATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGtt                                                         >  2:1020704/1‑100 (MQ=255)
                                        ccAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCt                                                      <  2:290635/100‑1 (MQ=255)
                                            cgacagcgaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCa                                                  <  1:2552193/100‑1 (MQ=255)
                                              acagcgaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGc                                                <  2:2599921/100‑1 (MQ=255)
                                               cagcgaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCt                                               >  1:1954049/1‑100 (MQ=255)
                                                agcgaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCtt                                              >  1:631550/1‑100 (MQ=255)
                                                    aTTCCACTGCTTTTTTGATGGTGTCTTCGCTTAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGcc                                          <  1:109409/100‑1 (MQ=255)
                                                     ttCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCt                                         <  2:2021961/100‑1 (MQ=255)
                                                       ccACTGCTTTTTTTTTTTTTTTTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTTGCAGCTTCGCCttt                                       >  2:2065750/1‑100 (MQ=255)
                                                        cACTGCTTTTTTTTTTTTTTTTTCACTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGTTTGCTGGCAGCTTCTCCttt                                       >  2:1108211/1‑99 (MQ=255)
                                                        cACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCttt                                       <  2:54519/99‑1 (MQ=255)
                                                        cACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTc                                      >  2:2880546/1‑100 (MQ=255)
                                                        cACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTc                                      >  1:13346/1‑100 (MQ=255)
                                                         aCTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTTTTTTTTTTTATGGTACGGGCGATCACTTGCTGTTTTCTGTCAGCTTCACCttt                                       >  2:1469055/1‑98 (MQ=255)
                                                         aCTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCg                                     <  1:1021398/100‑1 (MQ=255)
                                                          cTGCTTTTTTTATTGTTTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCg                                     >  2:1554288/1‑99 (MQ=255)
                                                           tGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCggn                                   <  1:468865/100‑2 (MQ=255)
                                                           tGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGc                                   <  1:1979020/100‑1 (MQ=255)
                                                           tGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGc                                   <  1:592659/100‑1 (MQ=255)
                                                           tGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGc                                   <  1:2261745/100‑1 (MQ=255)
                                                            gCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCg                                  >  1:2075942/1‑100 (MQ=255)
                                                            gCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCg                                  >  1:218944/1‑100 (MQ=255)
                                                            gCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCg                                  >  1:2544360/1‑100 (MQ=255)
                                                            gCTTTTTTGATGGTGACTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCg                                  <  2:1372388/100‑1 (MQ=255)
                                                             cTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTTTTTTTTTTTATAGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCg                                  >  2:2149607/1‑99 (MQ=255)
                                                             cTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGa                                 >  2:2166329/1‑100 (MQ=255)
                                                              ttttttGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGaa                                >  1:549185/1‑100 (MQ=255)
                                                               tttttGATGGTGTCTTCGCTAAGTTCCAGCTGTTTTTTCAGTTTACGAAGGGTGCGGGCGATCACTTCAGGGTTGCTGGCAGCTTCGCCTTTCAGCGaaa                               >  2:1490041/1‑100 (MQ=255)
                                                                 tttGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTTTTTTTTTTTAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGGCTTTCGGCGGAATg                             >  2:1541927/1‑100 (MQ=255)
                                                                    gATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCtg                          >  2:2792168/1‑100 (MQ=255)
                                                                    gATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCtg                          <  1:710325/100‑1 (MQ=255)
                                                                    gATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCGATAGGCGaagctgac                          <  1:2053894/100‑7 (MQ=255)
                                                                    gATGGTGTCTTCGATCAGTTCCGACTGTTTTTACAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTAGCAGCTTCGGCTTTCGGCGACATGc                            >  2:1490933/1‑98 (MQ=255)
                                                                     aTGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCtgt                         <  1:216747/100‑1 (MQ=255)
                                                                     aTGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCtgt                         <  2:13346/100‑1 (MQ=255)
                                                                      tGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCtgt                         <  2:3032412/99‑1 (MQ=255)
                                                                      tGGTGTCTTCGCTCAGTTCCGGCTGTTTTTACAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCAAAATGCtgtg                        >  2:1993890/1‑100 (MQ=255)
                                                                       ggTGTCTTCGCTCAGTTCCTGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTgg                       <  1:2497596/100‑1 (MQ=255)
                                                                       ggTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTgg                       <  1:772434/100‑1 (MQ=255)
                                                                       ggTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTgg                       <  1:582647/100‑1 (MQ=255)
                                                                        gtgtCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGt                      <  1:1851456/100‑1 (MQ=255)
                                                                         tgtCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTa                     >  1:182653/1‑100 (MQ=255)
                                                                          gtCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTAt                    <  2:333619/100‑1 (MQ=255)
                                                                           tCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTAtt                   >  1:1359504/1‑100 (MQ=255)
                                                                            cTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGTTCTTg                  >  2:1058866/1‑100 (MQ=255)
                                                                            cTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTg                  >  1:917058/1‑100 (MQ=255)
                                                                            cTTCGCTCAGTTCCGGCTGGTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTg                  >  2:2263852/1‑100 (MQ=255)
                                                                             ttCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGc                 <  1:281623/100‑1 (MQ=255)
                                                                                gCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTcg              >  1:174716/1‑100 (MQ=255)
                                                                                 cTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGTGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTcgg             <  2:2808004/100‑1 (MQ=255)
                                                                                 cTCAGTTCCGGCTGTTTTTTATGTTTTCGCAAGGAGCGGGCGATCACTTCCGGATTGCTAGCAGCTTCGACTTTCGGCGAAATGCTGTGatatatct                >  2:1855045/1‑93 (MQ=255)
                                                                                   cAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTcggc            >  1:576574/1‑99 (MQ=255)
                                                                                     gTTCCGGCTGTTTTTTCAGTTTACGCCGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTAGGCGAAATGCTGTGGTATTGCTCGGCGgt         >  1:895066/1‑100 (MQ=255)
                                                                                          ggCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTg     >  2:2248834/1‑99 (MQ=255)
                                                                                           gCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTgg    <  1:2620708/99‑1 (MQ=255)
                                                                                            cTGTTTTTTCTGTTTATTCAGAGAACGAGCGACCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGTCGCAATTCTTTTGTATTACTCGGCGGttttaacc  >  2:2996179/1‑96 (MQ=255)
                                                                                            cTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGcc  >  1:1097824/1‑100 (MQ=255)
                                                                                            cTGTTTGTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGc   <  2:1624745/99‑1 (MQ=255)
                                                                                               |                                                                                                
CGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCC  >  NC_000913/4152665‑4152856

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