Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,493,194 G→C 7.1% S242T (AGT→ACT)  trg → methyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,493,1940GC 7.1% 497.8 / 22.5 140S242T (AGT→ACT) trgmethyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor
Reads supporting (aligned to +/- strand):  ref base G (64/66);  new base C (5/5);  total (69/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 140 total aligned reads displayed.

CCATCAGCGCACCCGCCTGGGTGGGATGTTCATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACC  >  NC_000913/1493048‑1493336
                                                                                                                                                  |                                                                                                                                              
ccATCAGCGCACCCGCCTGGGTGGGATGTTCATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTg                                                                                                                                              >  4:844577/1‑149 (MQ=255)
  aTCAGCGCACCCGCCTGGGTGGGATGTTCATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCACTGGc                                                                                                                                            >  3:1150755/1‑149 (MQ=255)
    cAGCGCACCCGCCTGGGTGGGATGTTCATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGa                                                                                                                                          >  1:1219092/1‑149 (MQ=255)
     aGCGCACCCGCCTGGGTGGGATGTTCATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGAt                                                                                                                                         <  4:230664/149‑1 (MQ=255)
       cgcACCCGCCTGGGTGGGATGTTCATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCt                                                                                                                                       <  2:1160738/149‑1 (MQ=255)
       cgcACCCGCCTGGGTGGGATGTTCATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCACTGGCGAttt                                                                                                                                       <  4:1150755/149‑3 (MQ=255)
          aCCCGCCTGGGTGGGATGTTCATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCACCATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGAc                                                                                                                                    <  2:693131/149‑1 (MQ=255)
                cTGGGTGGGATGTTCATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGatgaa                                                                                                                              <  2:281372/149‑1 (MQ=255)
                    gtgggATGTTCATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAatgat                                                                                                                          >  4:948251/1‑149 (MQ=255)
                    gtgggATGTTCATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAatgat                                                                                                                          <  1:665051/149‑1 (MQ=255)
                    gtgggATGTTCATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAatgat                                                                                                                          <  3:948251/149‑1 (MQ=255)
                          tGTTCATGATTGGCGCGTTTGTGCTTGCCCTGGTCACGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACcg                                                                                                                    <  2:996589/149‑1 (MQ=255)
                              cATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGa                                                                                                                                     >  3:1219861/1‑128 (MQ=255)
                              cATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGa                                                                                                                                     <  2:118125/128‑1 (MQ=255)
                              cATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGa                                                                                                                                     <  4:1219861/128‑1 (MQ=255)
                              cATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGa                                                                                                                                     >  1:118125/1‑128 (MQ=255)
                              cATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACcggcgg                                                                                                                >  3:883162/1‑149 (MQ=255)
                                 gATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTc                                                                                                             >  2:644941/1‑149 (MQ=255)
                                      gcgcGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAAt                                                                                                        >  3:1072852/1‑149 (MQ=255)
                                      gcgcGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAAt                                                                                                        >  4:96061/1‑149 (MQ=255)
                                      gcgcGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAAt                                                                                                        >  2:1288743/1‑149 (MQ=255)
                                        gcgTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGa                                                                                                      >  3:1313084/1‑149 (MQ=255)
                                          gTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGaaa                                                                                                    <  1:940106/149‑1 (MQ=255)
                                           tttGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAAt                                                                                                   <  3:663059/149‑1 (MQ=255)
                                                gCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTc                                                                                              <  4:1072852/149‑1 (MQ=255)
                                                 cTTGCCCTGGTCATGACGCGGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGatg                                                                                                                                <  1:250955/114‑1 (MQ=255)
                                                      ccTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTaa                                                                                        <  1:1288743/149‑1 (MQ=255)
                                                      ccTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTaa                                                                                        <  1:897621/149‑1 (MQ=255)
                                                             aTGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCa                                                                                 >  2:40036/1‑149 (MQ=255)
                                                               gACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCACTGGCGAtttt                                                                                                                                      >  4:227995/1‑91 (MQ=255)
                                                               gACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCACTGGCGAtttt                                                                                                                                      <  3:227995/94‑4 (MQ=255)
                                                                aCGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCAttt                                                                              >  4:562579/1‑149 (MQ=255)
                                                                aCGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCAttt                                                                              >  2:656811/1‑149 (MQ=255)
                                                                  gCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTAc                                                                            <  1:40036/149‑1 (MQ=255)
                                                                   cTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACa                                                                           <  3:96061/149‑1 (MQ=255)
                                                                    tGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAg                                                                          >  1:643370/1‑149 (MQ=255)
                                                                    tGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAg                                                                          >  2:822482/1‑149 (MQ=255)
                                                                     gATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGa                                                                                                                                     <  3:323954/89‑1 (MQ=255)
                                                                     gATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGa                                                                                                                                     >  4:323954/1‑89 (MQ=255)
                                                                      aTAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCa                                                                        >  2:766031/1‑149 (MQ=255)
                                                                        aaCATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATc                                                                                                                                        >  1:1238377/1‑83 (MQ=255)
                                                                        aaCATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATc                                                                                                                                        <  2:1238377/83‑1 (MQ=255)
                                                                         aCATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAAt                                                                     <  1:822482/149‑1 (MQ=255)
                                                                              tATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgcagc                                                                <  4:883162/149‑1 (MQ=255)
                                                                                       aCGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTgg                                                       <  1:180550/149‑1 (MQ=255)
                                                                                         gTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAAc                                                                                                                      <  4:1428442/84‑1 (MQ=255)
                                                                                         gTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAAc                                                                                                                      >  3:1428442/1‑84 (MQ=255)
                                                                                         gTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCtt                                                                                          >  3:920036/1‑112 (MQ=255)
                                                                                         gTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTgggg                                                     >  2:342629/1‑149 (MQ=255)
                                                                                             gATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAgtcgtc                                                                                  >  4:617103/1‑116 (MQ=255)
                                                                                              aTCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGatgaa                                                                                                                              <  3:815652/71‑1 (MQ=255)
                                                                                              aTCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTc                                                                                              >  2:461019/1‑103 (MQ=255)
                                                                                              aTCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACa                                                                           <  3:247396/122‑1 (MQ=255)
                                                                                              aTCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGAc                                                >  3:812498/1‑149 (MQ=255)
                                                                                               tCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACa                                               >  4:110596/1‑149 (MQ=255)
                                                                                                 gTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCAc                                                          >  3:548573/1‑136 (MQ=255)
                                                                                                 gTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGt                                             >  1:1116604/1‑149 (MQ=255)
                                                                                                  tCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGAt                                                   <  4:725699/142‑1 (MQ=255)
                                                                                                  tCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGAt                                                   >  3:481275/1‑142 (MQ=255)
                                                                                                  tCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGAt                                                   <  4:481275/142‑1 (MQ=255)
                                                                                                  tCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGAt                                                   >  3:725699/1‑142 (MQ=255)
                                                                                                    aTTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCaaa                                                                      >  1:446521/1‑121 (MQ=255)
                                                                                                    aTTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAg                                           <  2:324006/148‑1 (MQ=255)
                                                                                                    aTTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAg                                           >  1:324006/1‑148 (MQ=255)
                                                                                                     ttCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCg                                                                                                 >  2:207788/1‑93 (MQ=255)
                                                                                                      tCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGa                                        <  4:1313084/149‑1 (MQ=255)
                                                                                                            aCTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTc                                  >  2:627319/1‑149 (MQ=255)
                                                                                                              tGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGa                                        >  3:1345839/1‑141 (MQ=255)
                                                                                                              tGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGa                                        <  4:1345839/141‑1 (MQ=255)
                                                                                                              tGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGa                                >  1:1370270/1‑149 (MQ=255)
                                                                                                              tGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGa                                <  1:584392/149‑1 (MQ=255)
                                                                                                                cAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAgg                                          >  3:1303087/1‑137 (MQ=255)
                                                                                                                cAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACa                              <  1:656811/149‑1 (MQ=255)
                                                                                                                cAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACa                              <  3:562579/149‑1 (MQ=255)
                                                                                                                         gCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAg                                              <  2:417464/124‑1 (MQ=255)
                                                                                                                         gCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAg                                              >  1:417464/1‑124 (MQ=255)
                                                                                                                             aaCGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTa                                            >  4:861400/1‑122 (MQ=255)
                                                                                                                             aaCGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGCGGATGACAGTa                                            <  3:861400/122‑1 (MQ=255)
                                                                                                                              aCGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAggg                           <  3:327841/138‑1 (MQ=255)
                                                                                                                              aCGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAggg                           >  4:327841/1‑138 (MQ=255)
                                                                                                                               cGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGAc                               >  1:559107/1‑133 (MQ=255)
                                                                                                                               cGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGAc                               <  2:559107/133‑1 (MQ=255)
                                                                                                                                 gATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgcag                                                                 <  2:640080/97‑1 (MQ=255)
                                                                                                                                 gATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgcag                                                                 >  1:640080/1‑97 (MQ=255)
                                                                                                                                 gATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGAttt             <  1:627319/149‑1 (MQ=255)
                                                                                                                                  aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACcggc                                                                                                                  <  1:1460609/47‑1 (MQ=255)
                                                                                                                                  aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACcggc                                                                                                                  >  2:1460609/1‑47 (MQ=255)
                                                                                                                                  aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAAt                                                                     <  3:633512/92‑1 (MQ=255)
                                                                                                                                  aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAAt                                                                     >  4:633512/1‑92 (MQ=255)
                                                                                                                                  aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGaa                    <  4:30488/141‑1 (MQ=255)
                                                                                                                                  aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGaa                    >  3:30488/1‑141 (MQ=255)
                                                                                                                                  aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTa            >  2:478757/1‑149 (MQ=255)
                                                                                                                                    tGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATc          <  4:812498/149‑1 (MQ=255)
                                                                                                                                       tcagaTCGCCACTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGAt                                                   >  4:1359267/5‑105 (MQ=255)
                                                                                                                                       tcagaTCGCCACTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGAt                                                   <  3:1359267/101‑1 (MQ=255)
                                                                                                                                          gaTCGCCACTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCa    >  4:1391712/2‑149 (MQ=255)
                                                                                                                                           aTCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCAc   >  2:197092/1‑149 (MQ=255)
                                                                                                                                           aTCGCCACTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGgcgc                        >  4:386759/1‑128 (MQ=255)
                                                                                                                                            tCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCAcc  >  2:942564/1‑149 (MQ=255)
                                                                                                                                               ccAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCAc                                                          <  4:1275233/90‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CCATCAGCGCACCCGCCTGGGTGGGATGTTCATGATTGGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACC  >  NC_000913/1493048‑1493336

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