Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,136,629 C→A 7.1% intergenic (‑258/+34) pitB ← / ← gss phosphate transporter/glutathionylspermidine amidase and glutathionylspermidine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,136,6290CA 7.1% 487.5 / 18.2 172intergenic (‑258/+34)pitB/gssphosphate transporter/glutathionylspermidine amidase and glutathionylspermidine synthetase
Reads supporting (aligned to +/- strand):  ref base C (84/74);  new base A (6/8);  total (90/82)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 182 total aligned reads displayed.

TGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATACC  >  NC_000913/3136481‑3136777
                                                                                                                                                    |                                                                                                                                                    
tgtgTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTAttct                                                                                                                                                     <  1:309318/150‑1 (MQ=255)
tgtgTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTAttct                                                                                                                                                     >  1:1553864/1‑150 (MQ=255)
  tgtATCAGACAATAAGGATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTAttcttc                                                                                                                                                   <  1:822776/150‑1 (MQ=255)
  tgtATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTAttcttc                                                                                                                                                   >  2:672288/1‑150 (MQ=255)
   gtATCAGACAATAAGAATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCACCGTACTGATTATTCTTCa                                                                                                                                                  >  1:842949/1‑150 (MQ=255)
     aTCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTATTAATa                                                                                                                                                >  1:1018414/1‑150 (MQ=255)
      tCATACAATAATCATTCATTAAAAAACTTTAGAAATAATCAAAACAAATGTATTCCGCAGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATaa                                                                                                                                               <  1:43457/150‑1 (MQ=255)
      tCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATaa                                                                                                                                               <  2:1405070/150‑1 (MQ=255)
        aGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATc                                                                                                                                             >  1:298659/1‑150 (MQ=255)
        aGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATc                                                                                                                                             <  1:1297832/150‑1 (MQ=255)
           cAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTATTCATAATAAGt                                                                                                                                          >  2:197859/1‑150 (MQ=255)
            aaTAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTAATGATTATTCTTCATAATCAGTa                                                                                                                                         >  2:814461/1‑150 (MQ=255)
             aTAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTAc                                                                                                                                        <  1:631955/150‑1 (MQ=255)
                aGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACggg                                                                                                                                     <  2:1219248/150‑1 (MQ=255)
                aGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTACggg                                                                                                                                     <  2:873979/150‑1 (MQ=255)
                 gCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACgggg                                                                                                                                    <  2:1508768/150‑1 (MQ=255)
                     tcattcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGTGGTTGc                                                                                                                                >  2:179702/1‑150 (MQ=255)
                      cattcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCt                                                                                                                               >  1:869566/1‑150 (MQ=255)
                      cattcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCt                                                                                                                               <  1:1625666/150‑1 (MQ=255)
                        ttcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTac                                                                                                                             >  2:41991/1‑150 (MQ=255)
                        ttcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTac                                                                                                                             >  2:532905/1‑150 (MQ=255)
                        ttcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTACggggattaatc                                                                                                                             <  1:601733/150‑8 (MQ=38)
                          caAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTacaa                                                                                                                           <  1:423991/150‑1 (MQ=255)
                           aaaaCCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAAt                                                                                                                                              <  2:329064/130‑1 (MQ=255)
                           aaaaCCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAAt                                                                                                                                              >  1:329064/1‑130 (MQ=255)
                            aaaCCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGTGGTTGCTacaaca                                                                                                                         >  1:669773/1‑150 (MQ=255)
                            aaaCCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTacaaca                                                                                                                         >  1:1004720/1‑150 (MQ=255)
                                cTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTa                                                                                                                     <  1:314340/150‑1 (MQ=255)
                                  ttAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTAtt                                                                                                                   <  2:1133770/150‑1 (MQ=255)
                                  ttAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTAtt                                                                                                                   <  1:856833/150‑1 (MQ=255)
                                   tAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTa                                                                                                                  >  2:349023/1‑150 (MQ=255)
                                    aGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAATACTGTATTTATACAACATGTATTAc                                                                                                                 >  1:100927/1‑150 (MQ=255)
                                       aataatCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACttt                                                                                                              <  1:789945/150‑1 (MQ=255)
                                         taatCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACttttt                                                                                                            <  1:1487636/150‑1 (MQ=255)
                                         taatCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACttttt                                                                                                            <  2:857791/150‑1 (MQ=255)
                                          aatCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTc                                                                                                           <  2:1202052/150‑1 (MQ=255)
                                           atCAAAACACATGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTca                                                                                                          <  2:1701360/150‑1 (MQ=255)
                                           atCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTca                                                                                                          <  2:806073/150‑1 (MQ=255)
                                           atCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTca                                                                                                          >  1:975040/1‑150 (MQ=255)
                                           atCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTca                                                                                                          <  1:387543/150‑1 (MQ=255)
                                           atCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTca                                                                                                          <  1:1561508/150‑1 (MQ=255)
                                            tCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTcac                                                                                                         <  1:525288/150‑1 (MQ=255)
                                             cAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTcacc                                                                                                        <  1:1601499/150‑1 (MQ=255)
                                              aaaaCACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTTATTAATCCCCTTACTGATTATTATTCATAATCAGTACGGGGTTGCTACAACATTTATTAATTTTTcacaa                                                                                                       >  1:1234041/1‑148 (MQ=255)
                                               aaaCACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTcaccac                                                                                                      <  1:1276421/150‑1 (MQ=255)
                                                aaCACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACcaca                                                                                                     <  1:720380/150‑1 (MQ=255)
                                                aaCACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACcaca                                                                                                     <  2:1553864/150‑1 (MQ=255)
                                                  cacaGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAAt                                                                                                   >  1:536426/1‑150 (MQ=255)
                                                           tCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTc                                                                                          <  2:298659/150‑1 (MQ=255)
                                                           tCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCAACACAATAAACTGtta                                                                                          >  2:1077221/1‑149 (MQ=255)
                                                            ccACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCa                                                                                         >  2:720744/1‑150 (MQ=255)
                                                             cACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCaa                                                                                        >  1:1380844/1‑150 (MQ=255)
                                                              aCTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAat                                                                                       <  2:751564/150‑1 (MQ=255)
                                                               cTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAata                                                                                      <  1:1344134/150‑1 (MQ=255)
                                                                tGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAatat                                                                                     <  2:1004720/150‑1 (MQ=255)
                                                                 ggTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATc                                                                                    >  2:1467805/1‑150 (MQ=255)
                                                                   ttAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATTTCAc                                                                                  >  1:644619/1‑150 (MQ=255)
                                                                   ttAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCAc                                                                                  <  2:536426/150‑1 (MQ=255)
                                                                   ttAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCAc                                                                                  >  2:894307/1‑150 (MQ=255)
                                                                             cacaTATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCtttttt                                                                        >  2:1578608/1‑150 (MQ=255)
                                                                             aacaTATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCtttttt                                                                        <  1:1178671/149‑1 (MQ=255)
                                                                              acaTATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTtg                                                                       >  1:617671/1‑150 (MQ=255)
                                                                               caTATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTtga                                                                      <  2:1225201/150‑1 (MQ=255)
                                                                               caTATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTtga                                                                      >  2:203778/1‑150 (MQ=255)
                                                                               caTATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTtga                                                                      >  2:677460/1‑150 (MQ=255)
                                                                                        tgaaTAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTg                                                             <  2:1380844/150‑1 (MQ=255)
                                                                                         gaaTAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGa                                                            >  2:810778/1‑150 (MQ=255)
                                                                                          aataaataTTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGAt                                                           <  1:412060/150‑1 (MQ=255)
                                                                                          aataaataTTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGAt                                                           >  1:1539020/1‑150 (MQ=255)
                                                                                               atatTGGAATATATAAATATTGAATATTTTGATTAATCCCAGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATtcatt                                                      <  1:179702/150‑2 (MQ=255)
                                                                                                tatTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACAAGTGATtcatca                                                     >  1:1001619/1‑150 (MQ=255)
                                                                                                  ttGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCAcc                                                   <  1:1400245/150‑1 (MQ=255)
                                                                                                  ttGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCAcc                                                   >  1:979366/1‑150 (MQ=255)
                                                                                                   tGGAATATATAAATATTGAATATTTTGTTTAATACCAGTACTTATTATTCTTCATAATCAGTACGGGGTTGCTAAAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATTACCAGTGATTCATCACcg                                                  <  2:591108/150‑1 (MQ=255)
                                                                                                     gAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACcgc                                                 <  1:254144/149‑1 (MQ=255)
                                                                                                     gAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACcgc                                                 <  1:1116876/149‑1 (MQ=255)
                                                                                                     gAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACcgc                                                 >  2:254144/1‑149 (MQ=255)
                                                                                                     gAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACcgc                                                 >  2:1116876/1‑149 (MQ=255)
                                                                                                         atataAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCaaa                                            >  2:1468853/1‑150 (MQ=255)
                                                                                                          tataAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAac                                           >  2:779916/1‑150 (MQ=255)
                                                                                                           ataAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAaca                                          <  2:1676714/150‑1 (MQ=255)
                                                                                                           ataAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAaca                                          <  2:1234041/150‑1 (MQ=255)
                                                                                                               atatTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACAACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTc                                      >  2:1324224/1‑150 (MQ=255)
                                                                                                                tatTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTcc                                     <  1:1685383/150‑1 (MQ=255)
                                                                                                                    gAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCAcc                                 >  2:571852/1‑150 (MQ=255)
                                                                                                                     aaTATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCg                                <  2:965967/150‑1 (MQ=255)
                                                                                                                     aaTATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCATTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCg                                >  1:74056/1‑150 (MQ=255)
                                                                                                                      atatTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGt                               <  1:202519/150‑1 (MQ=255)
                                                                                                                            tGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTgc                         >  2:1690013/1‑150 (MQ=255)
                                                                                                                              aTTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCCCCACCATTAACGGTTCATTATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTgccg                       >  1:1496069/1‑150 (MQ=255)
                                                                                                                                tAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTgccgcc                     <  2:62078/150‑1 (MQ=255)
                                                                                                                                 aaTCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCa                    >  2:391909/1‑150 (MQ=255)
                                                                                                                                 aaTCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCa                    >  1:999450/1‑150 (MQ=255)
                                                                                                                                       cGTACTGATTATGAAGAATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTg              >  1:640992/1‑150 (MQ=255)
                                                                                                                                          aCTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAg           >  2:1254750/1‑150 (MQ=255)
                                                                                                                                          aCTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAg           >  1:906694/1‑150 (MQ=255)
                                                                                                                                            tGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGt         >  1:1096970/1‑150 (MQ=255)
                                                                                                                                             gATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTa        <  1:1583728/150‑1 (MQ=255)
                                                                                                                                                 attCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATa    <  2:1023346/150‑1 (MQ=255)
                                                                                                                                                   tcttcATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATAcc  >  1:964229/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATACC  >  NC_000913/3136481‑3136777

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