Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,957,124 +G intergenic (+58/+26) 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,1231.G98.9% 561.0 / ‑0.8 183intergenic (+57/+27)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (50/131);  total (52/131)
Fisher's exact test for biased strand distribution p-value = 7.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 204 total aligned reads displayed.

GACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTC  >  NC_000913/3957033‑3957210
                                                                                                 |                                                                                       
gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                          <  1:1013524/93‑1 (MQ=21)
gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                          <  1:1783949/93‑1 (MQ=21)
 aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                                                        <  2:622719/93‑1 (MQ=21)
  cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                                                       >  2:3191891/1‑92 (MQ=21)
  cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                                                       <  1:2237195/93‑2 (MQ=21)
  cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCACG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                                                       <  2:4658949/93‑2 (MQ=21)
     atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcct                                                                                    <  2:1523584/93‑5 (MQ=21)
         gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                                                                 <  2:1923959/93‑1 (MQ=18)
          aagcGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                                                                <  1:3253288/90‑1 (MQ=18)
          aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                                                                <  2:2000249/93‑1 (MQ=18)
          aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctca                                                                                <  1:2069487/93‑2 (MQ=18)
            aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                              <  1:2586138/93‑1 (MQ=18)
              gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                            <  1:4455735/93‑1 (MQ=18)
              gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                            <  1:3093559/93‑1 (MQ=18)
                 ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCC‑‑ACCCCCTCCCCCCGGGGGGGGGGGGAAAAACCCAAAGCCTTCTCCTCGAGaa                                                                         >  1:1952279/1‑93 (MQ=12)
                      gTTGCGGTTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACCCTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                    <  1:3432653/93‑1 (MQ=17)
                      gTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                    <  1:4627217/93‑1 (MQ=18)
                      gTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCC‑‑ACCCTCTCCCCCCGGGGGAGGGGGGAAACAACCCAAGCCTTCTCCTCGAGAAGGCCt                                                                    >  2:6379163/1‑93 (MQ=16)
                        tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                                                                  <  2:2910210/93‑1 (MQ=18)
                          cGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                                                <  2:2130620/93‑1 (MQ=18)
                          cGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                                                <  1:1302114/93‑1 (MQ=18)
                          cGGGTTAAGTGCGCGCTGATGCCCTCACCCCC‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                                                >  2:1003920/1‑93 (MQ=18)
                            ggTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                                                              <  2:3128559/93‑1 (MQ=18)
                               tAAGTGCGCGCTGATGCCCTCCCCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                                                           <  2:5460211/93‑1 (MQ=18)
                                aaGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                                                          <  1:2658925/93‑1 (MQ=18)
                                 aGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                                                         <  1:2281406/93‑1 (MQ=18)
                                 aGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                                                         <  1:5691133/93‑1 (MQ=18)
                                  gTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                                                        <  1:3979240/93‑1 (MQ=18)
                                      gcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt                                                    <  2:3741859/93‑1 (MQ=21)
                                       cgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTa                                                   <  2:764208/93‑1 (MQ=21)
                                        gcTGATGCCGTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTAc                                                  <  1:1044885/93‑1 (MQ=21)
                                        gcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTAc                                                  <  2:83492/93‑1 (MQ=21)
                                          tGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg                                                <  1:908507/93‑1 (MQ=25)
                                          cgATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg                                                <  2:3119394/92‑1 (MQ=25)
                                           gATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                                               <  1:555552/93‑1 (MQ=25)
                                            aTGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCa                                              <  1:1230155/93‑1 (MQ=25)
                                            aTGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCa                                              <  2:6324766/93‑1 (MQ=25)
                                             tGCCCTCACCCCGCCC‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAc                                             <  1:3453432/93‑1 (MQ=25)
                                              gCCCTCACCCCG‑‑ACCCTCTACCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                                            >  2:5520604/1‑93 (MQ=25)
                                               cccTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt                                           <  1:3684870/93‑1 (MQ=25)
                                                ccTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt                                          <  2:5280271/93‑1 (MQ=25)
                                                ccTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt                                          <  2:5104223/93‑1 (MQ=25)
                                                  tCACCCCG‑‑ACGCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTGGCGGTACAGCACCTTaa                                        <  2:3182267/93‑1 (MQ=37)
                                                  tCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAAGCCTTCTCCTGGAGAAGGGCTTGCTATTAGTTGCGGTACAGCACCTTaa                                        >  1:6123546/1‑93 (MQ=25)
                                                  tCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACAATCAAGGCCTTCTCCTGGAGAAAGCCTTGCTATTAGTTTCGGTACAGCACCTTaa                                        >  2:5642319/1‑93 (MQ=25)
                                                   cACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                       <  2:2853495/93‑1 (MQ=37)
                                                           gaccctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgtg                                 <  1:1824227/89‑1 (MQ=38)
                                                             ccctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat                               >  2:702516/3‑93 (MQ=255)
                                                              cctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGata                              <  1:6223647/92‑1 (MQ=255)
                                                              cctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGata                              >  1:2582342/2‑93 (MQ=255)
                                                               ctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGatat                             <  1:6244084/93‑1 (MQ=255)
                                                               ctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGatat                             >  1:2905703/1‑93 (MQ=255)
                                                               ctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGatat                             >  1:5625701/1‑93 (MQ=255)
                                                                tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATc                            <  1:2769854/93‑1 (MQ=255)
                                                                  tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCg                          >  2:2142715/1‑93 (MQ=255)
                                                                    ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGaa                        <  2:324246/93‑1 (MQ=255)
                                                                     cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                       >  2:5312705/1‑93 (MQ=255)
                                                                     cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                       <  1:5541453/93‑1 (MQ=255)
                                                                      acaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACt                      <  2:379790/93‑1 (MQ=255)
                                                                      acaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACt                      >  1:431313/1‑93 (MQ=255)
                                                                       caGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTg                     <  1:693811/93‑1 (MQ=255)
                                                                        agggagagggagAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                    >  2:4572181/1‑93 (MQ=255)
                                                                        agggagagggagAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                    <  2:4154943/93‑1 (MQ=255)
                                                                        agggagagggagAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                    <  1:4982155/93‑1 (MQ=255)
                                                                         gggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTggg                   <  1:1339427/93‑1 (MQ=255)
                                                                         gggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTggg                   <  2:2372566/93‑1 (MQ=255)
                                                                         gggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTggg                   >  1:5802113/1‑93 (MQ=255)
                                                                         gggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTggg                   <  1:3164676/93‑1 (MQ=255)
                                                                          ggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt                  <  2:4895567/93‑1 (MQ=255)
                                                                          ggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt                  <  1:334185/93‑1 (MQ=255)
                                                                           gagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtg                 <  1:1438588/93‑1 (MQ=255)
                                                                             gagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGtgtg               >  1:3111704/1‑93 (MQ=255)
                                                                             gagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGtgtg               >  2:438591/1‑93 (MQ=255)
                                                                             gagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGtgtg               >  1:5119969/1‑93 (MQ=255)
                                                                              agggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGc              <  2:1765105/93‑1 (MQ=255)
                                                                               gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATCTCCGAACTGGGTGTGCa             >  2:2699521/1‑93 (MQ=255)
                                                                               gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCa             <  2:6278330/93‑1 (MQ=255)
                                                                                ggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAg            <  2:3673147/93‑1 (MQ=255)
                                                                                ggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAg            <  2:5712953/93‑1 (MQ=255)
                                                                                ggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAg            >  2:3928783/1‑93 (MQ=255)
                                                                                  agaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCg          >  2:4471239/1‑93 (MQ=255)
                                                                                   gaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg         >  2:117353/1‑93 (MQ=255)
                                                                                   gaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg         <  1:4195062/93‑1 (MQ=255)
                                                                                   gaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg         <  1:1494056/93‑1 (MQ=255)
                                                                                    aaaaCACTCAAGGCCTTTTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCggg        >  2:2621518/1‑93 (MQ=255)
                                                                                    aaaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCggg        <  1:4876554/93‑1 (MQ=255)
                                                                                    aaaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCggg        >  1:6086207/1‑93 (MQ=255)
                                                                                     aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc       <  2:4734772/93‑1 (MQ=255)
                                                                                     aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc       >  2:2216629/1‑93 (MQ=255)
                                                                                     aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc       <  2:6080620/93‑1 (MQ=255)
                                                                                     aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc       >  2:3865531/1‑93 (MQ=255)
                                                                                      aaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGcc      <  2:1269961/93‑1 (MQ=255)
                                                                                      aaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGcc      <  1:5420713/93‑1 (MQ=255)
                                                                                       acacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACTTTAATGATGTGATATCCGAACTGGGTGTGCAGCGCGCCg     >  2:4859304/1‑93 (MQ=255)
                                                                                       acacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCg     >  1:3887205/1‑93 (MQ=255)
                                                                                          acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGt   <  2:591680/93‑1 (MQ=255)
                                                                                          acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGt   >  1:5057743/1‑93 (MQ=255)
                                                                                          acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGt   <  2:3607545/93‑1 (MQ=255)
                                                                                          acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGAAGCGGGCCGGt   >  1:3801826/1‑93 (MQ=255)
                                                                                           cTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTc  <  2:675546/93‑1 (MQ=255)
                                                                                                 |                                                                                       
GACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTC  >  NC_000913/3957033‑3957210

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