Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,957,092 +CC intergenic (+26/+58) ilvC → / ← ppiC ketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,957,0911.C94.4% 218.0 / 11.3 108intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base C (54/48);  total (58/50)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
*NC_0009133,957,0912.C94.4% 215.8 / 11.3 108intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base C (54/47);  total (58/50)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

AACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTT  >  NC_000913/3957017‑3957155
                                                                             ||                                                                   
aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCACTCACCCCG‑‑A‑‑c                                                                    >  1:5579258/1‑76 (MQ=255)
 aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑cc                                                                   <  1:5063904/76‑2 (MQ=255)
 aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑cc                                                                   <  1:4062646/76‑2 (MQ=255)
 aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCC‑‑A‑‑cc                                                                   >  2:2937085/1‑75 (MQ=255)
  cTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑ccc                                                                  >  1:1655960/1‑74 (MQ=255)
  cTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑ccc                                                                  <  2:1947871/76‑3 (MQ=255)
   tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑cccc                                                                 >  2:1019727/1‑73 (MQ=255)
    gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccctc                                                                >  1:3375651/1‑76 (MQ=255)
    gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccctc                                                                >  1:2443665/1‑76 (MQ=255)
      gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccctctc                                                              >  2:5006736/1‑74 (MQ=255)
       cgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCC‑‑A‑‑Ccctctcc                                                             >  2:2294856/1‑73 (MQ=38)
          cTATATGACAGATATGATTCGCATTGCTGTTGCGGGTTAAGTGCGCGCTCATGCCCTCGCCCCG‑‑ACCCTCTCCcac                                                            <  1:3104484/76‑1 (MQ=255)
           tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCcaca                                                           >  2:5185836/1‑76 (MQ=25)
           tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCcaca                                                           >  1:3407468/1‑76 (MQ=25)
             tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagg                                                         <  1:5150382/76‑1 (MQ=25)
             tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagg                                                         >  1:1756682/1‑76 (MQ=25)
              atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCATaggg                                                        <  2:2911017/76‑1 (MQ=25)
               tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACaggga                                                       >  1:4032077/1‑76 (MQ=25)
               tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCG‑‑ACCCTCTCCCACagggg                                                       >  1:396461/1‑75 (MQ=25)
                 aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggaga                                                     <  1:568118/76‑1 (MQ=21)
                  cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACcgggagag                                                    >  1:2660181/1‑76 (MQ=21)
                  cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagag                                                    >  2:360313/1‑76 (MQ=21)
                  cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagag                                                    <  1:2373068/76‑1 (MQ=21)
                   aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑TCCCTCTCCCACagggagagg                                                   <  2:5078886/76‑1 (MQ=18)
                   aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagagg                                                   <  1:170736/76‑1 (MQ=21)
                    gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagaggg                                                  >  1:838738/1‑76 (MQ=21)
                    gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagaggg                                                  >  1:540589/1‑76 (MQ=21)
                    gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagaggg                                                  >  1:2320191/1‑76 (MQ=21)
                       atGAAAGGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAgggagagggaga                                               <  2:154338/76‑1 (MQ=18)
                        tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGaa                                              <  1:2230881/76‑1 (MQ=18)
                         gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGaaa                                             <  2:1252669/76‑1 (MQ=18)
                         gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGaaa                                             >  2:2744983/1‑76 (MQ=18)
                          aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACCGGGAGAGGGAGaaaa                                            >  1:3572536/1‑76 (MQ=18)
                          aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGTGAGGGAGaaaa                                            >  2:2085213/1‑76 (MQ=18)
                           aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAac                                           <  1:5241498/76‑1 (MQ=18)
                           aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAac                                           <  1:4353170/76‑1 (MQ=18)
                            aCGTATAGCTGTTGCGGGTTAAGTCCGCGCTGATGCCCTCCCCCCGCCC‑‑CTCTCCCACAGGGAGAGGGAGAAAaca                                          <  1:751280/76‑1 (MQ=14)
                             cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAacac                                         <  2:2637887/76‑1 (MQ=17)
                              gTATTGCTGTTGCGGGTTAAATGCGCGCCGATGCCCTCACCCCG‑‑ACCCTCTCACACAGGGAGAGGGAGAAAacaac                                        >  2:181878/1‑74 (MQ=17)
                               tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑CCCCTCTCCCACAGGGAGAGGGAGAAACCACTc                                       <  2:5416311/76‑1 (MQ=14)
                               tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACCGGGCGAGGGCGAAAACACTc                                       >  2:360567/1‑76 (MQ=17)
                                 ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagg                                                         >  1:4672413/1‑56 (MQ=17)
                                 ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                   >  2:2183171/1‑76 (MQ=17)
                                 ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                   >  1:4098605/1‑76 (MQ=17)
                                   gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGTGGGAGAAAACACTCAAgg                                >  1:667915/1‑75 (MQ=17)
                                   gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                >  2:3331788/1‑75 (MQ=17)
                                   gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                >  1:2672247/1‑75 (MQ=17)
                                    cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggc                               >  1:2894758/1‑74 (MQ=17)
                                     tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcc                              >  2:3383084/1‑73 (MQ=17)
                                     tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcc                              >  2:4053990/1‑73 (MQ=17)
                                     tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcc                              <  1:2029390/76‑4 (MQ=17)
                                     tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcc                              <  2:2672247/76‑4 (MQ=17)
                                        tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                            <  2:5180557/76‑1 (MQ=14)
                                         gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                           <  1:1389653/76‑1 (MQ=14)
                                         gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                           <  1:3918899/76‑1 (MQ=14)
                                         gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                           >  1:3556923/1‑76 (MQ=14)
                                          cGGGTTAAGTGCGCGCTGATGGCCCCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGGAAACACACAAAGCCGtctc                          >  1:3648481/1‑76 (MQ=2)
                                           gggTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCCCCGGGGGCGGGAGGAAACACTCAAAGCCCTCTcc                         >  1:2448633/1‑76 (MQ=9)
                                           gggTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACaggga                                                       >  1:1013914/1‑48 (MQ=11)
                                           gggTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                         >  1:3834939/1‑76 (MQ=14)
                                            ggTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                        <  2:5838529/76‑1 (MQ=14)
                                            ggTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                        <  1:494626/76‑1 (MQ=14)
                                               tAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgta                     <  1:2826526/76‑3 (MQ=14)
                                               tAAGGGCGCGCTGATGCCCTCACCCCC‑‑ACCCTCTCCCACAGGGAGAGGGGGACAACACTCACAGCGTTCTCCTggc                     >  2:63938/1‑75 (MQ=2)
                                                aattGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgag                    <  2:2108225/73‑1 (MQ=14)
                                                aaGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGaag                    >  1:5059748/1‑76 (MQ=14)
                                                 aGTGCGCGCTGATGCCCTCACCCCG‑‑ACGCTGTCGCGCGGGGGGAGGGAGAAAACACTCAGGGCCTTCTCCTGgaga                   <  2:453276/76‑1 (MQ=9)
                                                 aGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgaga                   <  2:184316/76‑1 (MQ=14)
                                                  gTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                  <  2:5185756/76‑1 (MQ=14)
                                                  gTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                  <  1:740381/76‑1 (MQ=14)
                                                   tGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgagacg                 >  2:1598999/1‑74 (MQ=14)
                                                   tGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                 <  2:4408198/76‑1 (MQ=17)
                                                    gcgcgcTGATGCCCTCACCCGT‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                <  1:2184182/76‑1 (MQ=17)
                                                     cgcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc               <  1:5745061/76‑1 (MQ=17)
                                                      gcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc              >  1:1500671/1‑76 (MQ=17)
                                                      gcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc              >  1:3271886/1‑76 (MQ=17)
                                                      gcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc              >  1:419565/1‑76 (MQ=17)
                                                      gcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc              >  1:4704333/1‑76 (MQ=17)
                                                       cgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt             >  2:1008446/1‑76 (MQ=17)
                                                       cgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt             >  2:4567663/1‑76 (MQ=17)
                                                       cgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt             <  2:4455165/76‑1 (MQ=17)
                                                        gcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCCTCTCCTGGCGAAGGccct            >  2:2552597/1‑74 (MQ=14)
                                                        gcTGATGCCCTCACCCAG‑‑ACCCTCTCCCACACGGAGAGGGAGAAAACACACAAGGCCTTCTCCGCCAGAAGGCCta            >  2:2485822/1‑75 (MQ=2)
                                                         cTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg           <  2:5258299/76‑1 (MQ=18)
                                                         cTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg           >  2:3072015/1‑76 (MQ=18)
                                                          tGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc          <  2:4587282/76‑1 (MQ=18)
                                                           gATGCCCTCACCCCGGCG‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTTGCt         <  2:5579256/76‑1 (MQ=17)
                                                           gATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt         >  2:4392837/1‑76 (MQ=18)
                                                           gATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt         <  1:3334456/76‑1 (MQ=18)
                                                           gATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt         >  2:4500325/1‑76 (MQ=18)
                                                            aTGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCtn        <  1:2937085/76‑2 (MQ=18)
                                                               cccTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGCAGGCCTTGCTATTa     >  2:5566711/1‑76 (MQ=18)
                                                               cccTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa     >  1:753237/1‑76 (MQ=21)
                                                               cccTCACCCCG‑‑ACCCTCGCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa     <  2:5179180/76‑1 (MQ=21)
                                                                 cTCACCCCG‑‑ACCCTCTCCCACAGGGAGCGGGAGAAAACACTCAACGCCTTCTCCTGGAGCAGGCCTTGCTACTAGt   >  2:2920327/1‑76 (MQ=18)
                                                                 cTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt   >  1:4638063/1‑76 (MQ=21)
                                                                 cTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt   >  2:5453732/1‑76 (MQ=21)
                                                                 cTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt   <  1:1019727/76‑1 (MQ=21)
                                                                  tCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt  >  1:4734886/1‑76 (MQ=21)
                                                                  tCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt  <  2:3556923/76‑1 (MQ=21)
                                                                             ||                                                                   
AACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTT  >  NC_000913/3957017‑3957155

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