Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 894,952 T→A 9.2% intergenic (+56/‑39) potF → / → potG putrescine ABC transporter periplasmic binding protein/putrescine ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913894,9520TA 9.2% 348.1 / 22.7 109intergenic (+56/‑39)potF/potGputrescine ABC transporter periplasmic binding protein/putrescine ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (47/52);  new base A (4/6);  total (51/58)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

Only 100 of 109 total aligned reads displayed.

CGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGAT  >  NC_000913/894807‑895098
                                                                                                                                                 |                                                                                                                                                   
cGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATg                                                                                                                                                  <  4:678485/149‑1 (MQ=255)
   gCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTg                                                                                                                                               <  3:251454/149‑1 (MQ=255)
      aaGCTGTTCACTCTGAAAGTGCAGGAGCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAAGAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTgcg‑c                                                                                                                                           <  2:683938/149‑1 (MQ=255)
       aGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATAATGGTgcgaa                                                                                                                                           <  1:114359/149‑3 (MQ=255)
       aGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATAATGGTgcgaa                                                                                                                                           <  3:161321/149‑3 (MQ=255)
       aGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATAATGGTgcgaa                                                                                                                                           <  4:265755/149‑3 (MQ=255)
       aGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATAATGGTgcgaa                                                                                                                                           <  1:342385/149‑3 (MQ=255)
        gCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATAATGGTgcgaat                                                                                                                                          <  2:687812/149‑4 (MQ=255)
            ttCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCAGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGcac                                                                                                                                     >  2:786212/1‑149 (MQ=255)
                 tctGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACAtt                                                                                                                                >  1:241555/1‑149 (MQ=255)
                 tctGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACAtt                                                                                                                                >  2:513244/1‑149 (MQ=255)
                 tctGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACAtt                                                                                                                                >  2:367149/1‑149 (MQ=255)
                 tctGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACAtt                                                                                                                                <  4:30990/149‑1 (MQ=255)
                 tctGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTACACACAtt                                                                                                                                >  2:132804/1‑149 (MQ=255)
                   tGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCa                                                                                                                              >  4:315031/1‑149 (MQ=255)
                     aaaGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAAt                                                                                                                            <  1:184590/149‑1 (MQ=255)
                     aaaGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAAt                                                                                                                            <  3:865961/149‑1 (MQ=255)
                     aaaGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAAt                                                                                                                            <  1:716645/149‑1 (MQ=255)
                         tGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCg                                                                                                                        >  3:531127/1‑149 (MQ=255)
                           cAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGga                                                                                                                      <  2:47401/149‑1 (MQ=255)
                            aGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCGCCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCggtg                                                                                                                     >  4:758011/1‑147 (MQ=255)
                            aGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGgag                                                                                                                     <  4:345701/149‑1 (MQ=255)
                             ggATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGgaga                                                                                                                    <  1:531737/149‑1 (MQ=255)
                             ggATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGgaga                                                                                                                    <  3:99361/149‑1 (MQ=255)
                                   gAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGcc                                                                                                              >  4:859816/1‑149 (MQ=255)
                                   gAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGcc                                                                                                              >  2:537272/1‑149 (MQ=255)
                                   gAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGcc                                                                                                              >  2:103271/1‑149 (MQ=255)
                                   gAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGcc                                                                                                              >  1:7005/1‑149 (MQ=255)
                                   gAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGcc                                                                                                              >  3:522251/1‑149 (MQ=255)
                                       aTCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATAATGGTgcga                                                                                                                                            <  1:683169/116‑2 (MQ=255)
                                       aTCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATAATGGTgcga                                                                                                                                            >  2:683169/1‑115 (MQ=255)
                                        tCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGaa                                                                                                         <  1:808739/149‑1 (MQ=255)
                                        tCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGaa                                                                                                         <  2:372395/149‑1 (MQ=255)
                                         cGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAAt                                                                                                        <  4:852280/149‑1 (MQ=255)
                                         cGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAAt                                                                                                        <  2:263700/149‑1 (MQ=255)
                                                    gcACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATccc                                                                                             <  1:305012/149‑1 (MQ=255)
                                                      aCCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATccccc                                                                                           >  1:809396/1‑147 (MQ=255)
                                                       cccGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCg                                                                                          >  1:211214/1‑149 (MQ=255)
                                                              gTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACAtt                                                                                                                                >  2:3322/1‑104 (MQ=255)
                                                              gTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACAtt                                                                                                                                <  1:3322/104‑1 (MQ=255)
                                                                ggACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGc                                                                                 >  2:486924/1‑149 (MQ=255)
                                                                  agcAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTGGATGGCGGAGGGGCGCACAACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAAGGCCGGAGAGCAGCCGTGAATAAAGCTATCCCTCGCCCACAGGCgg                                                                               >  4:353610/3‑148 (MQ=255)
                                                                    cAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCTGCGGGCGagg                                                                             <  1:552813/149‑3 (MQ=255)
                                                                     aaaGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGaaaa                                                                            >  1:671335/1‑149 (MQ=255)
                                                                          gAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGCCGAAAACCCGt                                                                       <  3:500922/149‑1 (MQ=255)
                                                                               gCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGgc                                                                  <  3:315031/149‑1 (MQ=255)
                                                                                cGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGgcg                                                                 >  2:237519/1‑149 (MQ=255)
                                                                                  gAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCt                                                               >  3:222153/1‑149 (MQ=255)
                                                                                  gAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCt                                                               <  1:786212/149‑1 (MQ=255)
                                                                                    aaaTAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGa                                                             <  1:651479/149‑1 (MQ=255)
                                                                                     aataatCCGCATTCGTAGATGCCGGAGGAGCGCACCACACCCGCTTGCAATTCTCACCATTATGGTGCG‑CTTGCACACATTCAATGCCAGAGTGCATCCGTGAATGACGCTATCCCGCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAc                                                            <  2:139754/149‑1 (MQ=255)
                                                                                     aataatCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAc                                                            >  4:254622/1‑149 (MQ=255)
                                                                                     aataatCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAc                                                            >  1:288014/1‑149 (MQ=255)
                                                                                         atCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAcgccg                                                        >  4:2803/1‑149 (MQ=255)
                                                                                          tCCGCAGTCGTGGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAcgccgc                                                       >  4:82139/1‑149 (MQ=255)
                                                                                             gCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTAt                                                    >  4:718822/1‑149 (MQ=255)
                                                                                               aGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTa                                                  <  2:241555/149‑1 (MQ=255)
                                                                                               aGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTa                                                  <  1:523243/149‑1 (MQ=255)
                                                                                               aGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTa                                                  <  3:27795/149‑1 (MQ=255)
                                                                                                    tAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAAt                                             <  4:567568/149‑1 (MQ=255)
                                                                                                      gATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCAGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATcc                                           <  4:635400/149‑1 (MQ=255)
                                                                                                          ccGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCaa                                       <  3:22378/149‑1 (MQ=255)
                                                                                                            ggAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAAcc                                     <  2:518630/149‑1 (MQ=255)
                                                                                                             gAGGGGCGCACCACACCCGCCGGCAATTCGCATCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGGTATCCCTCGCCCGCAGGCGAAAACCCCTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCt                                    <  3:742834/149‑1 (MQ=255)
                                                                                                             gAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCt                                    >  3:379728/1‑149 (MQ=255)
                                                                                                                gggCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGACATCCGCAACCTGAc                                 <  2:809396/149‑1 (MQ=255)
                                                                                                                gggCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGAc                                 <  1:589739/149‑1 (MQ=255)
                                                                                                                gggCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTACACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGAc                                 <  1:132804/149‑1 (MQ=255)
                                                                                                                  gtgcgccacACGCCCCGGCAAGTCGCACCATTATGGGCCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCa                               >  2:352133/6‑149 (MQ=255)
                                                                                                                  gcgcACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCg                                                                        >  3:172716/1‑108 (MQ=255)
                                                                                                                  gcgcACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCg                                                                        <  4:172716/108‑1 (MQ=255)
                                                                                                                   cgcACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCaa                              >  2:300611/1‑149 (MQ=255)
                                                                                                                   cgcACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCaa                              >  2:379587/1‑149 (MQ=255)
                                                                                                                   cgcACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCaa                              <  2:80602/149‑1 (MQ=255)
                                                                                                                     caccacACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAAt                            <  1:300611/149‑1 (MQ=255)
                                                                                                                     caccacACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAAt                            <  1:379587/149‑1 (MQ=255)
                                                                                                                      accacACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGCAATCCGCCACCTGACCAAATc                           >  1:660322/1‑149 (MQ=255)
                                                                                                                      accacACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATc                           >  1:437613/1‑149 (MQ=255)
                                                                                                                      accacACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATc                           <  2:205369/149‑1 (MQ=255)
                                                                                                                           acCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACg                      >  4:701243/1‑149 (MQ=255)
                                                                                                                            cccGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGa                     <  1:367149/149‑1 (MQ=255)
                                                                                                                             ccgccgGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGaa                                                                                                         <  1:221513/64‑1 (MQ=255)
                                                                                                                             ccgccgGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGaa                                                                                                         >  2:221513/1‑64 (MQ=255)
                                                                                                                             ccgccgGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGAt                    >  4:872970/1‑149 (MQ=255)
                                                                                                                                  gtgcaagcGCACCATAATGGTGCGAATTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTc                >  2:209817/8‑149 (MQ=255)
                                                                                                                                  gtgcaagcGCACCATAATGGTGCGAATTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTc                >  3:642361/8‑149 (MQ=255)
                                                                                                                                  ggCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCa               >  2:861647/1‑149 (MQ=255)
                                                                                                                                       ttCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATg          >  3:731767/1‑149 (MQ=255)
                                                                                                                                       ttCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATg          >  3:500022/1‑149 (MQ=255)
                                                                                                                                        tCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCATCATGc         <  2:489425/149‑1 (MQ=255)
                                                                                                                                        gcGCACCATAATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCGGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGc         >  4:408439/2‑149 (MQ=255)
                                                                                                                                           cACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATAGTCAACATGCGGt      >  2:668983/1‑149 (MQ=255)
                                                                                                                                              cATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGga   <  1:144064/149‑1 (MQ=255)
                                                                                                                                               aTTATGGTGCG‑CTTGCGCACCTTCAATGCCGGAGAGAAGACGTGAATGACGCTATCCCTCGCCCGCAGGCGTAAACCCGTAAGGCGGTGACGCCGCTATTAGAAATCCGCAACCTGACCAACTTCTACGATGGTCAACATGCGGTGgat  >  2:880334/1‑149 (MQ=255)
                                                                                                                                               aTTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAGCCTGACCAAAGCCTACGATGGTCAACATGCGGTGgat  <  4:830199/149‑1 (MQ=255)
                                                                                                                                               aTTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGgat  >  2:241559/1‑149 (MQ=255)
                                                                                                                                               aTTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGgat  >  1:33591/1‑149 (MQ=255)
                                                                                                                                               aTTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGgat  <  2:676973/149‑1 (MQ=255)
                                                                                                                                               aTTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGgat  <  3:484588/149‑1 (MQ=255)
                                                                                                                                               aTTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGgat  >  3:433019/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCG‑CTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGAT  >  NC_000913/894807‑895098

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