Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,476,406 A→G D124G (GAC→GGC)  yjgD → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,476,4060AG99.6% 795.9 / ‑4.6 231D124G (GAC→GGC) yjgDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (101/129);  total (102/129)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

Only 100 of 247 total aligned reads displayed.

TTGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGT  >  NC_000913/4476316‑4476495
                                                                                          |                                                                                         
ttGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGgcggcg                                                                                         <  1:5434677/93‑1 (MQ=255)
 tGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGcggcga                                                                                        <  2:5395811/93‑1 (MQ=255)
 tGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGcggcga                                                                                        <  2:3874258/93‑1 (MQ=255)
   aacaacTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATg                                                                                      >  1:2870176/1‑93 (MQ=255)
    acaacTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGa                                                                                     >  2:3168468/1‑93 (MQ=255)
     caacTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGaa                                                                                    >  2:4317715/1‑93 (MQ=255)
      aacTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGGTGAAg                                                                                   >  2:6124155/1‑93 (MQ=255)
       acTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGTGGCGATGAAGa                                                                                  <  1:3647463/93‑1 (MQ=255)
       acTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGa                                                                                  >  1:3708463/1‑93 (MQ=255)
       acTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGa                                                                                  >  1:6170533/1‑93 (MQ=255)
          gatgaCGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGAttt                                                                               >  2:82743/1‑93 (MQ=255)
           atgaCGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGAtttt                                                                              >  1:6088489/1‑93 (MQ=255)
             gaCGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGt                                                                            >  1:2777559/1‑93 (MQ=255)
              aCGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTc                                                                           <  1:500847/93‑1 (MQ=255)
              aCGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTc                                                                           <  1:6015662/93‑1 (MQ=255)
              aCGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTc                                                                           >  1:1894259/1‑93 (MQ=255)
              aCGCTGACAGAGAAATTTGACGTTGAGTACGACGGATTGGGCACTTTCTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTc                                                                           <  2:952615/93‑1 (MQ=255)
               cGCTGGCAAAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCg                                                                          >  1:6293814/1‑93 (MQ=255)
                gCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGa                                                                         <  2:2428002/93‑1 (MQ=255)
                gCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGa                                                                         <  2:5019342/93‑1 (MQ=255)
                 cTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGACGATCCTAACGGCGAAGATGGCGGGGATGAACATTTTGTCggg                                                                        >  2:4015580/1‑91 (MQ=255)
                 cTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGAt                                                                        <  2:4079813/93‑1 (MQ=255)
                 cTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGAt                                                                        <  1:6007011/93‑1 (MQ=255)
                     cAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAg                                                                    >  2:5507034/1‑93 (MQ=255)
                      agagAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGa                                                                   >  2:2723369/1‑93 (MQ=255)
                      agagAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGa                                                                   >  1:5559045/1‑93 (MQ=255)
                      agagAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGa                                                                   <  1:1142216/93‑1 (MQ=255)
                       gagaAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGAc                                                                  >  2:5581338/1‑93 (MQ=255)
                        agaAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACg                                                                 >  2:5932748/1‑93 (MQ=255)
                         gaAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGa                                                                >  2:3657028/1‑93 (MQ=255)
                         gaAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGa                                                                <  1:4911029/93‑1 (MQ=255)
                          aaaTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGAt                                                               <  1:4727189/93‑1 (MQ=255)
                          aaaTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGAt                                                               <  1:560318/93‑1 (MQ=255)
                           aaTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATg                                                              <  2:1897165/93‑1 (MQ=255)
                           aaTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCAAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATg                                                              <  1:4532411/93‑1 (MQ=255)
                            aTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGa                                                             <  1:2210236/93‑1 (MQ=255)
                            aTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGa                                                             <  1:3033234/93‑1 (MQ=255)
                            aTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGa                                                             <  1:660315/93‑1 (MQ=255)
                             tttGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGAc                                                            <  2:3800325/93‑1 (MQ=255)
                               tGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACgg                                                          <  1:4263155/93‑1 (MQ=255)
                               tGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACgg                                                          >  2:3748562/1‑93 (MQ=255)
                                gACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGa                                                         >  2:3097698/1‑93 (MQ=255)
                                gACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGa                                                         >  2:5720988/1‑93 (MQ=255)
                                   gTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGtt                                                      >  2:3177559/1‑93 (MQ=255)
                                   gTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGtt                                                      <  2:1636798/93‑1 (MQ=255)
                                    ttGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTc                                                     <  1:3588859/93‑1 (MQ=255)
                                     tGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCg                                                    <  2:5247661/93‑1 (MQ=255)
                                         tACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCAc                                                <  2:2219879/93‑1 (MQ=255)
                                           cgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTa                                              <  2:2785126/93‑1 (MQ=255)
                                            gacgGATGGGCCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTaa                                             <  1:5489065/93‑1 (MQ=255)
                                                gATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTaa                                         >  2:5799676/1‑93 (MQ=255)
                                                gATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTaa                                         <  2:473886/93‑1 (MQ=255)
                                                gATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTaa                                         >  1:3802380/1‑93 (MQ=255)
                                                 aTGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAg                                        <  1:4821822/93‑1 (MQ=255)
                                                 aTGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAg                                        <  1:1249823/93‑1 (MQ=255)
                                                 aTGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAg                                        <  1:411411/93‑1 (MQ=255)
                                                  tGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGc                                       <  1:1743979/93‑1 (MQ=255)
                                                   ggggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCg                                      >  1:2661140/1‑93 (MQ=255)
                                                    gggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGt                                     <  2:6439133/93‑1 (MQ=255)
                                                      gCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAg                                   <  1:1693474/93‑1 (MQ=255)
                                                       cACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAgc                                  <  2:1259531/93‑1 (MQ=255)
                                                       cACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAgc                                  <  1:6126948/93‑1 (MQ=255)
                                                        acttacttTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAgcg                                 >  1:4497440/1‑93 (MQ=255)
                                                         cttacttTGAAGATCCTAACGGCGAAGATGGCGGCGGTGAAGATTTTGTCGATGAAGACGACGACGGAGATCGCCACTAAATAAGCGTAgcgc                                >  1:3283976/1‑93 (MQ=255)
                                                         cttacttTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAgcgc                                >  1:5863811/1‑93 (MQ=255)
                                                            acttTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATc                             <  2:301895/93‑1 (MQ=255)
                                                            acttTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATc                             <  1:2940080/93‑1 (MQ=255)
                                                             cttTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCa                            <  1:888402/93‑1 (MQ=255)
                                                              tttGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAg                           >  2:4620682/1‑93 (MQ=255)
                                                              tttGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAg                           <  2:2106290/93‑1 (MQ=255)
                                                               ttGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAgg                          <  2:3201726/93‑1 (MQ=255)
                                                                  aaGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGccc                       <  2:5902716/93‑1 (MQ=255)
                                                                  aaGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGccc                       <  1:1953220/93‑1 (MQ=255)
                                                                    gATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAg                     >  2:5155101/1‑93 (MQ=255)
                                                                    gATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAg                     >  1:3186712/1‑93 (MQ=255)
                                                                    gATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAg                     <  2:1669566/93‑1 (MQ=255)
                                                                     aTCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAgg                    >  2:4779314/1‑93 (MQ=255)
                                                                     aTCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAgg                    >  2:6007192/1‑93 (MQ=255)
                                                                     aTCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAgg                    >  1:1096919/1‑93 (MQ=255)
                                                                     aTCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAgg                    >  2:2333887/1‑93 (MQ=255)
                                                                      tCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGt                   >  2:5752185/1‑93 (MQ=255)
                                                                      tCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGt                   >  1:5872060/1‑93 (MQ=255)
                                                                       ccTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTa                  >  2:2688587/1‑93 (MQ=255)
                                                                       ccTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTa                  <  1:642382/93‑1 (MQ=255)
                                                                        cTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAg                 >  2:454109/1‑93 (MQ=255)
                                                                           aCGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTc              >  2:2825816/1‑93 (MQ=255)
                                                                           aCGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTc              >  2:6478432/1‑93 (MQ=255)
                                                                           aCGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTc              <  2:122867/93‑1 (MQ=255)
                                                                           aCGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTc              >  2:5095004/1‑93 (MQ=255)
                                                                            cGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTcc             <  1:2102746/93‑1 (MQ=255)
                                                                            cGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTcc             >  1:5892948/1‑93 (MQ=255)
                                                                            cGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTcc             >  1:136574/1‑93 (MQ=255)
                                                                              gCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCGGGCCCAGGTAGTTCCTg           >  2:630179/1‑93 (MQ=255)
                                                                              gCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTg           >  2:1036073/1‑93 (MQ=255)
                                                                                  aGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAAc       >  1:3982678/1‑93 (MQ=255)
                                                                                    aTGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCt     <  1:3441799/93‑1 (MQ=255)
                                                                                     tGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTg    >  2:109638/1‑93 (MQ=255)
                                                                                     tGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTg    <  2:1779785/93‑1 (MQ=255)
                                                                                     tGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTg    <  2:4300785/93‑1 (MQ=255)
                                                                                       gcggcgATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGt  <  1:1422586/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGT  >  NC_000913/4476316‑4476495

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