Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,358,800 G→A 5.7% P130S (CCG→TCG)  puuP ← putrescine importer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,358,8000GA 5.7% 222.5 / 2.5 73P130S (CCG→TCG) puuPputrescine importer
Reads supporting (aligned to +/- strand):  ref base G (28/39);  new base A (2/2);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

ACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGA  >  NC_000913/1358652‑1358948
                                                                                                                                                    |                                                                                                                                                    
aCTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCATGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCggg                                                                                                                                                     <  2:1217886/150‑1 (MQ=255)
   aCCAGTAATAAAAACACCACCATGATGGAGATTTGCCCCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAag                                                                                                                                                  <  2:216879/150‑1 (MQ=255)
    ccAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAaga                                                                                                                                                 <  2:209178/150‑1 (MQ=255)
     cAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAagag                                                                                                                                                <  1:384061/150‑1 (MQ=255)
     cAGGAAGAGAAACACCACCATGAAGGAGATTTGCACCAGGACAAACAGGGCATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAagag                                                                                                                                                <  1:1372160/150‑1 (MQ=255)
             tAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGgaga                                                                                                                                        <  1:1076893/150‑1 (MQ=255)
              aaaCACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGgagag                                                                                                                                       <  2:531745/150‑1 (MQ=255)
               aaCACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGCGAAGAGGGCGgagaga                                                                                                                                      >  1:573079/1‑150 (MQ=255)
                   ccaccaTGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTAAGATTCGCGGCGGTTAAAATGGCTACGAAGGTTACCCCCAACACCAACGGCGGAACTTCCGAGAAGAGGGCGGAGagataga                                                                                                                                  >  2:470851/1‑150 (MQ=255)
                   ccaccaTGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTATTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCCCCCACCCCCACGGCGGCACTTCCGAGAAGCGGGCGgagcgcgaga                                                                                                                                  >  2:463870/1‑143 (MQ=255)
                        atgatgGAGATTTGAACCAGTACAAACAGGGTATTGAAGTTAGCGACCCGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTc                                                                                                                             <  1:924472/150‑1 (MQ=255)
                         tgatgGAGATTTTCACCAGTACAAAGAGGGTATTGAAGGTAGCGACCAGGTGGACGCTCTTCAGATTCGCGGCGGTTAAATTGACGACGATGGTTACAACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCg                                                                                                                            <  1:694517/150‑1 (MQ=255)
                         tgatgGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCg                                                                                                                            <  1:1144534/150‑1 (MQ=255)
                           atgGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGcc                                                                                                                          <  2:1202792/150‑1 (MQ=255)
                            tgGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTTCCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCca                                                                                                                         >  1:908044/1‑150 (MQ=255)
                            tgGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCca                                                                                                                         >  2:473431/1‑150 (MQ=255)
                             ggAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCcaa                                                                                                                        >  1:416162/1‑150 (MQ=255)
                              gagaTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAACTGGCGACGAAGGGTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCcaac                                                                                                                       <  2:1647374/150‑1 (MQ=255)
                              gagaTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACTACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCcaac                                                                                                                       >  2:1612970/1‑150 (MQ=255)
                              gagaTGTGCACCAGCACAAACAGCGTATTGAAGTAAGCTACCCGGTTGACGCTCTGCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCCCGGCGGCACTTCCGGGAAGAGGGAGGAGAGATAGAGTTTCGCcaac                                                                                                                       <  2:101419/150‑1 (MQ=255)
                               agaTTTGCACCAGTACAAACACGGTATTGAACTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCcaaca                                                                                                                      <  2:1252027/150‑1 (MQ=255)
                                   tttccccAGTCCAAACATCGTATTGAAGTTAGCGACCAGGTTGACTCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAATGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGAc                                                                                                                  <  2:998999/145‑1 (MQ=255)
                                      cACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGtt                                                                                                               >  1:70114/1‑150 (MQ=255)
                                          aGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAatc                                                                                                           <  2:869729/150‑1 (MQ=255)
                                              cAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCg                                                                                                       >  2:637179/1‑150 (MQ=255)
                                                  cAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCaa                                                                                                   >  1:1497990/1‑150 (MQ=255)
                                                   aGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCaaa                                                                                                  >  1:555663/1‑150 (MQ=255)
                                                    gggTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCTGCaaaa                                                                                                 >  2:362824/1‑150 (MQ=255)
                                                      gTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAag                                                                                               <  1:1315220/150‑1 (MQ=255)
                                                      gTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAag                                                                                               <  2:1663607/150‑1 (MQ=255)
                                                        agtGACGTTAGCGACCAGGTGGCCGCTCTTCAGAGTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAagag                                                                                             <  2:614844/148‑1 (MQ=255)
                                                        aTTGAATTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAagag                                                                                             <  2:1279434/150‑1 (MQ=255)
                                                         ttGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAgaga                                                                                            >  2:125525/1‑150 (MQ=255)
                                                           gAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCAACACCAACGGCTGCACTTCCTGGAAGAGGGAGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATa                                                                                          <  1:37837/150‑1 (MQ=255)
                                                               ttAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATcc                                                                                      <  2:240227/150‑1 (MQ=255)
                                                                  gCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCcagc                                                                                   <  1:1288378/150‑1 (MQ=255)
                                                                    gACCATGTAGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGCTACGAAGGTTACCACCCACACCCACGGCGGCTCTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACTTTAATCATAGGCAAAAAGAGATAATCcagaag                                                                                 >  2:1445490/1‑150 (MQ=255)
                                                                     aCCTGGTTGACGCTATTCAGATTCGCGGCGGGTAAAATGGCGACGATTGTTACCACCCACACACACGGCGGCACTTCCGGGAATATGGCGGAGATATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCagcagc                                                                                <  1:153288/150‑1 (MQ=255)
                                                                     aCCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCagcagc                                                                                <  1:417092/150‑1 (MQ=255)
                                                                         ggTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATg                                                                            >  1:1603152/1‑150 (MQ=255)
                                                                           ttGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGAc                                                                          <  1:1418057/150‑1 (MQ=255)
                                                                              aCGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAg                                                                       >  2:375612/1‑150 (MQ=255)
                                                                              aCGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAg                                                                       <  1:473431/150‑1 (MQ=255)
                                                                                gCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGcc                                                                     <  2:169189/150‑1 (MQ=255)
                                                                                   ctTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGAc                                                                  >  2:1833685/1‑150 (MQ=255)
                                                                                    ttCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTACGGGAAGAGGGAGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGCTAATCCAGCAGCGATGACCAGCCGAcc                                                                 >  1:1238973/1‑150 (MQ=255)
                                                                                            cgcgGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATcc                                                         <  1:1740822/150‑1 (MQ=255)
                                                                                            cgcgGCGGGTAAAATGGCGACGAACGTGACCCCCCACACCCACCGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACCAGACGTTAATCATCGGCAAAACGCGCTAATCCAGCAGCCATGACCAGCCGCCCATAAATcc                                                         <  2:1191773/150‑1 (MQ=255)
                                                                                             gcggcgGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCg                                                        <  1:1775251/150‑1 (MQ=255)
                                                                                              cggcggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGa                                                       <  1:637179/150‑1 (MQ=255)
                                                                                              cggcggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGAGAATCCAGCAGCGATGACCAGCCGACCATAAATCCGa                                                       >  1:1133198/1‑150 (MQ=255)
                                                                                                  ggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTg                                                   <  1:1288103/150‑1 (MQ=255)
                                                                                                      aaaaTGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCggg                                               <  2:640934/150‑1 (MQ=255)
                                                                                                      aaaaTGGCGACGAAGGTTACCACCAACACCCACGGCGGCACTTCCGAGAAGAGGGCGGAGAGATAGATTTTCGACAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCggg                                               <  1:470851/150‑1 (MQ=255)
                                                                                                      aaaaTGGCGACGAAGGGTACCACCCACACCCCCGGCGGCACTTCCGAGCAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCggg                                               <  1:463870/150‑1 (MQ=255)
                                                                                                        aaTGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGtt                                             >  2:1561870/1‑150 (MQ=255)
                                                                                                             cgacgaAGGTTACCACCCACACCCACGGCGGCACTTCCGGGTAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCg                                        >  1:1178419/1‑150 (MQ=255)
                                                                                                                   aGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAATACGTAAACCATCGGCAATAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCCACGTGCGGGTTAATCGACtttt                                  <  2:914705/150‑1 (MQ=255)
                                                                                                                    ggTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTg                                 >  2:988971/1‑150 (MQ=255)
                                                                                                                    ggTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTT‑ATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTgc                                >  2:129587/1‑150 (MQ=255)
                                                                                                                      gtACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCGGACCATAAATCCGACGTGCGGGTTAATCGACTTTTgcg                               <  2:1791576/149‑1 (MQ=255)
                                                                                                                       tACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTgcgc                              <  2:404367/150‑1 (MQ=255)
                                                                                                                          cacccacacccaCGGCGGCACTTCCGGGAAGAGGGCGGAGAGAGAGATTGGCGCCAACACGCCGGTAAGCCTCGGCAAAAAGAGATAAGCCCGCCGCGATGACCAGCCGACCATAAAGCCGCCTTGCGGGTGAATCGACGTTTGCGCGTa                           >  1:1297548/1‑150 (MQ=255)
                                                                                                                             ccacacccacGTCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGt                        <  1:1446235/150‑1 (MQ=255)
                                                                                                                              cccccccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCAACGGCAACAAGAGATAATCCAGCAGCGATGCCCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGta                       >  2:644997/5‑150 (MQ=255)
                                                                                                                              cacacccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGta                       <  2:1094733/150‑1 (MQ=255)
                                                                                                                                  cccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAgg                   >  2:1066886/1‑150 (MQ=255)
                                                                                                                                   ccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGCTTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTGAATCGACTTTTGCGCGTAGGTATCGGc                  >  1:177907/1‑150 (MQ=255)
                                                                                                                                    cacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGcc                 <  1:125525/150‑1 (MQ=255)
                                                                                                                                        gcggcACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAAc             >  1:1738368/1‑150 (MQ=255)
                                                                                                                                              cTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAAGCCAGCAGCGATGCCCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCt       >  2:302352/1‑150 (MQ=255)
                                                                                                                                              cTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGACAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTGTTGCGCGTAGGTATAGGCCGAACCGGCCt       >  1:1287073/1‑150 (MQ=255)
                                                                                                                                               ttCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTc      >  1:1549327/1‑150 (MQ=255)
                                                                                                                                                 ccGTGAAGATGGCGGAGATAAAGATTTTCCCCAACAAGACCTTAATCATCGGCAAAAAGAGATAATCCAGCATCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCg    <  2:908044/150‑1 (MQ=255)
                                                                                                                                                   gggAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCTGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGGTTCGGCCtatacc  <  2:70114/150‑6 (MQ=255)
                                                                                                                                                   gggAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGa  >  1:1026612/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGA  >  NC_000913/1358652‑1358948

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