Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,499,347 +C intergenic (+190/+324) yjgZ → / ← yi41 hypothetical protein/IS4 hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,499,3431.C97.5% 626.7 / 7.4 205intergenic (+186/+328)yjgZ/yi41hypothetical protein/IS4 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (4/1);  new base C (114/83);  total (121/84)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 234 total aligned reads displayed.

ACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGT  >  NC_000913/4499251‑4499430
                                                                                             |                                                                                       
acGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGActagctag                                                                                          >  1:3125488/1‑93 (MQ=255)
  gcgGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGcc                                                                                       >  1:1416319/1‑93 (MQ=255)
    gggCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGcccc                                                                                     >  2:1377618/1‑93 (MQ=255)
     ggCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGccccc                                                                                    >  1:3233986/1‑92 (MQ=255)
      gCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGccccct                                                                                   <  2:200671/93‑3 (MQ=255)
        tGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGccccctac                                                                                 >  2:2780326/1‑89 (MQ=255)
        tGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGccccctac                                                                                 <  1:965148/93‑5 (MQ=255)
          aCAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGc                                                                                <  1:4027152/93‑1 (MQ=255)
            aaTGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCaa                                                                              <  1:3490566/93‑1 (MQ=255)
              tGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAAcc                                                                            >  2:4247482/1‑93 (MQ=255)
                ggTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCaa                                                                          >  1:4299181/1‑93 (MQ=255)
                ggTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCaa                                                                          <  2:4116140/93‑1 (MQ=255)
                 gTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCaaa                                                                         >  2:1962697/1‑93 (MQ=255)
                  tCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCaaat                                                                        >  2:2292623/1‑93 (MQ=255)
                  tCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCaaat                                                                        <  2:2664337/93‑1 (MQ=255)
                  gcAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCaaat                                                                        <  1:639509/92‑1 (MQ=255)
                    aGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCaaataa                                                                      >  2:3287913/1‑93 (MQ=255)
                     gAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCaaataaa                                                                     >  2:58297/1‑93 (MQ=255)
                      aaTATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCaaataaat                                                                    <  2:2629084/93‑1 (MQ=255)
                      aaTATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCaaataaat                                                                    >  2:3440938/1‑93 (MQ=255)
                       atatGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATc                                                                   >  2:1231501/1‑93 (MQ=255)
                       atatGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATc                                                                   >  2:151507/1‑93 (MQ=255)
                       atatGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATc                                                                   <  2:1775218/93‑1 (MQ=255)
                        tatGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTGCCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCa                                                                  >  2:1729749/1‑93 (MQ=255)
                        tatGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCa                                                                  <  1:4481389/93‑1 (MQ=255)
                        tatGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCa                                                                  >  1:3914691/1‑93 (MQ=255)
                        tatGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCa                                                                  >  1:1555902/1‑93 (MQ=255)
                        tatGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCa                                                                  >  2:2447067/1‑93 (MQ=255)
                         atGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCaa                                                                 >  2:2990665/1‑93 (MQ=255)
                         atGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCaa                                                                 <  1:2893641/93‑1 (MQ=255)
                          tGAGTTAGGCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAAc                                                                >  2:4577202/1‑93 (MQ=255)
                          tGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAAc                                                                <  1:3514823/93‑1 (MQ=255)
                          tGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAAc                                                                >  1:2226492/1‑93 (MQ=255)
                            aGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCa                                                              >  1:1873450/1‑93 (MQ=255)
                            aGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCa                                                              >  2:3725487/1‑93 (MQ=255)
                            aGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCa                                                              >  1:3837510/1‑93 (MQ=255)
                            aGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCa                                                              >  2:1050065/1‑93 (MQ=255)
                             gTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCAt                                                             >  1:3259896/1‑93 (MQ=255)
                                 gCCCGGGGTTATTCCTGGGGTCAGTACTTCAGATCTACCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTc                                                         <  2:4369864/93‑1 (MQ=255)
                                 gCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTc                                                         >  1:1248362/1‑93 (MQ=255)
                                 gCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTc                                                         >  2:1013267/1‑93 (MQ=255)
                                 gCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACAATCGTc                                                         >  1:482389/1‑93 (MQ=255)
                                  cccGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCa                                                        <  2:3545613/93‑1 (MQ=255)
                                    cGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCaga                                                      >  2:2092419/1‑93 (MQ=255)
                                    cGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCaga                                                      <  1:1213558/93‑1 (MQ=255)
                                    cGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCaga                                                      <  1:3973900/93‑1 (MQ=255)
                                     ggCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCagag                                                     <  1:2483854/93‑1 (MQ=255)
                                     ggCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCagag                                                     <  2:1951279/93‑1 (MQ=255)
                                        gttgttCCTGGGTTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGc                                                  <  1:371970/93‑1 (MQ=255)
                                        gTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCCACCAAATAAATCAACCATCGTCAGAGTGc                                                  >  2:1594854/1‑93 (MQ=255)
                                        gTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGc                                                  >  2:477840/1‑93 (MQ=255)
                                        gTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGc                                                  >  1:2377908/1‑93 (MQ=255)
                                        gTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGc                                                  >  2:1465715/1‑93 (MQ=255)
                                         ttATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCa                                                 >  2:4308461/1‑93 (MQ=255)
                                          tATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCaa                                                <  1:1634134/93‑1 (MQ=255)
                                           aTTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAAt                                               <  1:3108692/93‑1 (MQ=255)
                                           aTTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAAt                                               >  1:2761383/1‑93 (MQ=255)
                                            ttCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATc                                              >  1:2018334/1‑93 (MQ=255)
                                              ccTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGt                                            >  2:2921779/1‑93 (MQ=255)
                                               cTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGtt                                           <  1:3022696/93‑1 (MQ=255)
                                                tGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTc                                          >  2:2499785/1‑93 (MQ=255)
                                                tGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTc                                          >  1:3062164/1‑93 (MQ=255)
                                                tGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTc                                          >  1:4623453/1‑93 (MQ=255)
                                                tGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTc                                          >  2:935906/1‑93 (MQ=255)
                                                tGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTc                                          <  1:3212818/93‑1 (MQ=255)
                                                  gggTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCat                                        >  1:3800969/1‑93 (MQ=255)
                                                   ggTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCata                                       <  1:2150283/93‑1 (MQ=255)
                                                   ggTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCata                                       >  2:3079466/1‑93 (MQ=255)
                                                     tCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCAtata                                     >  2:3262920/1‑93 (MQ=255)
                                                      taGTACTTCAGAGCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATaca                                   <  1:1460660/92‑1 (MQ=255)
                                                       aGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATaca                                   <  2:3951387/93‑1 (MQ=255)
                                                        gTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATacac                                  >  2:523210/1‑93 (MQ=255)
                                                        gTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATacac                                  <  2:2413073/93‑1 (MQ=255)
                                                        gTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATacac                                  <  2:3973954/93‑1 (MQ=255)
                                                        gTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATacac                                  >  1:3787210/1‑93 (MQ=255)
                                                           cTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGt                               >  1:4667206/1‑93 (MQ=255)
                                                           cTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGt                               <  2:8221/93‑1 (MQ=255)
                                                            ttCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTc                              >  1:341002/1‑93 (MQ=255)
                                                            ttCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTc                              >  2:1492160/1‑93 (MQ=255)
                                                             tCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCg                             >  1:771989/1‑93 (MQ=255)
                                                               aGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGt                           <  1:2510719/93‑1 (MQ=255)
                                                               aGATCTAGCCTCCGGTACCAGACTAGCTAGACCCCTACGCAACCAAATAAATCAACCATAGTCAGAGTGCAATCGTTCATATACACTGTCGGt                           >  1:2317736/1‑93 (MQ=255)
                                                                gATCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTg                          >  2:4664324/1‑93 (MQ=255)
                                                                 aTCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGa                         >  1:2644698/1‑93 (MQ=255)
                                                                 aTCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGa                         <  2:4189420/93‑1 (MQ=255)
                                                                  tCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGAt                        >  2:387334/1‑93 (MQ=255)
                                                                  tCTAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGAt                        >  2:4070156/1‑93 (MQ=255)
                                                                    tAGCCTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATgc                      <  2:3185484/93‑1 (MQ=255)
                                                                    tAGCCTCCGCTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATgc                      <  1:1014859/93‑1 (MQ=255)
                                                                       ccTCCGGTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCt                   >  1:4102646/1‑93 (MQ=255)
                                                                            ggTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTaaaaa              >  2:1126095/1‑93 (MQ=255)
                                                                             gTACCAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTaaaaaa             <  2:4642276/93‑1 (MQ=255)
                                                                              tatcAGAGTAGATAGCCCCCCACGCAACCAAATAAATCAACCATCGTCAGAGTGCACTCGTTCATATACACTGTCGGTGATGCGCTAAAAAAt            <  1:3632926/90‑1 (MQ=255)
                                                                                 cAGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCg         >  1:812978/1‑93 (MQ=255)
                                                                                  aGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGt        <  1:3572989/93‑1 (MQ=255)
                                                                                  aGACTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGt        >  2:2176009/1‑93 (MQ=255)
                                                                                    aCTAGCTAGCCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTa      <  1:46021/93‑1 (MQ=255)
                                                                                      tagctagcCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAAt    <  1:4554161/93‑1 (MQ=255)
                                                                                       agctagcCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATg   <  1:257051/93‑1 (MQ=255)
                                                                                        gctagcCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGt  >  1:1776580/1‑93 (MQ=255)
                                                                                             |                                                                                       
ACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGT  >  NC_000913/4499251‑4499430

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