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.7% 1043.2 / ‑5.0 336H137R (CAC→CGC) b3879putative aldose‑1‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (175/159);  total (176/160)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

Only 100 of 376 total aligned reads displayed.

TTGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCTGGTG  >  NC_000913/4067458‑4067637
                                                                                          |                                                                                         
ttGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCgtgcgt                                                                                         <  1:4607035/93‑1 (MQ=255)
ttGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCgtgcgt                                                                                         >  2:2917237/1‑93 (MQ=255)
 tGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGAGATCTCAACGGAAAGACCcgtgcgtg                                                                                        >  1:3009049/1‑93 (MQ=255)
     ggTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGtt                                                                                    >  1:167123/1‑93 (MQ=255)
     ggTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGtt                                                                                    >  2:6231627/1‑93 (MQ=255)
      gTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTg                                                                                   >  2:5851241/1‑93 (MQ=255)
       ttAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGa                                                                                  >  1:754251/1‑93 (MQ=255)
        tAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAg                                                                                 >  2:4151012/1‑93 (MQ=255)
         aaaTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGc                                                                                >  2:3562491/1‑93 (MQ=255)
           aTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGa                                                                              >  2:100233/1‑93 (MQ=255)
            tAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGaa                                                                             <  1:1208423/93‑1 (MQ=255)
            tAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGaa                                                                             >  1:2823649/1‑93 (MQ=255)
            tAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGaa                                                                             >  2:4731082/1‑93 (MQ=255)
             aaGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAAt                                                                            >  2:6942402/1‑93 (MQ=255)
              aGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATa                                                                           >  2:5486076/1‑93 (MQ=255)
              aGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATa                                                                           >  1:4090446/1‑93 (MQ=255)
              aGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATa                                                                           <  1:2026553/93‑1 (MQ=255)
              aGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATa                                                                           <  1:6961416/93‑1 (MQ=255)
              aGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATa                                                                           <  2:5326987/93‑1 (MQ=255)
                gATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATaaa                                                                         <  2:3660653/93‑1 (MQ=255)
                 aTGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAaac                                                                        <  2:1843139/93‑1 (MQ=255)
                 aTGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAaac                                                                        >  1:3668325/1‑93 (MQ=255)
                  tGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAaaca                                                                       >  1:6462901/1‑93 (MQ=255)
                    aaTACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTTCGTGCGTTGAGCGAATAaacaac                                                                     >  2:4314167/1‑93 (MQ=255)
                    aaTACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAaacaac                                                                     <  1:5257116/93‑1 (MQ=255)
                     aTACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACt                                                                    <  1:1925084/93‑1 (MQ=255)
                      tACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACtt                                                                   >  2:3964713/1‑93 (MQ=255)
                        cccACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtgt                                                                 <  1:5887467/93‑1 (MQ=255)
                        cccACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtgt                                                                 >  2:3468286/1‑93 (MQ=255)
                        cccACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtgt                                                                 >  1:6643797/1‑93 (MQ=255)
                         ccACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtgtg                                                                >  1:4116814/1‑93 (MQ=255)
                         ccACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtgtg                                                                >  1:2140693/1‑93 (MQ=255)
                          cACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGa                                                               >  2:5352687/1‑93 (MQ=255)
                          cACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGa                                                               >  2:4870728/1‑93 (MQ=255)
                           cgcccATAAGGCGCCGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGaa                                                              <  1:7336580/90‑1 (MQ=255)
                           aCGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGaa                                                              >  2:2787552/1‑93 (MQ=255)
                           aCGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGaa                                                              <  2:66496/93‑1 (MQ=255)
                           aCGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAAAAACTTGTGaa                                                              >  2:6044192/1‑93 (MQ=255)
                            cGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAg                                                             <  2:2011135/93‑1 (MQ=255)
                             gCCATAATGCGCTGCGACAGTGCCGATATTCTGGCTTGCAATCTCAACGGAAAGACCCGTGCGCGCGTTGAGCGAATACACCACTTGTGAAGc                                                            >  2:7936386/1‑93 (MQ=255)
                             gCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGc                                                            <  1:228003/93‑1 (MQ=255)
                              ccATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGcc                                                           >  1:1526506/1‑93 (MQ=255)
                               cATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCa                                                          <  1:5973902/93‑1 (MQ=255)
                               cATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCa                                                          <  2:90729/93‑1 (MQ=255)
                               cATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCa                                                          >  1:3817803/1‑93 (MQ=255)
                                aTAAGGCGCTGCGACAGTGCCGATATTCTGGTTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAg                                                         <  1:2885271/93‑1 (MQ=255)
                                 tAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGc                                                        >  1:3111792/1‑93 (MQ=255)
                                 tAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGc                                                        >  1:1794589/1‑93 (MQ=255)
                                 tAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGc                                                        >  1:1014842/1‑93 (MQ=255)
                                  aaGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGACAGCa                                                       >  2:373782/1‑93 (MQ=255)
                                   aGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCAt                                                      >  1:5136583/1‑93 (MQ=255)
                                   aGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCAt                                                      <  1:1376414/93‑1 (MQ=255)
                                   aGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCAt                                                      <  2:1127908/93‑1 (MQ=255)
                                      cgcTGCGACAGGGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaaa                                                   <  1:3394414/93‑1 (MQ=255)
                                       gcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAAc                                                  >  1:1116719/1‑93 (MQ=255)
                                       gcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAAc                                                  <  1:4395785/93‑1 (MQ=255)
                                        cTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACg                                                 <  1:6140933/93‑1 (MQ=255)
                                         tGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACgg                                                <  2:7324714/93‑1 (MQ=255)
                                          gCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGa                                               <  1:2602067/93‑1 (MQ=255)
                                            gACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATa                                             <  2:2617144/93‑1 (MQ=255)
                                              cAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGc                                           <  1:2515096/93‑1 (MQ=255)
                                              cAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGc                                           >  2:5709863/1‑93 (MQ=255)
                                               aGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGcc                                          <  1:5520566/93‑1 (MQ=255)
                                                gTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCg                                         >  1:3299732/1‑93 (MQ=255)
                                                 tGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGt                                        >  2:4598334/1‑93 (MQ=255)
                                                   ccGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTaa                                      <  2:6360524/93‑1 (MQ=255)
                                                    cGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAAc                                     <  1:3245824/93‑1 (MQ=255)
                                                       tatTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTgg                                  <  2:2317265/93‑1 (MQ=255)
                                                       tatTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTgg                                  <  1:4834043/93‑1 (MQ=255)
                                                        atTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTggg                                 >  1:2341668/1‑93 (MQ=255)
                                                        atTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTggg                                 <  2:1266511/93‑1 (MQ=255)
                                                         ttCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGggc                                >  2:7617176/1‑93 (MQ=255)
                                                             ggCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCa                            >  1:5730591/1‑93 (MQ=255)
                                                               cTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGa                          >  1:5076503/1‑93 (MQ=255)
                                                                ttGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGaa                         >  1:7785176/1‑93 (MQ=255)
                                                                 tGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGaaa                        <  1:857685/93‑1 (MQ=255)
                                                                 tGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGaaa                        <  2:7782084/93‑1 (MQ=255)
                                                                 tGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGaaa                        <  2:1605173/93‑1 (MQ=255)
                                                                   cGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAg                      >  1:3451426/1‑93 (MQ=255)
                                                                    gATCTCAACGGAAAGACCCGTGCGTGCGTTGAGGGAATAAACAACTTGTGCAGACAGGATAAACGGATAGTCGAAACCGGGCGGCAGAAAAGc                     >  1:4101987/1‑93 (MQ=255)
                                                                      tctcAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTg                   >  1:6855459/1‑93 (MQ=255)
                                                                      actcAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTg                   <  2:6835532/92‑1 (MQ=255)
                                                                       ctcAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGt                  >  2:3779801/1‑93 (MQ=255)
                                                                         cAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCa                >  2:5552068/1‑93 (MQ=255)
                                                                         cAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCa                <  1:7925305/93‑1 (MQ=255)
                                                                          aaCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAg               >  2:754931/1‑93 (MQ=255)
                                                                          aaCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAg               >  1:3534182/1‑93 (MQ=255)
                                                                           aCGGAAATACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAgg              <  2:2215315/93‑1 (MQ=255)
                                                                           aCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAgg              >  2:4285786/1‑93 (MQ=255)
                                                                            cgaaaaGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAggg             <  2:6308884/90‑1 (MQ=255)
                                                                             ggAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGt            >  1:7714782/1‑93 (MQ=255)
                                                                              gAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTa           >  2:5095356/1‑93 (MQ=255)
                                                                               aaaGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTaa          <  1:4182726/93‑1 (MQ=255)
                                                                                aaGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAAc         >  2:706312/1‑93 (MQ=255)
                                                                                aaGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAAc         <  1:4906067/93‑1 (MQ=255)
                                                                                 aGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACg        >  1:5617439/1‑93 (MQ=255)
                                                                                  gACCCGTGCGTGCGTTGAGCGAATAAACCACTTGTGAAACCAGCATCAACGGATAGCCGTAACTGGTCCGCAGACAAGCTGTCAGCGTacagc       >  1:5389375/1‑89 (MQ=255)
                                                                                    cccGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtg     >  1:6272610/1‑93 (MQ=255)
                                                                                    cccGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtg     <  2:3219791/93‑1 (MQ=255)
                                                                                      cgtgcgtgCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtgg  >  2:6518380/1‑92 (MQ=255)
                                                                                          |                                                                                         
TTGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCTGGTG  >  NC_000913/4067458‑4067637

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