Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,149,302 T→C Q459R (CAG→CGG)  ppc ← phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,149,3020TC93.4% 661.2 / 24.4 212Q459R (CAG→CGG) ppcphosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base T (8/6);  new base C (86/112);  total (94/118)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

AGCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGG  >  NC_000913/4149212‑4149391
                                                                                          |                                                                                         
agcagGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCgg                                                                                         >  1:2711100/1‑93 (MQ=255)
agcagGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCgg                                                                                         <  2:2669034/93‑1 (MQ=255)
 gcagGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGc                                                                                        <  1:2934913/93‑1 (MQ=255)
  cagGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCa                                                                                       <  2:2308664/93‑1 (MQ=255)
   agGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAg                                                                                      >  2:4154554/1‑93 (MQ=255)
   agGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAg                                                                                      <  2:363519/93‑1 (MQ=255)
    ggTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAgg                                                                                     >  1:3248279/1‑93 (MQ=255)
    ggTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAgg                                                                                     >  2:5739559/1‑93 (MQ=255)
     gTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGt                                                                                    <  2:4162512/93‑1 (MQ=255)
     gTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGt                                                                                    >  1:1781747/1‑93 (MQ=255)
      tGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTa                                                                                   <  1:5236565/93‑1 (MQ=255)
      tGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTa                                                                                   >  1:5797542/1‑93 (MQ=255)
        gACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATc                                                                                 <  1:4964535/93‑1 (MQ=255)
        gACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATc                                                                                 <  2:246838/93‑1 (MQ=255)
           aGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAg                                                                              <  1:1037793/93‑1 (MQ=255)
           aGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAg                                                                              >  1:1496759/1‑93 (MQ=255)
           aGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGACTTGCGGTGCTTCGGCAATCACCCGGCAGGTATGGAg                                                                              >  1:2356231/1‑93 (MQ=255)
            gCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGc                                                                             >  2:2713142/1‑93 (MQ=255)
            gCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGc                                                                             <  1:3317874/93‑1 (MQ=255)
            gCCAGTACGTCGGACGGCGTTTTCGCCATCGAAATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGc                                                                             <  1:5966870/93‑1 (MQ=255)
             ccAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCa                                                                            <  1:625120/93‑1 (MQ=255)
             ccAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCa                                                                            <  1:4248974/93‑1 (MQ=255)
              cAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCAc                                                                           <  2:1145931/93‑1 (MQ=255)
                gTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACtt                                                                         <  1:1608883/93‑1 (MQ=255)
                gTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACtt                                                                         <  2:2641272/93‑1 (MQ=255)
                 tACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTc                                                                        <  2:3298094/93‑1 (MQ=255)
                 tACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTc                                                                        >  1:4829318/1‑93 (MQ=255)
                    gTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTcgcg                                                                     >  2:3917501/1‑93 (MQ=255)
                    gTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGGATCGAGCACTTcgcg                                                                     >  2:2996814/1‑93 (MQ=255)
                     tCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTcgcgc                                                                    >  1:5897780/1‑93 (MQ=255)
                         gcGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGttt                                                                <  2:5801487/92‑1 (MQ=255)
                             cGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGgc                                                            >  2:31591/1‑93 (MQ=255)
                             cGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGgc                                                            <  2:5299473/93‑1 (MQ=255)
                              gTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGgcg                                                           <  2:776548/93‑1 (MQ=255)
                               ttttCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGTGTTTCGgcgc                                                          >  1:868184/1‑93 (MQ=255)
                               ttttCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGgcgc                                                          <  1:5424890/93‑1 (MQ=255)
                               ttttCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGgcgc                                                          <  2:2588015/93‑1 (MQ=255)
                               ttttCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGgcgc                                                          <  1:4163364/93‑1 (MQ=255)
                               ttttCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGgcgc                                                          <  1:1449939/93‑1 (MQ=255)
                                   cGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTgg                                                      <  2:1476009/93‑1 (MQ=255)
                                   cGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTgg                                                      >  1:3265720/1‑93 (MQ=255)
                                    gCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGt                                                     >  1:4549399/1‑93 (MQ=255)
                                    gCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGt                                                     <  2:2393900/93‑1 (MQ=255)
                                     ccATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGtt                                                    <  1:4524994/93‑1 (MQ=255)
                                     ccATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGtt                                                    <  2:5448643/93‑1 (MQ=255)
                                       aTCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGc                                                  >  2:6061557/1‑93 (MQ=255)
                                        tCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGcc                                                 >  1:2468374/1‑93 (MQ=255)
                                        tCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGcc                                                 >  2:5587808/1‑93 (MQ=255)
                                          gagaTCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAg                                               >  2:4496738/1‑93 (MQ=255)
                                          gagaTCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAg                                               <  1:4356494/93‑1 (MQ=255)
                                           agaTCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATTACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTCGTTGCCAGt                                              >  1:3888609/1‑93 (MQ=255)
                                            gaTCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCCGtt                                             >  2:1464710/1‑93 (MQ=255)
                                            gaTCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGtt                                             <  2:3067748/93‑1 (MQ=255)
                                             aTCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTg                                            >  1:2052163/1‑93 (MQ=255)
                                             aTCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTAGCGCGTTTCGGCGCTTGGTTGCCAGTTg                                            >  1:21956/1‑93 (MQ=255)
                                              tCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTgc                                           >  1:3107161/1‑93 (MQ=255)
                                              tCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTgc                                           >  1:138735/1‑93 (MQ=255)
                                               cACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTgcg                                          <  2:5228130/93‑1 (MQ=255)
                                                aCGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTgcgc                                         >  1:3048036/1‑93 (MQ=255)
                                                 cGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGGATCGAGCACTTCGCGCGTTTCGGCGCGTGGTTGCCAGTTGcgcg                                        >  2:3033165/1‑93 (MQ=255)
                                                  gTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGTCAGTTGCGCgg                                       >  1:2826195/1‑93 (MQ=255)
                                                   tAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGc                                      >  2:2742889/1‑93 (MQ=255)
                                                    aGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCa                                     >  2:3651228/1‑93 (MQ=255)
                                                       gggcAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGaa                                  <  2:4756640/92‑1 (MQ=255)
                                                       cggcAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGaa                                  <  2:1860964/93‑1 (MQ=255)
                                                        ggcAATGGAGCCTTGCGGTGCTTCGGCAATCGCCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAg                                 <  1:1863774/93‑1 (MQ=255)
                                                         gcAATGGAGCCTTGCGGTGCTTCGGCAATCCCCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGc                                <  2:212350/93‑1 (MQ=255)
                                                         gcAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGc                                <  2:4696854/93‑1 (MQ=255)
                                                         gcAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGc                                >  1:3970542/1‑93 (MQ=255)
                                                         gcAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGc                                <  1:3729462/93‑1 (MQ=255)
                                                          cAATGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCg                               <  2:3370981/93‑1 (MQ=255)
                                                           aaTGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCgg                              <  2:6028164/93‑1 (MQ=255)
                                                            aTGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGa                             <  2:3749345/93‑1 (MQ=255)
                                                             tGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGAc                            >  2:1962859/1‑93 (MQ=255)
                                                             tGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGAc                            <  2:1594043/93‑1 (MQ=255)
                                                             tGGAGCCTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCACTTGGTTGCCAGTTGCGCGGCAGAAGCGGAc                            >  1:5693217/1‑93 (MQ=255)
                                                                   cTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTgg                      >  1:3226556/1‑93 (MQ=255)
                                                                   cTTGCGGTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTgg                      <  2:2851067/93‑1 (MQ=255)
                                                                         gTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttca                <  1:5050511/93‑1 (MQ=255)
                                                                         gTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttca                <  2:3463649/93‑1 (MQ=255)
                                                                         gTGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGaggtca                >  1:2110887/1‑93 (MQ=255)
                                                                          tGCTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttcag               <  1:1074191/93‑1 (MQ=255)
                                                                          tGCTTCGGCAATCACCCGGCAGGCATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttcag               <  1:3483129/93‑1 (MQ=255)
                                                                            cTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttcagtt             <  2:4977125/93‑1 (MQ=255)
                                                                            cTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttcagtt             <  2:2196965/93‑1 (MQ=255)
                                                                            cTTCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttcagtt             >  2:4904567/1‑93 (MQ=255)
                                                                             ttCGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGagttcagttc            <  1:745290/93‑1 (MQ=255)
                                                                               cGGCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTcgc          <  1:1242828/93‑1 (MQ=255)
                                                                                ggCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTcgcg         >  2:851401/1‑93 (MQ=255)
                                                                                ggCAATCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTcgcg         <  2:3881955/93‑1 (MQ=255)
                                                                                   aaTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGAt      <  2:3923460/93‑1 (MQ=255)
                                                                                   aaTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGAt      >  2:1247798/1‑93 (MQ=255)
                                                                                   aaTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGAt      >  1:99048/1‑93 (MQ=255)
                                                                                   aaTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGAt      >  1:5111456/1‑93 (MQ=255)
                                                                                   aaTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGAt      <  1:1939374/93‑1 (MQ=255)
                                                                                    aTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATc     <  2:2800849/93‑1 (MQ=255)
                                                                                    aTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATc     <  2:1217086/93‑1 (MQ=255)
                                                                                    aTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATc     >  1:4021549/1‑93 (MQ=255)
                                                                                    aTCACCCGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATc     <  1:1667407/93‑1 (MQ=255)
                                                                                     tCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGTGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCa    <  2:2159592/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGG  >  NC_000913/4149212‑4149391

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: