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% 825.0 / ‑5.2 257H137R (CAC→CGC) b3879putative aldose‑1‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (122/132);  total (124/133)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 283 total aligned reads displayed.

TTGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCTGGTG  >  NC_000913/4067458‑4067637
                                                                                          |                                                                                         
  gCAGGTTAAATAAGGATGAATAACCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGGAAGACCCgtgcgtgc                                                                                       >  1:4823323/1‑93 (MQ=255)
   cAGGTTAAATAAGGATGAATATCCACGCCATAAGGCGCCGCGACAGTGCCGATATTCTGGCTTGGGATCTCAAAGGAAAGACCCGtgcgtgcg                                                                                      >  2:4495444/1‑93 (MQ=255)
   cAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGtgcgtgcg                                                                                      <  2:4294916/93‑1 (MQ=255)
   cAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACACGtgcgtgcg                                                                                      >  2:1841546/1‑93 (MQ=255)
   aaGGTTAAATAAGTATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGtgcgtgcg                                                                                      <  1:5454883/92‑1 (MQ=255)
    aGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCCACGGAAAGACCCGTgcgtgcgt                                                                                     >  2:1090755/1‑93 (MQ=255)
     ggTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGtt                                                                                    >  1:3976153/1‑93 (MQ=255)
       ttAAATAAGGATGAATACCCCCGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGa                                                                                  <  1:4299213/93‑1 (MQ=255)
       ttAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGa                                                                                  <  2:4830275/93‑1 (MQ=255)
       ttAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGa                                                                                  >  2:6322022/1‑93 (MQ=255)
       ttAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGa                                                                                  >  1:3365184/1‑93 (MQ=255)
       ttAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGa                                                                                  >  1:3766992/1‑93 (MQ=255)
        tAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGGGTTGAg                                                                                 >  2:440164/1‑93 (MQ=255)
        tAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAg                                                                                 <  2:1360657/93‑1 (MQ=255)
         aaaTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGc                                                                                <  2:1475321/93‑1 (MQ=255)
          aaTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCg                                                                               >  1:2032782/1‑93 (MQ=255)
          aaTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCg                                                                               >  2:1278895/1‑93 (MQ=255)
          aaTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCg                                                                               >  1:5400556/1‑93 (MQ=255)
          aaTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCg                                                                               <  2:644575/93‑1 (MQ=255)
           aTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGa                                                                              >  1:1162130/1‑93 (MQ=255)
           aTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGa                                                                              >  2:5735391/1‑93 (MQ=255)
           aTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGa                                                                              >  2:6603790/1‑93 (MQ=255)
           aTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGa                                                                              >  1:170583/1‑93 (MQ=255)
              aGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATa                                                                           >  1:5884801/1‑93 (MQ=255)
              aGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATa                                                                           >  2:3044187/1‑93 (MQ=255)
              aGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATa                                                                           >  1:909304/1‑93 (MQ=255)
               ggATGAATCCCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATaa                                                                          <  2:3746295/93‑1 (MQ=255)
                gATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATaaa                                                                         >  2:6626219/1‑93 (MQ=255)
                 aTGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATaaaa                                                                        >  2:5284320/1‑92 (MQ=255)
                 aTGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAaac                                                                        <  1:2623435/93‑1 (MQ=255)
                  tGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATCaaca                                                                       >  2:6265176/1‑93 (MQ=255)
                     aTACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACt                                                                    >  1:6230814/1‑93 (MQ=255)
                      tACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACtt                                                                   >  1:4240370/1‑93 (MQ=255)
                      tACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACtt                                                                   >  2:2280548/1‑93 (MQ=255)
                       aCCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtg                                                                  >  1:361159/1‑93 (MQ=255)
                        cccACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtgt                                                                 <  1:6590194/93‑1 (MQ=255)
                          cACGCCATAAGGCGCTGCGATAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGa                                                               >  1:6455349/1‑93 (MQ=255)
                           aCGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGaa                                                              >  2:2501385/1‑93 (MQ=255)
                           aCGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCATGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGaa                                                              <  1:2317257/93‑1 (MQ=255)
                            cGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAg                                                             >  2:4041549/1‑93 (MQ=255)
                             gCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGc                                                            <  2:4181065/93‑1 (MQ=255)
                                   aGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCAt                                                      <  1:3540386/93‑1 (MQ=255)
                                   aGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCAt                                                      >  1:5762522/1‑93 (MQ=255)
                                   aGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCAt                                                      >  1:5697407/1‑93 (MQ=255)
                                   aGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCAt                                                      >  1:3150111/1‑93 (MQ=255)
                                    ggCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATa                                                     >  1:3484437/1‑93 (MQ=255)
                                    ggCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCCAATAAACAACTTGTGAAGCCAGCATa                                                     >  1:1558823/1‑93 (MQ=255)
                                     gcgcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaa                                                    <  2:5544422/93‑1 (MQ=255)
                                     gcgcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaa                                                    >  2:6390948/1‑93 (MQ=255)
                                      cgcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaaa                                                   <  1:5200053/93‑1 (MQ=255)
                                      cgcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaaa                                                   <  1:5604700/93‑1 (MQ=255)
                                      cgcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaaa                                                   <  1:1316982/93‑1 (MQ=255)
                                       gcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAAc                                                  <  1:6081957/93‑1 (MQ=255)
                                       gcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAAc                                                  <  1:2480911/93‑1 (MQ=255)
                                        cTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACg                                                 <  1:4677728/93‑1 (MQ=255)
                                          gCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGa                                               >  2:3063030/1‑93 (MQ=255)
                                           cGACAGTACCCATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTGGAGCGAATAAACATCTTGTGAAGCCAGCATAAACGGAt                                              <  2:768952/93‑1 (MQ=255)
                                            gACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATa                                             <  1:805063/93‑1 (MQ=255)
                                            gACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATa                                             <  2:3433808/93‑1 (MQ=255)
                                            gACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATa                                             <  1:6663610/93‑1 (MQ=255)
                                             aCAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAg                                            <  2:4595029/93‑1 (MQ=255)
                                               aGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGcc                                          <  1:1767614/93‑1 (MQ=255)
                                                gTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCg                                         <  2:6050097/93‑1 (MQ=255)
                                                  gCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTa                                       >  2:5224221/1‑93 (MQ=255)
                                                  gCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTa                                       <  2:2210505/93‑1 (MQ=255)
                                                  gCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTa                                       >  2:2717855/1‑93 (MQ=255)
                                                    cGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAAc                                     <  2:3509758/93‑1 (MQ=255)
                                                    cGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAAc                                     <  2:1994503/93‑1 (MQ=255)
                                                    cGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAAc                                     >  1:5054783/1‑93 (MQ=255)
                                                    cGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAAc                                     >  2:4783317/1‑93 (MQ=255)
                                                     gATATTCTGGCTTGCGATCTTAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACt                                    >  2:3922076/1‑93 (MQ=255)
                                                     gATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACt                                    >  2:1629374/1‑93 (MQ=255)
                                                     gATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACt                                    >  2:721144/1‑93 (MQ=255)
                                                       tatTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTgg                                  >  1:2993603/1‑93 (MQ=255)
                                                            tGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGggcggc                             >  1:2061151/1‑93 (MQ=255)
                                                             ggCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCa                            <  2:1010324/93‑1 (MQ=255)
                                                                ttGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGaa                         >  2:2067663/1‑93 (MQ=255)
                                                                ttGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGaa                         >  1:3260719/1‑93 (MQ=255)
                                                                      actcAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCTAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTg                   <  1:5438962/92‑1 (MQ=255)
                                                                       ctcAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGt                  <  2:2773544/93‑1 (MQ=255)
                                                                         cAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAACCAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCa                <  2:816367/93‑1 (MQ=255)
                                                                         cAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCa                <  2:6137115/93‑1 (MQ=255)
                                                                           aCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAgg              >  1:4160271/1‑93 (MQ=255)
                                                                           aCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAgg              >  1:3189810/1‑93 (MQ=255)
                                                                           aCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAgg              >  1:1455422/1‑93 (MQ=255)
                                                                           aCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGACAGCCGTAACTGGGCGGCAGAAAAGCTGTCAgg              <  1:1009861/93‑1 (MQ=255)
                                                                            cGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAggg             <  1:2531415/93‑1 (MQ=255)
                                                                             ggAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGt            <  1:5381139/93‑1 (MQ=255)
                                                                              gAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTa           >  2:6207844/1‑93 (MQ=255)
                                                                              gAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTa           <  1:4352572/93‑1 (MQ=255)
                                                                               aaaGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTaa          <  2:407656/93‑1 (MQ=255)
                                                                                aaGACCCGTGCGTGCGTTGAGCGAATAAACAAATTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAAc         >  2:562984/1‑93 (MQ=255)
                                                                                 aGACCCGTGCGTGCGTTGAGCGACTAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTTGGCGGCAGAAAAGCTGTCAGGGTAACg        >  1:841705/1‑93 (MQ=255)
                                                                                  gACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGc       >  2:4620982/1‑93 (MQ=255)
                                                                                  gACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGc       <  2:1181954/93‑1 (MQ=255)
                                                                                   aCCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCt      <  2:1552533/93‑1 (MQ=255)
                                                                                    cccGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtg     <  1:4523815/93‑1 (MQ=255)
                                                                                     ccGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtgg    >  1:443606/1‑93 (MQ=255)
                                                                                     ccGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtgg    <  1:5981336/93‑1 (MQ=255)
                                                                                       gtgcgtgcGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtggtg  >  2:883968/1‑93 (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: