Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,548,005 C→T G701D (GGT→GAT)  malP ← maltodextrin phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,548,0050CT99.5% 782.5 / NA 212G701D (GGT→GAT) malPmaltodextrin phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (99/112);  total (99/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

ATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTACCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAGATATTTTCT  >  NC_000913/3547915‑3548094
                                                                                          |                                                                                         
aTACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCaccttatct                                                                                         <  2:2993677/93‑1 (MQ=255)
 tACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCaccttatctt                                                                                        <  1:4835752/93‑1 (MQ=255)
 tACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCaccttatctt                                                                                        <  1:3771620/93‑1 (MQ=255)
   cTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTttc                                                                                      <  2:2730790/93‑1 (MQ=255)
   cTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTttc                                                                                      <  2:4835534/93‑1 (MQ=255)
   cTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTttc                                                                                      <  1:2598564/93‑1 (MQ=255)
    tgtgCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTttct                                                                                     <  2:2243241/93‑1 (MQ=255)
     gtgCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTAtctttctt                                                                                    >  1:3492395/1‑93 (MQ=255)
       gcagcaTCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTcc                                                                                  >  1:4399665/1‑93 (MQ=255)
       gcagcaTCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTcc                                                                                  >  1:3251384/1‑93 (MQ=255)
         agcaTCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGc                                                                                >  2:2161378/1‑93 (MQ=255)
          gcaTCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGcc                                                                               <  1:3661666/93‑1 (MQ=255)
           caTCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCa                                                                              <  2:5560055/93‑1 (MQ=255)
             tCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCAtt                                                                            >  1:3816876/1‑93 (MQ=255)
             tCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCAtt                                                                            >  2:4198716/1‑93 (MQ=255)
                ggTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCa                                                                         >  2:3420754/1‑93 (MQ=255)
                 gTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCAc                                                                        >  1:4310474/1‑93 (MQ=255)
                 gTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCAc                                                                        <  2:5643656/93‑1 (MQ=255)
                 gTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCAc                                                                        <  1:3208630/93‑1 (MQ=255)
                 gTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCAc                                                                        <  2:2815297/93‑1 (MQ=255)
                   cGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCg                                                                      <  2:5254376/93‑1 (MQ=255)
                    gAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCgg                                                                     >  2:5738085/1‑93 (MQ=255)
                     aaGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCggg                                                                    >  2:896965/1‑93 (MQ=255)
                     aaGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCggg                                                                    >  1:2202152/1‑93 (MQ=255)
                      aGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTGCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGt                                                                   <  1:4975875/93‑1 (MQ=255)
                       ggCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTc                                                                  >  2:3528788/1‑93 (MQ=255)
                        gCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCg                                                                 >  2:2475229/1‑93 (MQ=255)
                         cATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGt                                                                <  1:216222/93‑1 (MQ=255)
                           tGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAg                                                              >  2:5390078/1‑93 (MQ=255)
                           tGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAg                                                              >  1:2043895/1‑93 (MQ=255)
                            gCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGc                                                             <  1:2291349/93‑1 (MQ=255)
                              ttATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCt                                                           >  1:4581198/1‑93 (MQ=255)
                               tATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCtt                                                          <  1:1554812/93‑1 (MQ=255)
                                aTCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCttt                                                         <  2:6060965/93‑1 (MQ=255)
                                aTCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCttt                                                         <  1:4409626/93‑1 (MQ=255)
                                aTCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCttt                                                         <  1:2459235/93‑1 (MQ=255)
                                aTCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCttt                                                         >  2:4658617/1‑93 (MQ=255)
                                 tCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTg                                                        <  1:3575902/93‑1 (MQ=255)
                                 tCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTg                                                        <  1:345856/93‑1 (MQ=255)
                                 tCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTg                                                        <  1:5903318/93‑1 (MQ=255)
                                  cGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTgg                                                       <  1:544648/93‑1 (MQ=255)
                                         gCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGGGTCCAGCCCTTAACCTTCCTTCGGCCATTTCACCGGGTCGTACCTTTTGGCAAGAAt                                                >  2:5178089/1‑93 (MQ=255)
                                           tGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCGTTTCGCCAGAAttc                                              >  1:2721506/1‑92 (MQ=255)
                                            gTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGc                                             <  2:4510363/93‑1 (MQ=255)
                                            gTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGc                                             >  1:1892568/1‑93 (MQ=255)
                                             tATTTACCGCTTGCCAGCTCTTTCAATCCTTCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGcc                                            <  2:3140806/93‑1 (MQ=255)
                                               tttACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCtt                                          >  1:1311586/1‑93 (MQ=255)
                                               tttACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCtt                                          >  1:2444086/1‑93 (MQ=255)
                                               tttACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCtt                                          >  1:4656515/1‑93 (MQ=255)
                                                ttACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTc                                         <  1:2884044/93‑1 (MQ=255)
                                                ttACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTc                                         <  2:5595834/93‑1 (MQ=255)
                                                 tACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCa                                        <  2:4426787/93‑1 (MQ=255)
                                                 tACCGCTTCCCAGCTCTTTCAACACTGCGTCCAGCACCTTATCTTTCTCCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCa                                        <  1:3080560/93‑1 (MQ=255)
                                                  aCCGCTTTCCAACTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCAc                                       >  1:3686419/1‑93 (MQ=255)
                                                   ccGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACt                                      <  2:3481808/93‑1 (MQ=255)
                                                     gCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTg                                    <  1:651100/93‑1 (MQ=255)
                                                       tttCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTggt                                  >  1:6005173/1‑91 (MQ=255)
                                                        ttCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTc                                 >  2:5946203/1‑93 (MQ=255)
                                                         tCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTcc                                <  2:196254/93‑1 (MQ=255)
                                                         tCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTcc                                <  1:4496229/93‑1 (MQ=255)
                                                            aGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACg                             <  2:4978293/93‑1 (MQ=255)
                                                            aGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACg                             <  1:314910/93‑1 (MQ=255)
                                                             gCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACgg                            >  2:4282070/1‑93 (MQ=255)
                                                             gCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACgg                            <  2:3936091/93‑1 (MQ=255)
                                                             gCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACgg                            <  1:551294/93‑1 (MQ=255)
                                                              ctctTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGt                           >  1:1832861/1‑93 (MQ=255)
                                                              ctctTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGt                           <  2:5873013/93‑1 (MQ=255)
                                                              ctctTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGt                           <  2:3377551/93‑1 (MQ=255)
                                                               tctTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTa                          >  1:1296744/1‑93 (MQ=255)
                                                                ctTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTAt                         <  1:4049198/93‑1 (MQ=255)
                                                                ctTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTAt                         <  1:5579598/93‑1 (MQ=255)
                                                                  ttCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGa                       <  2:107249/93‑1 (MQ=255)
                                                                   tCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGAc                      <  1:1116680/93‑1 (MQ=255)
                                                                   tCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGAc                      >  2:1830186/1‑93 (MQ=255)
                                                                    cAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGAcc                     >  2:5323901/1‑93 (MQ=255)
                                                                     aaTACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCa                    >  1:4223709/1‑93 (MQ=255)
                                                                      aTACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCaa                   >  1:275124/1‑93 (MQ=255)
                                                                       tACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCaaa                  >  1:1067301/1‑93 (MQ=255)
                                                                        aCTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCaaaa                 <  1:475747/93‑1 (MQ=255)
                                                                         cTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCaaaaa                <  2:3704229/93‑1 (MQ=255)
                                                                         cTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCaaaaa                <  1:767629/93‑1 (MQ=255)
                                                                          tGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAAt               >  2:2692122/1‑93 (MQ=255)
                                                                           gCGTCCAGCCCCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATa              <  2:1216114/93‑1 (MQ=255)
                                                                             gTCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATaaa            <  1:1695665/93‑1 (MQ=255)
                                                                              tCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAg           >  2:612358/1‑93 (MQ=255)
                                                                              tCCAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAg           >  2:2602941/1‑93 (MQ=255)
                                                                                cAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAACATAAAGat         >  2:1365965/1‑93 (MQ=255)
                                                                                cAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAGat         <  2:2147828/93‑1 (MQ=255)
                                                                                cAGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAGat         <  1:4003549/93‑1 (MQ=255)
                                                                                 aGCACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAGata        <  2:4104172/93‑1 (MQ=255)
                                                                                   cACCTTATCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAGATAtt      <  2:2503120/93‑1 (MQ=255)
                                                                                    accttatcttTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAGATAttt     <  1:3273501/93‑1 (MQ=255)
                                                                                    accttatcttTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAGATAttt     <  2:4711743/93‑1 (MQ=255)
                                                                                    accttatcttTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAGATAttt     >  1:1590258/1‑93 (MQ=255)
                                                                                     ccttatcttTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAGATAtttt    <  2:2230771/93‑1 (MQ=255)
                                                                                     ccttatcttTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAGATAtttt    >  1:5771431/1‑93 (MQ=255)
                                                                                     ccttatcttTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAGATAtttt    >  1:5323941/1‑93 (MQ=255)
                                                                                      cttatcttTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAGATATTttc   <  1:1206065/93‑1 (MQ=255)
                                                                                      cttatcttTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAGATATTttc   >  2:72837/1‑93 (MQ=255)
                                                                                       ttatcttTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAGATATTttct  <  2:1549161/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTACCTTTCTTCCGCCATTTCACCGGGTCGTAGCCTTTGGCCAGAATTGCCTTCACTTGTTCCACGGTATGACCAAAAATAAAGATATTTTCT  >  NC_000913/3547915‑3548094

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: