Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,067,548 T→C H137R (CAC→CGC)  b3879 ← putative aldose‑1‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,067,5480TC99.6% 753.8 / ‑4.4 234H137R (CAC→CGC) b3879putative aldose‑1‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (116/117);  minor base A (0/1);  total (116/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

Only 100 of 266 total aligned reads displayed.

TTGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCTGGTG  >  NC_000913/4067458‑4067637
                                                                                          |                                                                                         
 tGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCcgtgcgtg                                                                                        <  1:5812906/93‑1 (MQ=255)
 tGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCcgtgcgtg                                                                                        <  2:1626795/93‑1 (MQ=255)
 tGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACACTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCcgtgcgtg                                                                                        <  1:4085148/93‑1 (MQ=255)
  gCAGGTAAAATAAGAATGAATCCCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCgtgcgtgc                                                                                       <  2:5422244/93‑1 (MQ=255)
  gCAGATTAAATAAAGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCgtgcgtgc                                                                                       <  1:4469907/93‑1 (MQ=255)
   cAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGtgcgtgcg                                                                                      >  2:3815597/1‑93 (MQ=255)
      gTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTg                                                                                   >  2:4225456/1‑93 (MQ=255)
      gTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTg                                                                                   >  1:457417/1‑93 (MQ=255)
       ttAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGGCAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGa                                                                                  <  1:2871499/93‑1 (MQ=255)
        tAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAg                                                                                 >  1:5477117/1‑93 (MQ=255)
         aaaTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGc                                                                                >  2:4859302/1‑93 (MQ=255)
            tAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGaa                                                                             >  2:1155220/1‑93 (MQ=255)
            tAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGaa                                                                             >  1:4639026/1‑93 (MQ=255)
             aaGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAAt                                                                            <  2:4500918/93‑1 (MQ=255)
               ggATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATaa                                                                          >  1:6094176/1‑93 (MQ=255)
               ggATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATaa                                                                          >  1:658877/1‑93 (MQ=255)
               ggATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATaa                                                                          >  2:4080442/1‑93 (MQ=255)
                gATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATaaa                                                                         >  1:3278857/1‑93 (MQ=255)
                  tGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAaaca                                                                       <  1:4256861/93‑1 (MQ=255)
                   gAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAaacaa                                                                      <  2:2634201/93‑1 (MQ=255)
                   gAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAaacaa                                                                      >  1:3423194/1‑93 (MQ=255)
                    aaTACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAaacaac                                                                     >  1:4945864/1‑93 (MQ=255)
                     aTACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACt                                                                    <  1:3464621/93‑1 (MQ=255)
                     aTACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACt                                                                    <  2:1527558/93‑1 (MQ=255)
                       aCCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtg                                                                  <  2:3497133/93‑1 (MQ=255)
                       aCCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtg                                                                  >  1:1847675/1‑93 (MQ=255)
                        cccACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtgt                                                                 <  1:256508/93‑1 (MQ=255)
                        cccACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtgt                                                                 <  2:4013210/93‑1 (MQ=255)
                        cccACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGCGCGTGCGTTGAGCGAATAAACAACTtgt                                                                 >  1:1366726/1‑93 (MQ=255)
                          cACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGa                                                               >  1:1665742/1‑93 (MQ=255)
                           aCGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGaa                                                              <  2:2499190/93‑1 (MQ=255)
                           aCGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGaa                                                              <  1:869245/93‑1 (MQ=255)
                            cGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAg                                                             <  2:4571891/93‑1 (MQ=255)
                             gCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGc                                                            <  1:1545275/93‑1 (MQ=255)
                              ccATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGcc                                                           >  1:5030657/1‑93 (MQ=255)
                               cATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCa                                                          >  1:4207638/1‑93 (MQ=255)
                               cATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCa                                                          >  1:3802843/1‑93 (MQ=255)
                               cATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCa                                                          <  2:275907/93‑1 (MQ=255)
                               cATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCa                                                          <  1:5205979/93‑1 (MQ=255)
                                  aaGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCa                                                       <  2:222548/93‑1 (MQ=255)
                                     gcgcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaa                                                    <  1:1234580/93‑1 (MQ=255)
                                     gcgcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaa                                                    <  2:5255394/93‑1 (MQ=255)
                                      cgcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaaa                                                   <  2:4575892/93‑1 (MQ=255)
                                      cgcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaaa                                                   <  1:5998462/93‑1 (MQ=255)
                                        cTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTTCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACg                                                 >  1:3612496/1‑93 (MQ=255)
                                         tGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACgg                                                >  2:883042/1‑93 (MQ=255)
                                         tGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACgg                                                >  1:1753535/1‑93 (MQ=255)
                                         tGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACgg                                                <  1:5421907/93‑1 (MQ=255)
                                          gCGACAGTGCAGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGTGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGa                                               <  1:4335478/93‑1 (MQ=255)
                                           cGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGAt                                              <  1:5576947/93‑1 (MQ=255)
                                           cGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGAt                                              <  1:4707761/93‑1 (MQ=255)
                                           cGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGAt                                              <  1:2312599/93‑1 (MQ=255)
                                            gACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATa                                             <  2:440216/93‑1 (MQ=255)
                                            gACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATa                                             <  2:2958132/93‑1 (MQ=255)
                                             aCAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAg                                            <  2:3367199/93‑1 (MQ=255)
                                             aCAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAg                                            >  2:6292501/1‑93 (MQ=255)
                                             aCAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAg                                            >  2:202272/1‑93 (MQ=255)
                                               aGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGcc                                          >  1:2043404/1‑93 (MQ=255)
                                                gTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCg                                         >  1:1969668/1‑93 (MQ=255)
                                                 tGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGt                                        >  1:5871884/1‑93 (MQ=255)
                                                 tGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGt                                        <  2:1572530/93‑1 (MQ=255)
                                                  gCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACCGATAGCCgca                                       >  2:5224180/1‑91 (MQ=255)
                                                    cGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAAc                                     <  2:5619220/93‑1 (MQ=255)
                                                       tatTCTGGCTTTCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTgg                                  <  2:1871789/93‑1 (MQ=255)
                                                       tatTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTgg                                  <  2:1746537/93‑1 (MQ=255)
                                                       tatTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTgg                                  >  2:6308247/1‑93 (MQ=255)
                                                       tatTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTgg                                  <  2:1199234/93‑1 (MQ=255)
                                                        atTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTggg                                 >  1:130846/1‑93 (MQ=255)
                                                         ttCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGggc                                >  1:4861171/1‑93 (MQ=255)
                                                            tGGCTTGCAATCTCAACGGTAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGggcggc                             >  2:6119471/1‑93 (MQ=255)
                                                              gCTTGCGATCTCAACGGAAAGACCCGTGCGGGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCCGCAg                           >  2:3290476/1‑93 (MQ=255)
                                                              gCTTGCGATCTCAACGGAAAGACCCGTGAGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAg                           <  2:5884714/93‑1 (MQ=255)
                                                                ttGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGCATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGaa                         >  1:6226903/1‑93 (MQ=255)
                                                                 tGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGaaa                        >  2:5673500/1‑93 (MQ=255)
                                                                  gCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGaaaa                       <  1:3091302/93‑1 (MQ=255)
                                                                  gCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGaaaa                       >  2:1347244/1‑93 (MQ=255)
                                                                    gATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAGCGGATAGCCGTAACTGGGCGGCAGAAAAGc                     <  1:1010432/93‑1 (MQ=255)
                                                                     aTCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCt                    <  1:2649182/93‑1 (MQ=255)
                                                                     aTCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCt                    <  2:5099459/93‑1 (MQ=255)
                                                                       ctcAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGt                  <  2:779614/93‑1 (MQ=255)
                                                                        tcAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTc                 >  1:2467869/1‑93 (MQ=255)
                                                                        tcAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTc                 <  1:4481468/93‑1 (MQ=255)
                                                                        tcAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTc                 >  2:5524861/1‑93 (MQ=255)
                                                                           aCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAgg              >  1:276585/1‑93 (MQ=255)
                                                                           aCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAgg              >  1:1113317/1‑93 (MQ=255)
                                                                            cGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAggg             >  2:180880/1‑93 (MQ=255)
                                                                             ggAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAggga            >  1:1210632/1‑92 (MQ=255)
                                                                             ggAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGt            >  1:377191/1‑93 (MQ=255)
                                                                              gAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTa           >  1:250897/1‑93 (MQ=255)
                                                                              gAAAGACCCTGGCTGCCTTTAGCCAATAAACCACTTGGGGAACCCGGAAAAACCGAATGGCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTaa          <  2:4158641/93‑1 (MQ=255)
                                                                               aaaGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTaa          <  2:598641/93‑1 (MQ=255)
                                                                               aaaGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTaa          <  1:5755259/93‑1 (MQ=255)
                                                                                aaGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAAc         <  2:4308102/93‑1 (MQ=255)
                                                                                  gACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGc       >  2:3167234/1‑93 (MQ=255)
                                                                                  gACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGc       <  2:138653/93‑1 (MQ=255)
                                                                                    cccGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtg     >  1:1401491/1‑93 (MQ=255)
                                                                                     ccGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtgg    <  2:3558049/93‑1 (MQ=255)
                                                                                     ccGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtgg    >  1:2912863/1‑93 (MQ=255)
                                                                                      cgtgcgtgCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtggt   <  2:2095880/93‑1 (MQ=255)
                                                                                       gtgcgtgcGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtggtg  <  2:5084865/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCTGGTG  >  NC_000913/4067458‑4067637

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