Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,747,858 A→G Y48C (TAC→TGC)  sgbE → putative epimerase/aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,747,8580AG99.8% 1325.4 / ‑5.3 432Y48C (TAC→TGC) sgbEputative epimerase/aldolase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (179/252);  total (179/253)
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.70e-01

Only 100 of 469 total aligned reads displayed.

CTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACAC  >  NC_000913/3747768‑3747947
                                                                                          |                                                                                         
 tCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTgcga                                                                                        <  2:3052729/93‑1 (MQ=255)
 tCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTgcga                                                                                        <  1:4842969/93‑1 (MQ=255)
 tCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTgcga                                                                                        <  2:2374576/93‑1 (MQ=255)
 tCACCATCTGGTGACGTTCACCTGGGGCAATGGCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTgcga                                                                                        <  1:1233184/93‑1 (MQ=255)
  cACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTgcgac                                                                                       >  2:4909724/1‑93 (MQ=255)
   aCCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTgcgacg                                                                                      <  1:6678630/93‑1 (MQ=255)
    ccATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTgcgacgg                                                                                     >  2:1496466/1‑92 (MQ=255)
    ccATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGt                                                                                     >  2:2697253/1‑93 (MQ=255)
     cATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGtg                                                                                    >  2:4454000/1‑93 (MQ=255)
      aTCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGtga                                                                                   >  2:1127956/1‑93 (MQ=255)
      aTCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGtga                                                                                   <  1:5123089/93‑1 (MQ=255)
      aTCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGtga                                                                                   <  2:2620897/93‑1 (MQ=255)
         tGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGtgatga                                                                                <  2:3289857/93‑1 (MQ=255)
           gTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGAcc                                                                              >  1:2175769/1‑93 (MQ=255)
            tGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGAccg                                                                             >  1:1624755/1‑93 (MQ=255)
             gACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGAccgc                                                                            >  1:563687/1‑93 (MQ=255)
                gTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCcga                                                                         >  1:2276798/1‑93 (MQ=255)
                gTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCcga                                                                         <  2:1565944/93‑1 (MQ=255)
                    aCCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGat                                                                     <  1:7532566/93‑1 (MQ=255)
                     ccTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGata                                                                    <  1:5038270/93‑1 (MQ=255)
                     ccTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGata                                                                    <  2:1297428/93‑1 (MQ=255)
                      cTGGGGCAATGTCAGCGCTGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGatat                                                                   <  1:681142/93‑1 (MQ=255)
                      cTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGatat                                                                   <  1:6295842/93‑1 (MQ=255)
                      cTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGatat                                                                   <  1:2118070/93‑1 (MQ=255)
                      cTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGatat                                                                   <  2:34478/93‑1 (MQ=255)
                       tGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATAtg                                                                  <  2:4598040/93‑1 (MQ=255)
                       tGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATAtg                                                                  >  2:6190451/1‑93 (MQ=255)
                         gggCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATAtggt                                                                <  2:6853895/93‑1 (MQ=255)
                           gCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATAtggtgg                                                              >  1:1087640/1‑93 (MQ=255)
                           gCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATAtggtgg                                                              <  1:4741204/93‑1 (MQ=255)
                            cAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATAtggtggt                                                             <  1:3654857/93‑1 (MQ=255)
                             aaTGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATAtggtggtg                                                            >  1:7640099/1‑93 (MQ=255)
                                gTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGtt                                                         <  1:4364813/93‑1 (MQ=255)
                                gTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGtt                                                         <  1:5824603/93‑1 (MQ=255)
                                gTCAGCGCGGGAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGtt                                                         <  2:4030157/93‑1 (MQ=255)
                                   aGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTgag                                                      >  1:391406/1‑93 (MQ=255)
                                   aGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTgag                                                      <  1:1002155/93‑1 (MQ=255)
                                   aGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTgag                                                      >  2:3984770/1‑93 (MQ=255)
                                    gcgcGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTgaga                                                     >  2:7288536/1‑93 (MQ=255)
                                     ggcgGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGAt                                                    <  1:5941327/92‑1 (MQ=255)
                                     cgcgGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGAt                                                    <  2:5903202/93‑1 (MQ=255)
                                      gcgGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATa                                                   <  2:2915362/93‑1 (MQ=255)
                                      gcgGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATa                                                   <  2:3612528/93‑1 (MQ=255)
                                       cgGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAg                                                  >  1:5317957/1‑93 (MQ=255)
                                        ggTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGc                                                 <  2:879623/93‑1 (MQ=255)
                                          tAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCa                                               <  2:665414/93‑1 (MQ=255)
                                          tAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCa                                               >  1:6460399/1‑93 (MQ=255)
                                          tAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCa                                               <  2:2214607/93‑1 (MQ=255)
                                           aGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAg                                              >  1:6127054/1‑93 (MQ=255)
                                            gACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGc                                             >  1:939604/1‑93 (MQ=255)
                                            gACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGc                                             >  1:5768246/1‑93 (MQ=255)
                                              cGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCgg                                           <  1:2922552/93‑1 (MQ=255)
                                               gAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCGAGCgga                                          >  2:1878158/1‑92 (MQ=255)
                                               gAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGt                                          >  1:694648/1‑93 (MQ=255)
                                                aaaCGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGtg                                         <  1:1413622/93‑2 (MQ=255)
                                                aaaCGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACAGCCGACGATATGGAGGTGTTTGAGATAGCCAGCGGTa                                         >  2:3849004/1‑93 (MQ=38)
                                                 aaCGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTaa                                        >  1:3082719/1‑93 (MQ=255)
                                                 aaCGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTaa                                        <  1:2354876/93‑1 (MQ=255)
                                                 aaCGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTaa                                        <  1:5568484/93‑1 (MQ=255)
                                                  aCGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAg                                       >  2:3736686/1‑93 (MQ=255)
                                                  aCGCGGCAATGGATGGTAATCAAACCTGCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAg                                       <  1:4517689/93‑1 (MQ=255)
                                                     cgGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAggtg                                    <  2:3601795/93‑1 (MQ=255)
                                                      ggCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAggtgg                                   >  2:6278009/1‑93 (MQ=255)
                                                       gtaaTGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAggtggt                                  <  1:7732100/91‑1 (MQ=255)
                                                        cAATGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGgtggtg                                 <  1:7396854/93‑1 (MQ=255)
                                                         aaTGGATGGTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGtggtgg                                <  2:5422845/93‑1 (MQ=255)
                                                          atggatggAAATCAAACCTTCCGGCGTCGAGGGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGCGGTGGa                               >  2:5170484/1‑93 (MQ=37)
                                                           tggatggTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGag                              >  1:2992103/1‑92 (MQ=255)
                                                            ggatggTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAg                             <  2:7205302/93‑1 (MQ=255)
                                                              atggTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGgc                           >  2:7551172/1‑93 (MQ=255)
                                                              atggTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGgc                           <  1:5416542/93‑1 (MQ=255)
                                                              atggTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGgc                           >  1:2035907/1‑93 (MQ=255)
                                                               tggTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGgca                          <  1:462138/93‑1 (MQ=255)
                                                                 gTAATCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGgcagc                        <  1:2491244/93‑1 (MQ=255)
                                                                    aTCAAACCTTCCGGCGTCGAGGCCGACGTGATGACCACCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCaaa                     <  1:4002597/93‑1 (MQ=255)
                                                                     tCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCaaaa                    >  1:5893337/1‑93 (MQ=255)
                                                                     tCAAACCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCaaaa                    >  2:165039/1‑93 (MQ=255)
                                                                       aaaCCTTCCGGCGTCGTGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGGGGTGGAAGGCAGCaaaaaa                  >  2:7030478/1‑93 (MQ=255)
                                                                       aaaCCTTCCGGCGTCGAGTGCGAAGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCaaaaaa                  <  2:4321172/93‑1 (MQ=255)
                                                                        aaCCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAAc                 <  2:7692296/93‑1 (MQ=255)
                                                                        aaCCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAAc                 <  2:5024030/93‑1 (MQ=255)
                                                                         aCCTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAAcc                <  2:5598801/93‑1 (MQ=255)
                                                                          ccTTCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAATGGGGGGGAAGGCGGCAAAAAAccc               >  2:7770699/1‑93 (MQ=38)
                                                                            ttCCGGCGTCGAGTGCGACGTGATGCCCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCctc             <  2:4733856/93‑1 (MQ=255)
                                                                            ttCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCctc             <  1:2640611/93‑1 (MQ=255)
                                                                            ttCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCctc             <  2:2138372/93‑1 (MQ=255)
                                                                            ttCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCctc             <  2:5957921/93‑1 (MQ=255)
                                                                            ttCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGGGGTGGTTGAGATAGCCAGCGGTAAGGTGGGGGAAGGCAGCAAAAAACActc             >  2:180663/1‑93 (MQ=255)
                                                                             tCCGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCctct            <  1:7286128/93‑1 (MQ=255)
                                                                              ccGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCtt           <  2:6419756/93‑1 (MQ=255)
                                                                               cGGCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCTAGGGGTAAGGTGGTCGAAGGCAGCAAAAAACCCTCTTc          <  2:5302101/93‑1 (MQ=255)
                                                                                ggCGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTcc         <  1:3208529/93‑1 (MQ=255)
                                                                                  cGTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCATCAAAAAACCCTCTTCCGa       <  2:7924367/93‑1 (MQ=255)
                                                                                   gTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGAt      <  1:3839515/93‑1 (MQ=255)
                                                                                   gTCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGAt      >  1:6576458/1‑93 (MQ=255)
                                                                                    tCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATa     >  2:5762611/1‑93 (MQ=255)
                                                                                    tCGAGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATa     >  1:3387960/1‑93 (MQ=255)
                                                                                     cGAGTGCGACGTGATGACCGCCGACGATATGGTGGAGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATac    >  1:4244291/1‑93 (MQ=255)
                                                                                       aGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATacac  >  1:1847405/1‑93 (MQ=255)
                                                                                       aGTGCGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATacac  <  2:6037073/93‑1 (MQ=255)
                                                                                          |                                                                                         
CTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGGCAGCAAAAAACCCTCTTCCGATACAC  >  NC_000913/3747768‑3747947

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: ‑