Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,782,932 T→C Y16H (TAT→CAT)  yibO → putative 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,782,9320TC99.5% 796.3 / ‑4.8 217Y16H (TAT→CAT) yibOputative 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (118/98);  total (119/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 250 total aligned reads displayed.

AGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACC  >  NC_000913/3782842‑3783021
                                                                                          |                                                                                         
agtaAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGatggccat                                                                                         >  2:5546503/1‑93 (MQ=255)
 gtaAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGatggccatg                                                                                        >  2:1803598/1‑93 (MQ=255)
 gtaAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGatggccatg                                                                                        <  2:4240425/93‑1 (MQ=255)
 gtaAAATTACACAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGatggccatg                                                                                        >  1:3811292/1‑93 (MQ=255)
    aaaTTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGatggccatggct                                                                                     <  2:2559630/93‑1 (MQ=255)
     aaTTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGAtggccatggcta                                                                                    <  1:5024337/93‑1 (MQ=255)
     aaTTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGAtggccatggcta                                                                                    <  1:1001967/93‑1 (MQ=255)
      aTTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATggccatggctat                                                                                   <  2:3945770/93‑1 (MQ=255)
      aTTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATggccatggctat                                                                                   <  1:244265/93‑1 (MQ=255)
      aTTACCGAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATggccatggctat                                                                                   <  2:1928602/93‑1 (MQ=255)
       ttACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATc                                                                                  >  1:4010991/1‑93 (MQ=255)
       ttACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATc                                                                                  >  2:3267476/1‑93 (MQ=255)
          cGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATcgcg                                                                               >  1:4594986/1‑93 (MQ=255)
          cGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATcgcg                                                                               <  1:3750699/93‑1 (MQ=255)
          cGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATcgcg                                                                               <  1:1120503/93‑1 (MQ=255)
            cAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCgaa                                                                             >  1:1324133/1‑93 (MQ=255)
            cAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCgaa                                                                             <  1:5040649/93‑1 (MQ=255)
               aTTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCgaagaa                                                                          <  1:5142578/93‑1 (MQ=255)
                ttttGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAAc                                                                         <  2:1533841/93‑1 (MQ=255)
                ttttGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAAc                                                                         <  1:2903215/93‑1 (MQ=255)
                ttttGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAAc                                                                         >  1:999621/1‑93 (MQ=255)
                  ttGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAAcag                                                                       <  2:4613807/93‑1 (MQ=255)
                    gACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAAccgca                                                                     <  1:6020230/93‑1 (MQ=255)
                    gACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAAcagca                                                                     <  1:2280479/93‑1 (MQ=255)
                    gACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAAcagca                                                                     <  2:1491245/93‑1 (MQ=255)
                     aCTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAAcagcag                                                                    <  2:321469/93‑1 (MQ=255)
                     aCTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAAcagcag                                                                    >  1:1186746/1‑93 (MQ=255)
                      ctctTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAgg                                                                   >  1:474752/1‑93 (MQ=255)
                          tGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATaa                                                               <  1:319950/93‑1 (MQ=255)
                          tGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATaa                                                               <  2:2078060/93‑1 (MQ=255)
                          tGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATaa                                                               <  1:295766/93‑1 (MQ=255)
                          tGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATaa                                                               <  2:2641104/93‑1 (MQ=255)
                           gAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAAc                                                              >  2:4055152/1‑93 (MQ=255)
                           gAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAAc                                                              >  1:2789226/1‑93 (MQ=255)
                             gTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGc                                                            >  2:436652/1‑93 (MQ=255)
                               aTGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCa                                                          >  2:6131410/1‑93 (MQ=255)
                               aTGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCa                                                          <  2:3112722/93‑1 (MQ=255)
                                tGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCAt                                                         >  1:721591/1‑93 (MQ=255)
                                   ggTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCAtttt                                                      <  2:5925680/93‑1 (MQ=255)
                                    gTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCAttttt                                                     <  2:5454327/93‑1 (MQ=255)
                                     ttGTCGCAATGTTGGTTTCTAAAAAGCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTa                                                    <  1:1465452/93‑1 (MQ=255)
                                     ttGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTa                                                    >  1:21433/1‑93 (MQ=255)
                                      tGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAg                                                   <  1:1247792/93‑1 (MQ=255)
                                      tGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAg                                                   <  2:5631544/93‑1 (MQ=255)
                                        tCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTg                                                 <  2:3546076/93‑1 (MQ=255)
                                          gCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCt                                               >  2:1066377/1‑93 (MQ=255)
                                            aaTGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTaa                                             >  1:3136532/1‑93 (MQ=255)
                                            aaTGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTaa                                             <  1:2991131/93‑1 (MQ=255)
                                              tGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTaaaa                                           <  1:4042254/93‑1 (MQ=255)
                                              tGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAAAAACAGGAGAACGCCATTTTTTTTTCTaaac                                           >  2:1271903/1‑92 (MQ=255)
                                               gTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAAc                                          >  1:4940587/1‑93 (MQ=255)
                                                ttGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAAcc                                         >  1:3244518/1‑93 (MQ=255)
                                                  ggTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAAcccc                                       >  2:4564750/1‑93 (MQ=255)
                                                   gTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCg                                      >  2:3351483/1‑93 (MQ=255)
                                                   gTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCg                                      >  1:258409/1‑93 (MQ=255)
                                                   gTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCg                                      >  2:3689909/1‑93 (MQ=255)
                                                   gTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCg                                      <  2:2788369/93‑1 (MQ=255)
                                                    tttCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCgg                                     >  1:1485959/1‑93 (MQ=255)
                                                     ttCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGt                                    <  2:3624136/93‑1 (MQ=255)
                                                       cTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTaa                                  <  2:5956518/93‑1 (MQ=255)
                                                       cTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTaa                                  <  1:4430672/93‑1 (MQ=255)
                                                       cTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTaa                                  <  2:3359930/93‑1 (MQ=255)
                                                       cTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTaa                                  <  1:4250904/93‑1 (MQ=255)
                                                        tAAAAAACCTATGGTACTGGGGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAAt                                 >  2:2382874/1‑93 (MQ=255)
                                                         aaaaaaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATg                                >  2:5854510/1‑93 (MQ=255)
                                                          aaaaaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATgg                               <  1:5705939/93‑1 (MQ=255)
                                                          aaaaaCCTATGGGCCTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATgg                               >  1:1501074/1‑93 (MQ=255)
                                                           aaaaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGa                              <  1:745132/93‑1 (MQ=255)
                                                           aaaaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGa                              <  1:1979101/93‑1 (MQ=255)
                                                            aaaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGAt                             >  2:4695027/1‑93 (MQ=255)
                                                            aaaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGAt                             <  1:1373683/93‑1 (MQ=255)
                                                             aaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATg                            >  1:5387355/1‑93 (MQ=255)
                                                               ccTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCa                          >  2:3837226/1‑93 (MQ=255)
                                                                cTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCAc                         >  1:144410/1‑93 (MQ=255)
                                                                cTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCAc                         <  2:5685865/93‑1 (MQ=255)
                                                                 tATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACt                        >  1:3859066/1‑93 (MQ=255)
                                                                  aTGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACtg                       >  1:1275286/1‑93 (MQ=255)
                                                                   tGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACtgt                      >  1:6059452/1‑93 (MQ=255)
                                                                     gTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTgg                    <  1:3626020/93‑1 (MQ=255)
                                                                      tACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTggg                   <  2:3030139/93‑1 (MQ=255)
                                                                      tACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTggg                   >  2:2008994/1‑93 (MQ=255)
                                                                        cTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGcc                 >  1:3519041/1‑93 (MQ=255)
                                                                        cTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGcc                 >  1:896621/1‑93 (MQ=255)
                                                                         tGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGGTGCACTGTGGGCCa                >  1:4144788/1‑93 (MQ=255)
                                                                         tGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCa                >  1:3032718/1‑93 (MQ=255)
                                                                          ggTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCaa               <  2:5274678/93‑1 (MQ=255)
                                                                           gTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAAt              >  1:2716227/1‑93 (MQ=255)
                                                                           gTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAAt              >  1:4646338/1‑93 (MQ=255)
                                                                             gATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCg            >  1:3970671/1‑93 (MQ=255)
                                                                             gATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCg            >  1:5557811/1‑93 (MQ=255)
                                                                               ttCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTc          >  1:4527396/1‑93 (MQ=255)
                                                                               ttCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTc          <  2:4741378/93‑1 (MQ=255)
                                                                                tCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTcc         <  1:2293999/93‑1 (MQ=255)
                                                                                tCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGAAGCACTGTGGGCCAATCGTcc         <  1:5866578/93‑1 (MQ=255)
                                                                                 cTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCg        <  1:2914464/93‑1 (MQ=255)
                                                                                 cTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCg        >  2:5151782/1‑93 (MQ=255)
                                                                                  tGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGc       >  1:679200/1‑93 (MQ=255)
                                                                                   ggATGGCCATGGCTATCGCGAAGAACAGCAGGAAAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCa      >  1:5237317/1‑93 (MQ=255)
                                                                                    gATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCAt     >  1:106572/1‑93 (MQ=255)
                                                                                      tggccatggctaTCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATAc   >  1:613807/1‑93 (MQ=255)
                                                                                          |                                                                                         
AGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACC  >  NC_000913/3782842‑3783021

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: