Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,263,892 A→G 5.2% W376R (TGG→CGG)  tdcC ← L‑threonine/L‑serine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,263,8920AG 5.2% 225.7 / 2.7 78W376R (TGG→CGG) tdcCL‑threonine/L‑serine transporter
Reads supporting (aligned to +/- strand):  ref base A (29/43);  new base G (2/2);  total (32/46)
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.64e-01

CGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGTTCCCAGATAGTGACCGAAGAAAGATTTGAAGATA  >  NC_000913/3263745‑3264040
                                                                                                                                                   |                                                                                                                                                    
cGTTATCCATACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATTCATCGGCAACAGGCACAGCAGGGATGCGATAATCGGGGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAg                                                                                                                                                     <  1:1013249/149‑1 (MQ=255)
 gTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAgg                                                                                                                                                    >  1:1001994/1‑149 (MQ=255)
  ttATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGt                                                                                                                                                   <  2:1001994/149‑1 (MQ=255)
    aTCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCCACAACCCAGGTgg                                                                                                                                                 >  2:383716/1‑149 (MQ=255)
    aTCCAGACGACCACGGTATATCGCCAGAGACTGAACTTTACGGATGGCATACATAGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGCCTACAATCATGTCAAGCATGTTCGGTTTTGCGTAGGCAACAACCCAGGTgg                                                                                                                                                 >  1:1027961/1‑149 (MQ=255)
        aGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCGGGTGGAGcc                                                                                                                                             >  1:132756/1‑149 (MQ=255)
        aGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCGGGTGGAGcc                                                                                                                                             >  1:662572/1‑149 (MQ=255)
        aGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGcc                                                                                                                                             >  1:96850/1‑149 (MQ=255)
         gacgacCACGGCAATTCGCCAGAGACGGCGCTTTACGGCTGGCATACCTCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGccc                                                                                                                                            <  1:508323/149‑1 (MQ=255)
           cgcccccGGTATTTCGCCACAGACGGCGCTTTAAGGATGGCAGACATCGGCAACACGCACAGCAGGGATGCGATAATCGCTGCGCCCATGGCTTCAATCAGGGCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCat                                                                                                                                          <  2:751889/143‑1 (MQ=255)
            gacCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCatg                                                                                                                                         >  1:878007/1‑149 (MQ=255)
              ccACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCatgat                                                                                                                                       <  2:1025955/149‑1 (MQ=255)
              ccACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGAGGAGCCCatgat                                                                                                                                       >  2:616593/1‑149 (MQ=255)
               cACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCatgatg                                                                                                                                      >  1:802209/1‑149 (MQ=255)
               cACGGTATTTCGCCAGAGAAGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCCGGTGGAGCCCatgatg                                                                                                                                      >  2:997523/1‑149 (MQ=255)
                    tATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCGTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGAt                                                                                                                                 <  2:279989/149‑1 (MQ=255)
                      tttCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGCTCa                                                                                                                               >  2:49924/1‑149 (MQ=255)
                      tttCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGAAAATCGGTGCGCCCCTGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACCACCCAGGTGGAGCCCATGATGAAGCTCa                                                                                                                               >  1:992072/1‑149 (MQ=255)
                         cGCCAGAGACGGCGCTTTACGGAGGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGc                                                                                                                            <  2:307327/149‑1 (MQ=255)
                                 aCGGCGCTTTACGGATGGCATACATCCGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAgtg                                                                                                                    <  2:701903/149‑1 (MQ=255)
                                   ggCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGtt                                                                                                                  <  1:1025758/149‑1 (MQ=255)
                                   ggCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGtt                                                                                                                  <  1:172360/149‑1 (MQ=255)
                                     cgcTTTCCGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCa                                                                                                                <  1:617657/149‑1 (MQ=255)
                                      gcTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAg                                                                                                               <  2:160324/149‑1 (MQ=255)
                                       cTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGt                                                                                                              <  1:90216/149‑1 (MQ=255)
                                        tttACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGtt                                                                                                             >  1:146853/1‑149 (MQ=255)
                                        gttACGGCTGCCATACATCGGCAACAGGCCCAGCAGGGATGCGATACTCTGTGCGCCCCTGGCGTCAATCAGGTCAAGGATGTGCGGGTTGTCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGtt                                                                                                             <  2:906253/148‑1 (MQ=255)
                                           aCGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTAc                                                                                                          <  2:894763/149‑1 (MQ=255)
                                                  gCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTGAAGGATGTTCGGTTTGGCGTAGGCAACAACCCTGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGa                                                                                                   <  1:616593/149‑1 (MQ=255)
                                                   cATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGac                                                                                                  >  2:612134/1‑149 (MQ=255)
                                                    aTACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGaca                                                                                                 >  2:158269/1‑149 (MQ=255)
                                                    aTACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGaca                                                                                                 >  1:417035/1‑149 (MQ=255)
                                                           ggCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGt                                                                                          <  1:373211/149‑1 (MQ=255)
                                                              aaCAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCGGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGtttt                                                                                       <  2:662572/149‑1 (MQ=255)
                                                              aaCAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCGGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGtttt                                                                                       <  2:132756/149‑1 (MQ=255)
                                                                   gCACAGCAGGTATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCCAGGATGTGCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGc                                                                                  <  1:184878/149‑1 (MQ=255)
                                                                    cacaGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGcc                                                                                 >  2:548819/1‑149 (MQ=255)
                                                                        ccagGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGCTCGTGGTGGCGTAGGCAACAACCCAGGTGGAGCCCATCATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTta                                                                             <  1:922933/148‑1 (MQ=255)
                                                                         cagGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTtat                                                                            >  2:999342/1‑149 (MQ=255)
                                                                            ggATGCGATAATCGGGGCGCCCATGGCTTCAAGCAGGTCAAGGATGTGCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGGCGCCTTTATAAc                                                                         <  2:594250/149‑1 (MQ=255)
                                                                            ggATGAGATAATCGGTGCGCCCATGGCTTCATTCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAAc                                                                         <  1:72474/149‑1 (MQ=255)
                                                                              aTGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTGAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCa                                                                       >  1:193507/1‑149 (MQ=255)
                                                                                gCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCaaa                                                                     <  2:762840/149‑1 (MQ=255)
                                                                                  gATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCTGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACt                                                                   >  2:474120/1‑149 (MQ=255)
                                                                                     aaTCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCa                                                                <  1:792844/149‑1 (MQ=255)
                                                                                       tCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAgg                                                              <  2:263889/149‑1 (MQ=255)
                                                                                        cGGAGCGCCCATGGCTTCAATCAGTTCAAGTATGTTCGGGTTCGCGTAGGCAACAACCCAGGTGGAGACCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGa                                                             <  2:495272/149‑1 (MQ=255)
                                                                                         ggTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGAc                                                            <  1:204329/149‑1 (MQ=255)
                                                                                          gTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGGTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGAcc                                                           <  1:544318/149‑1 (MQ=255)
                                                                                             cgcCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAgg                                                        >  2:882275/1‑149 (MQ=255)
                                                                                             cgcCCATGGCGTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAgg                                                        <  1:112978/149‑1 (MQ=255)
                                                                                               cccATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGcc                                                      >  1:407570/1‑149 (MQ=255)
                                                                                                ccATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACCCTGTAGTTTTGTCGCATTTATAACCACACTTCAGGACCAGGCCa                                                     >  2:804892/1‑149 (MQ=255)
                                                                                                ccATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAggata                                                     >  2:24215/1‑146 (MQ=255)
                                                                                                ccATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAggata                                                     <  1:24215/149‑4 (MQ=255)
                                                                                                    ggCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGTCAACAACCCAGGTGGAGACCATGATTAAGATCATGCTGATAGTGTTCAGTTTACCCCGCGCCACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCCTTcc                                                 >  2:968514/1‑148 (MQ=255)
                                                                                                    ggCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCAttaa                                                 >  1:474301/1‑147 (MQ=255)
                                                                                                     gCTTCAATCAGGTGAAGGATGTTCGGGGTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCATTTTACCCAGCGACACTTTAGTTTTTTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAg                                                <  2:488925/149‑1 (MQ=255)
                                                                                                     gCTTCAATCAGGTCAAGGATGTTGGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAg                                                <  2:446784/149‑1 (MQ=255)
                                                                                                       ttCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGTAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTAACCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGAc                                              >  1:675244/1‑149 (MQ=255)
                                                                                                         cAATCAGGTCAAGGTTGTTCGGGTTGGCGTAGGCTACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACTCTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCt                                            <  2:96515/149‑1 (MQ=255)
                                                                                                             cAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCa                                        <  1:809453/149‑1 (MQ=255)
                                                                                                               tgTCAAGGATGTTCGGGTTGGCGTATGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGc                                      <  1:599111/148‑1 (MQ=255)
                                                                                                                  cAAGGATGTTCGGGTTTGCGTAGGCAACAACCCAGGGGTAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGtt                                   <  1:162769/149‑1 (MQ=255)
                                                                                                                  cAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGtt                                   <  1:50071/149‑1 (MQ=255)
                                                                                                                  cAAGGATGTTCGGGGTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGtt                                   <  2:893527/149‑1 (MQ=255)
                                                                                                                   aaGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTGTTGTCGCCTTTATAACCAAACTTCAGGCCGAGGCCATTCAGACCTTCGAGCGTTc                                  >  2:356883/1‑149 (MQ=255)
                                                                                                                          ttCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGGTCCCAGATa                           >  1:532505/1‑149 (MQ=255)
                                                                                                                               gTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGGTCCCAGATAGTGAc                      >  1:400310/1‑149 (MQ=255)
                                                                                                                                     gTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGTTCCCAGATAGTGACCgaaga                <  1:158269/149‑1 (MQ=255)
                                                                                                                                          caacaaCCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGTTCCCAGATAGTGACCGAAGAAAGAt           <  1:169049/149‑1 (MQ=255)
                                                                                                                                          caacaaCCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTATAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGTTCCCAGATAGTGACCGAAGAAAGAt           <  2:186206/149‑1 (MQ=255)
                                                                                                                                            acaacCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGTTCCCAGATAGTGACCGAGGAAAGAttt         >  1:441410/1‑149 (MQ=255)
                                                                                                                                            acaacCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGTTCCCAGATAGTGACCGAAGAAAGAttt         <  2:345688/149‑1 (MQ=255)
                                                                                                                                            acaacCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGTTCCCAGATAGTGACCGAAGAAAGAttt         <  2:146853/149‑1 (MQ=255)
                                                                                                                                             aaacCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGTTCCCAGATAGTGACCGAAGAAAGATTTg        <  1:488270/148‑1 (MQ=255)
                                                                                                                                              aacCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGTTCCCAGATAGTGACCGAAGAAAGATTTGa       <  2:417035/149‑1 (MQ=255)
                                                                                                                                              aacCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGTTCCCAGATAGTGACCGAAGAAAGATTTGa       >  1:380163/1‑149 (MQ=255)
                                                                                                                                                   aGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGTTCCCAGATAGTGACCGAAGAAAGATTTGAAGATa  <  1:997523/149‑1 (MQ=255)
                                                                                                                                                   aGGTGGAGCCCATGATGAAGATCATGCTGATAGAGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGTTCCCAGATAGTGACAGAAGAAAGATTTGAAGATa  >  1:329232/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGTTATCCAGACGACCACGGTATTTCGCCAGAGACGGCGCTTTACGGATGGCATACATCGGCAACAGGCACAGCAGGGATGCGATAATCGGTGCGCCCATGGCTTCAATCAGGTCAAGGATGTTCGGGTTGGCGTAGGCAACAACCCAGGTGGAGCCCATGATGAAGATCATGCTGATAGTGTTCAGTTTACCCAGCGACACTTTAGTTTTGTCGCCTTTATAACCAAACTTCAGGACCAGGCCATTCAGACCTTCCAGCGTTCCCAGATAGTGACCGAAGAAAGATTTGAAGATA  >  NC_000913/3263745‑3264040

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑