Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 916,383 A→T 5.1% intergenic (‑336/‑90) aqpZ ← / → ybjD aquaporin Z/putative OLD family ATP‑dependent endonuclease; DUF2813 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913916,3830AT 5.1% 394.5 / 5.1 127intergenic (‑336/‑90)aqpZ/ybjDaquaporin Z/putative OLD family ATP‑dependent endonuclease; DUF2813 family protein
Reads supporting (aligned to +/- strand):  ref base A (54/66);  new base T (3/4);  total (57/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

TGACAGTGTAGTGCACGCAAAAGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGAT  >  NC_000913/916235‑916531
                                                                                                                                                    |                                                                                                                                                    
tGACAGTGTAGTGCACGCAAAAGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCa                                                                                                                                                      <  2:623117/149‑1 (MQ=255)
  aCAGTGTAGTGCACGCAAAAGATTTAATCCTTTAGGCGTAATATAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATg                                                                                                                                                    <  1:1046220/149‑1 (MQ=255)
   ctgtgtAGTGCACGCAAAAGTTTTAATCCTTTAGGAGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATATCGCGATTATTCGCTCATTTTTCAGACATTTGCCATGc                                                                                                                                                   <  1:298768/147‑1 (MQ=255)
   cAGTGTAGTGCACGCAAAAGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATATCGCGATTATTCGCTCATTTTTCAGACATTTGCCATGc                                                                                                                                                   <  1:856847/149‑1 (MQ=255)
   cAGTGTAGTGCACGCAAAAGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGc                                                                                                                                                   >  2:281592/1‑149 (MQ=255)
    aGTGTAGTGCACGCAAAAGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCt                                                                                                                                                  >  2:1213569/1‑149 (MQ=255)
    aGTGTAGTGCACGCAAAAGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCt                                                                                                                                                  <  2:1306489/149‑1 (MQ=255)
        tAGTGCACGCAAATGATTTAATGCTTTAGGCGTAATAAATAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTaaa                                                                                                                                              <  1:905373/149‑1 (MQ=255)
        tAGTGCACGCAAAAGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTaaa                                                                                                                                              <  1:1489275/149‑1 (MQ=255)
           tGCACGCAAAAGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAAtgt                                                                                                                                           <  2:1425212/149‑1 (MQ=255)
             cACGCAAAAGATTTAATCCTTTTTGTGTAATAATAAATAATTTTTGATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGa                                                                                                                                         <  2:134576/149‑1 (MQ=255)
             cACGCAAAAGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGa                                                                                                                                         <  1:1595066/149‑1 (MQ=255)
             cACGCAAAAGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTAGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGa                                                                                                                                         >  1:101036/1‑149 (MQ=255)
             cACGCAAAAGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAAAGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGa                                                                                                                                         >  2:549855/1‑149 (MQ=255)
              aCGCAAAAGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGAt                                                                                                                                        >  1:415758/1‑149 (MQ=255)
               cGCAAAAGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATg                                                                                                                                       >  2:222747/1‑149 (MQ=255)
                gCAAAAGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGAttt                                                                                                                                      >  2:1167045/1‑147 (MQ=255)
                  aaaaGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGtca                                                                                                                                    >  2:737830/1‑149 (MQ=255)
                   aaaGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCAAGCTAATTATTTGATTTAGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGtcat                                                                                                                                   >  2:972352/1‑149 (MQ=255)
                    aaGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATGATGCTAATTATTTGATTTTGTTTTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGtcatc                                                                                                                                  <  1:1437452/149‑1 (MQ=255)
                    aaGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGtcatc                                                                                                                                  <  2:467806/149‑1 (MQ=255)
                      gATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCAc                                                                                                                                <  2:1485462/149‑1 (MQ=255)
                        tttAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCTTGCTAATACGTGATGACATCAcat                                                                                                                              >  2:1306385/1‑147 (MQ=255)
                        tttAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGt                                                                                                                              >  1:827922/1‑149 (MQ=255)
                        tttAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGt                                                                                                                              >  2:673179/1‑149 (MQ=255)
                         ttAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTa                                                                                                                             >  2:1151061/1‑149 (MQ=255)
                               cTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGc                                                                                                                        >  2:404344/1‑148 (MQ=255)
                               cTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGc                                                                                                                        <  1:404344/148‑1 (MQ=255)
                                   aGGCGTAATTAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGc                                                                                                                   <  1:1589840/149‑1 (MQ=255)
                                       gTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGATGTAGTATAACGAGATTATTAGCTCATTTTTCAGACATTTGCCATGCTTAAATTTGATGTCATAACTTATTAGCAAGGCCttt                                                                                                               >  1:53864/1‑149 (MQ=255)
                                       gTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCttt                                                                                                               >  1:1202393/1‑149 (MQ=255)
                                        taataaAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTc                                                                                                              <  2:1340264/149‑1 (MQ=255)
                                         aataaAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTcc                                                                                                             <  1:1433141/149‑1 (MQ=255)
                                          ataaAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACTCGATTATTCGCTCATTTTTCAGACATTTGCCTTGCTAATACGTGATGACATCACATTTAAGCATGGcaaatgtc                                                                                                            <  2:56940/149‑8 (MQ=25)
                                          ataaAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCTTGCTAATACGTGATGACATCACATTTAAGCATGGcaaatgtc                                                                                                            <  1:1306385/149‑8 (MQ=37)
                                          ataaAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTccc                                                                                                            <  2:390839/149‑1 (MQ=255)
                                           taaAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCg                                                                                                           <  1:188326/149‑1 (MQ=255)
                                             aaaaaTAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGtt                                                                                                         >  1:1539460/1‑149 (MQ=255)
                                              aaaaTAATTTATCATGGTTATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATATCGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTta                                                                                                        <  2:463333/149‑1 (MQ=255)
                                              aaaaTAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTta                                                                                                        >  2:1513225/1‑149 (MQ=255)
                                               aaaTAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTtat                                                                                                       >  2:565411/1‑149 (MQ=255)
                                                aataatTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTtata                                                                                                      >  2:468061/1‑149 (MQ=255)
                                                  taatTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACt                                                                                                    <  1:1408093/149‑1 (MQ=255)
                                                    atTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGc                                                                                                  <  2:702545/149‑1 (MQ=255)
                                                     tttATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGcc                                                                                                 >  1:1597161/1‑149 (MQ=255)
                                                     tttATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGcc                                                                                                 >  2:948463/1‑149 (MQ=255)
                                                            tGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTaa                                                                                          >  2:43740/1‑149 (MQ=255)
                                                             gatAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGAGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTaaa                                                                                         <  2:593993/147‑1 (MQ=255)
                                                             gCTAATTATTTTAATTTGTTGTTTTTGCAGACTTTTCAGCAAGAGGGAGTAAAACGCGATTATTCGCTCAATTTTCACACATTTGCCATGGTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTaaa                                                                                         <  2:263268/149‑1 (MQ=255)
                                                             gCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTaaa                                                                                         >  2:890120/1‑149 (MQ=255)
                                                             gCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTaaa                                                                                         <  1:396783/149‑1 (MQ=255)
                                                              cTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAg                                                                                        <  1:803092/149‑1 (MQ=255)
                                                              cTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAg                                                                                        <  2:1392791/149‑1 (MQ=255)
                                                              cTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAg                                                                                        >  2:728492/1‑149 (MQ=255)
                                                               tAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAgg                                                                                       >  2:1446035/1‑149 (MQ=255)
                                                                 attattTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGc                                                                                                                   >  1:102810/1‑119 (MQ=255)
                                                                 attattTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGc                                                                                                                   <  2:102810/119‑1 (MQ=255)
                                                                      ttGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTc                                                                                <  1:856189/149‑1 (MQ=255)
                                                                      ttGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTc                                                                                >  2:838022/1‑149 (MQ=255)
                                                                      ttGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTc                                                                                <  1:1521211/149‑1 (MQ=255)
                                                                                  tttttGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTATTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTgg                                                                    >  2:118563/1‑149 (MQ=255)
                                                                                   ttttGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGa                                                                   <  2:415758/149‑1 (MQ=255)
                                                                                    tttGTAGACTTATCAGCATGATGTTGTATAACGCTATTATTCGCTCATTTTTCAGTCTTTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCGTTTCCCGTTGTACTGCCAGCGTAAAGGATAAGTCACATATTTCTggtg                                                                  <  1:590793/149‑3 (MQ=255)
                                                                                    tttGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAg                                                                  <  1:114972/149‑1 (MQ=255)
                                                                                        cAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGa                                                              >  2:812499/1‑149 (MQ=255)
                                                                                          gACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGata                                                            <  2:181497/149‑1 (MQ=255)
                                                                                          gACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGata                                                            <  1:707322/149‑1 (MQ=255)
                                                                                             ttATCAGCAAGAGGTAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTTCAATGCTTAAATGTGATGTCATGTCGTATTAGCAATGCCTTTCCCGTTATACTGCCAGCGTAATGGATAAGTCACTTATTTCTGGAGGGGATATGa                                                         <  1:729846/149‑1 (MQ=255)
                                                                                               aTCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGAtt                                                       <  1:1213569/149‑1 (MQ=255)
                                                                                               aTCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGAtt                                                       >  1:1561920/1‑149 (MQ=255)
                                                                                                tCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTc                                                      >  2:1221934/1‑149 (MQ=255)
                                                                                                  aGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCtt                                                    >  1:977082/1‑149 (MQ=255)
                                                                                                    cAAGAGTGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGa                                                  <  1:589736/149‑1 (MQ=255)
                                                                                                    cAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGa                                                  >  1:502150/1‑149 (MQ=255)
                                                                                                     aaGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAg                                                 >  1:1053021/1‑149 (MQ=255)
                                                                                                        agGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAgcgc                                              <  1:281592/149‑1 (MQ=255)
                                                                                                          ggAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGt                                            >  2:782758/1‑149 (MQ=255)
                                                                                                            aGTATAACGCGATTGTTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTg                                          <  1:972674/149‑1 (MQ=255)
                                                                                                                taACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAAt                                      >  1:1306408/1‑149 (MQ=255)
                                                                                                                 aaCGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAAtt                                     <  1:942380/149‑1 (MQ=255)
                                                                                                                   cgcgATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATtgt                                   <  1:1167045/149‑1 (MQ=255)
                                                                                                                   cgcgATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATtgt                                   >  1:1508133/1‑149 (MQ=255)
                                                                                                                   cgcgATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATtgt                                   <  2:1482326/149‑1 (MQ=255)
                                                                                                                     gcgatattCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTgg                                 >  1:1555903/5‑149 (MQ=255)
                                                                                                                     cgATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTgg                                 <  1:737830/149‑1 (MQ=255)
                                                                                                                     cgATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTgg                                 <  2:1539460/149‑1 (MQ=255)
                                                                                                                     cgATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTgg                                 <  1:948463/149‑1 (MQ=255)
                                                                                                                              gCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTcgcg                        >  2:1203132/1‑149 (MQ=255)
                                                                                                                              gCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGACTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTcgcg                        <  2:1201218/149‑1 (MQ=255)
                                                                                                                                 cATTTTTCTGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTa                     <  2:53864/149‑1 (MQ=255)
                                                                                                                                  aTTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTAt                    >  2:174779/1‑149 (MQ=255)
                                                                                                                                      ttCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGAGGGTTTTCGCGGTATCAAc                >  1:841583/1‑149 (MQ=255)
                                                                                                                                        cAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCg              >  2:708268/1‑149 (MQ=255)
                                                                                                                                         aGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGt             <  1:118563/149‑1 (MQ=255)
                                                                                                                                           ggaaaggCCTTGCTAATACGTGATGACATCACATTTAAGCATGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGttt           >  1:943938/7‑149 (MQ=37)
                                                                                                                                              tttGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTc        <  1:461409/149‑1 (MQ=255)
                                                                                                                                               aggCCTTGCTAATACGTGATGACATCACATTTAAGCATGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGttt           >  1:186481/3‑145 (MQ=37)
                                                                                                                                               aggCCTTGCTAATACGTGATGACATCACATTTAAGCATGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGttt           <  2:186481/143‑1 (MQ=37)
                                                                                                                                                tggcATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGt      <  1:425452/146‑1 (MQ=255)
                                                                                                                                                 gCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGtt     <  1:33667/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGACAGTGTAGTGCACGCAAAAGATTTAATCCTTTAGGCGTAATAAAAAATAATTTATCATGCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGAT  >  NC_000913/916235‑916531

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