Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,221,392 T→C R373R (CGT→CGC b3076 → evolved beta‑D‑galactosidase alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,221,3920TC84.8% 633.0 / 107.5 283R373R (CGT→CGCb3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (24/19);  new base C (115/125);  total (139/144)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.56e-01

Only 100 of 285 total aligned reads displayed.

CAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCAC  >  NC_000913/3221302‑3221482
                                                                                          |                                                                                          
   cGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGc                                                                                       <  1:973038/93‑1 (MQ=255)
    gATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCt                                                                                      >  1:4819815/1‑93 (MQ=255)
      tcatcGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCa                                                                                    >  2:4154377/1‑93 (MQ=255)
        atcGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACt                                                                                  <  2:5304335/93‑1 (MQ=255)
         tcGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTa                                                                                 >  2:3456151/1‑93 (MQ=255)
           gCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTAcc                                                                               >  1:772320/1‑93 (MQ=255)
           gCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTAcc                                                                               <  2:2643844/93‑1 (MQ=255)
            cAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTAccc                                                                              <  2:6800361/93‑1 (MQ=255)
             aaaGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCg                                                                             <  1:1274367/93‑1 (MQ=255)
               aGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGaa                                                                           <  1:71766/93‑1 (MQ=255)
               aGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGaa                                                                           >  1:2235844/1‑93 (MQ=255)
               aGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGaa                                                                           <  2:2233672/93‑1 (MQ=255)
               aGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGaa                                                                           <  1:3082697/93‑1 (MQ=255)
                ggCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAAc                                                                          <  1:2708746/93‑1 (MQ=255)
                ggCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAAc                                                                          <  1:101754/93‑1 (MQ=255)
                 gCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACg                                                                         >  2:6993302/1‑93 (MQ=255)
                     cgcCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATcc                                                                     <  1:2874038/93‑1 (MQ=255)
                     cgcCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATcc                                                                     <  2:4641784/93‑1 (MQ=255)
                     cgcCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATcc                                                                     >  2:663514/1‑93 (MQ=255)
                       ccGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCcgc                                                                   >  1:3738685/1‑93 (MQ=255)
                        cGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCcgcg                                                                  >  2:278537/1‑93 (MQ=255)
                         gTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGt                                                                 >  1:3188486/1‑93 (MQ=255)
                          ttGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGtt                                                                >  1:4507393/1‑93 (MQ=255)
                          ttGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGtt                                                                >  2:6539879/1‑93 (MQ=255)
                          ttGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGtt                                                                <  2:5178996/93‑1 (MQ=255)
                          ttGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGtt                                                                <  2:4722787/93‑1 (MQ=255)
                          ttGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGtt                                                                >  1:113419/1‑93 (MQ=255)
                          ttGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGtt                                                                >  1:1638057/1‑93 (MQ=255)
                          ttGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGtt                                                                >  1:1887528/1‑93 (MQ=255)
                            ggAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGtttt                                                              <  1:427703/93‑1 (MQ=255)
                              aaTGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTa                                                            <  2:2939984/93‑1 (MQ=255)
                               aTGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTAc                                                           >  2:1456201/1‑93 (MQ=255)
                               aTGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTAc                                                           >  2:4226202/1‑93 (MQ=255)
                                tGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACg                                                          >  2:7260390/1‑93 (MQ=255)
                                 ggATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGa                                                         >  2:6577188/1‑93 (MQ=255)
                                  gATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGaa                                                        <  1:3922223/93‑1 (MQ=255)
                                   aTCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAAc                                                       >  2:5072934/1‑93 (MQ=255)
                                   aTCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAAc                                                       >  2:213585/1‑93 (MQ=255)
                                    tCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACt                                                      >  1:5012908/1‑93 (MQ=255)
                                    tCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACt                                                      >  1:4097563/1‑93 (MQ=255)
                                    tCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACt                                                      >  1:5854127/1‑93 (MQ=255)
                                    tCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACt                                                      <  1:5942181/93‑1 (MQ=255)
                                     cgcgTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtg                                                     <  2:151035/93‑1 (MQ=255)
                                     cgcgTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtg                                                     <  2:3745427/93‑1 (MQ=255)
                                     cgcgTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtg                                                     <  1:7247913/93‑1 (MQ=255)
                                      gcgTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgt                                                    >  2:6079601/1‑93 (MQ=255)
                                      gcgTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgt                                                    >  2:4561337/1‑93 (MQ=255)
                                       cgTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgtg                                                   <  1:6834899/93‑1 (MQ=255)
                                       cgTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgtg                                                   <  1:6171283/93‑1 (MQ=255)
                                        gTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgtgt                                                  >  2:4041031/1‑93 (MQ=255)
                                        gTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgtgt                                                  <  2:786831/93‑1 (MQ=255)
                                          cGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGa                                                <  2:6003753/93‑1 (MQ=255)
                                           gagaAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGat                                               >  1:1744167/1‑93 (MQ=255)
                                            agaAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGata                                              <  2:2825918/93‑1 (MQ=255)
                                             gaAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGatat                                             >  1:3357010/1‑93 (MQ=255)
                                             gaAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGatat                                             >  1:281365/1‑93 (MQ=255)
                                              aaaGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATc                                            <  2:3301257/93‑1 (MQ=255)
                                               aaGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCt                                           >  1:5456846/1‑93 (MQ=255)
                                                aGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTa                                          <  2:2341280/93‑1 (MQ=255)
                                                aGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTa                                          >  1:2469379/1‑93 (MQ=255)
                                                 gATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTAc                                         >  2:4322710/1‑93 (MQ=255)
                                                 gATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTAc                                         >  2:4424019/1‑93 (MQ=255)
                                                  aTCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACg                                        <  1:2758902/93‑1 (MQ=255)
                                                   tctcCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACgg                                       <  2:1600229/93‑1 (MQ=255)
                                                   tctcCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACgg                                       >  1:542368/1‑93 (MQ=255)
                                                     tcCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGcc                                     >  2:6814032/1‑93 (MQ=255)
                                                      ccAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCt                                    >  2:180783/1‑93 (MQ=255)
                                                       cAGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTg                                   >  2:3068676/1‑93 (MQ=255)
                                                        aGTTGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGt                                  >  1:3481710/1‑93 (MQ=255)
                                                          ttGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGttt                                >  1:26233/1‑93 (MQ=255)
                                                          ttGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGttt                                <  1:2992360/93‑1 (MQ=255)
                                                          ttGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGttt                                >  2:3547830/1‑93 (MQ=255)
                                                          ttGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGttt                                <  1:4180470/93‑1 (MQ=255)
                                                          ttGATGAAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCATGctt                                >  1:6969452/1‑93 (MQ=255)
                                                           tgatgaAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTtg                               <  1:5069585/93‑1 (MQ=255)
                                                           tgatgaAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTtg                               <  1:3649966/93‑1 (MQ=255)
                                                            gatgaAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTtgt                              <  2:5880972/93‑1 (MQ=255)
                                                              tgaAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGa                            >  2:7137357/1‑93 (MQ=255)
                                                               gaAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGAt                           <  1:3591814/93‑1 (MQ=255)
                                                               gaAGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGAt                           <  1:6660913/93‑1 (MQ=255)
                                                                aaGCAGCACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATg                          >  2:1280344/1‑93 (MQ=255)
                                                                 agcagcACAATATCAACTCCGTGCGTCCCGCTCACTACCCGAGCGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATgg                         <  1:5339234/93‑1 (MQ=255)
                                                                 agcagcACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATgg                         <  2:2693684/93‑1 (MQ=255)
                                                                 agcagcACAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATgg                         <  2:615550/93‑1 (MQ=255)
                                                                  gcagcaCAATATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGc                        <  1:4973312/93‑1 (MQ=255)
                                                                      cacaATATCAACTCCGTGCGCACCGCTAACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGaa                    >  2:6430220/1‑93 (MQ=255)
                                                                         aaTATCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAAcc                 >  1:6294486/1‑93 (MQ=255)
                                                                           tatCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGa               <  2:3379499/93‑1 (MQ=255)
                                                                            atCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGAc              <  1:846383/93‑1 (MQ=255)
                                                                             tCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACg             >  1:1512061/1‑93 (MQ=255)
                                                                             tCAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACg             <  1:5689398/93‑1 (MQ=255)
                                                                              cAACTCCGTGCGCGCCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGGGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGt            <  1:4756583/93‑1 (MQ=255)
                                                                              cAACTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGt            <  2:3676296/93‑1 (MQ=255)
                                                                               aaCTCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTc           <  1:4637817/93‑1 (MQ=255)
                                                                                  tCCGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGaa        >  1:4341675/1‑93 (MQ=255)
                                                                                   ccGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAAt       >  2:5751611/1‑93 (MQ=255)
                                                                                   ccGTGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAAt       <  2:6784819/93‑1 (MQ=255)
                                                                                      tGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGc    >  2:855828/1‑93 (MQ=255)
                                                                                      tGCGCACCGCTCACTACCCGAACGATCCGCGTTTTTAAGAACTGTGAGATATCCACGGCCTGTTTGTCATGGCGGTAACCCGCGGCGAATCGc    >  1:2043122/1‑93 (MQ=255)
                                                                                       gcgcACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCa   >  1:5762750/1‑93 (MQ=255)
                                                                                          |                                                                                          
CAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCAC  >  NC_000913/3221302‑3221482

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