Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,356,445 A→T 7.2% intergenic (‑34/+25) dtpC ← / ← cadA dipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,356,4450AT 7.2% 651.9 / 32.2 195intergenic (‑34/+25)dtpC/cadAdipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand):  ref base A (93/88);  new base T (8/6);  total (101/94)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

Only 100 of 201 total aligned reads displayed.

GCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAG  >  NC_000913/4356297‑4356592
                                                                                                                                                    |                                                                                                                                                   
 cTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCaa                                                                                                                                                    <  4:114079/149‑1 (MQ=255)
      gggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCAc                                                                                                                                               <  2:1163192/149‑1 (MQ=255)
      gggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCAc                                                                                                                                               <  1:1274089/149‑1 (MQ=255)
         tGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTc                                                                                                                                            <  3:512292/149‑1 (MQ=255)
         tGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTc                                                                                                                                            <  4:745629/149‑1 (MQ=255)
         tGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTc                                                                                                                                            <  2:912899/149‑1 (MQ=255)
             aTAGAGAATGAGTAAGGCTCGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                                                                          <  2:105443/147‑1 (MQ=255)
             aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCtt                                                                                                                                        >  2:1216074/1‑149 (MQ=255)
               agagAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGt                                                                                                                                      >  4:1239424/1‑149 (MQ=255)
               agagAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGt                                                                                                                                      >  4:619527/1‑149 (MQ=255)
                 agaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTAc                                                                                                                                    >  4:1287735/1‑149 (MQ=255)
                  gaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACg                                                                                                                                   >  1:534082/1‑149 (MQ=255)
                   aaTGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGa                                                                                                                                  >  3:666799/1‑149 (MQ=255)
                   aaTGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGa                                                                                                                                  >  1:388454/1‑149 (MQ=255)
                      gAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCt                                                                                                                               >  3:562373/1‑149 (MQ=255)
                            ggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAattat                                                                                                                         >  1:436277/1‑149 (MQ=255)
                                 gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttttt                                                                                                                    <  2:134344/149‑1 (MQ=255)
                                        gTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttct                                                                                                             >  3:1285704/1‑149 (MQ=255)
                                        gTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttct                                                                                                             >  4:88210/1‑149 (MQ=255)
                                         tAAAAACTGAATTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttctt                                                                                                            >  3:842012/1‑149 (MQ=255)
                                           aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTtcttct                                                                                                          >  3:362665/1‑149 (MQ=255)
                                           aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTtcttct                                                                                                          <  3:1239424/149‑1 (MQ=255)
                                             aaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                                                                        >  1:1170123/1‑149 (MQ=255)
                                              aCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                                                                       <  1:360405/149‑1 (MQ=255)
                                                tGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCaa                                                                                                     >  2:404425/1‑149 (MQ=255)
                                                   aGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAc                                                                                                  >  1:126090/1‑149 (MQ=255)
                                                   aGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAc                                                                                                  <  1:1216074/149‑1 (MQ=255)
                                                        tCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCAc                                                                                                                                               <  3:1124180/99‑1 (MQ=255)
                                                        tCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAattatt                                                                                                                        <  1:208898/122‑1 (MQ=255)
                                                        tCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTaa                                                                                             >  1:1123426/1‑149 (MQ=255)
                                                        tCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTaa                                                                                             >  3:939429/1‑149 (MQ=255)
                                                        tCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTaa                                                                                             >  4:535934/1‑149 (MQ=255)
                                                         cccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAAc                                                                                            <  2:1298428/149‑1 (MQ=255)
                                                         cccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAAc                                                                                            >  4:241223/1‑149 (MQ=255)
                                                         cccAGACTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAAc                                                                                            <  2:131687/149‑1 (MQ=255)
                                                          ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGt                                                                                                                                      <  1:1307682/106‑1 (MQ=255)
                                                          ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACg                                                                                           >  1:1000903/1‑149 (MQ=255)
                                                          ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACg                                                                                           >  2:3230/1‑149 (MQ=255)
                                                          ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACg                                                                                           >  3:75752/1‑149 (MQ=255)
                                                          ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACg                                                                                           >  4:1384882/1‑149 (MQ=255)
                                                               tttGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGtata                                                                                      <  3:377174/149‑1 (MQ=255)
                                                                  ggATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAat                                                                                                                            >  2:510672/1‑108 (MQ=255)
                                                                    aTCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTa                                                                                                                                     <  3:887368/97‑1 (MQ=255)
                                                                    aTCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGc                                                                                 >  3:339719/1‑149 (MQ=255)
                                                                    aTCGCCACGATATAGTATATCGCGCGCGGCTGAGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTaa                                                                                             <  4:529094/137‑1 (MQ=255)
                                                                     tCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGcc                                                                                <  2:388454/149‑1 (MQ=255)
                                                                      cGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCa                                                                               <  2:1418494/149‑1 (MQ=255)
                                                                       gCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttt                                                                                                                      >  4:1340683/1‑109 (MQ=255)
                                                                         cACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                                                                            <  4:362665/149‑1 (MQ=255)
                                                                         cACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                                                                            >  3:1329809/1‑149 (MQ=255)
                                                                                  gTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACgg                                                                   <  1:501939/149‑1 (MQ=255)
                                                                                     tatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGt                                                                  <  4:1153425/147‑1 (MQ=255)
                                                                                     tatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGt                                                                  >  3:1153425/1‑147 (MQ=255)
                                                                                     tatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTAt                                                                <  2:999533/149‑1 (MQ=255)
                                                                                      atCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGcc                                                                                >  4:106658/1‑132 (MQ=255)
                                                                                                gCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAat                                                     <  4:842012/149‑1 (MQ=255)
                                                                                                  tgtgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTccc                                                                                                                                          <  2:581880/62‑1 (MQ=255)
                                                                                                     gAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCgg                                                <  1:907173/149‑1 (MQ=255)
                                                                                                      aGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGt                                               <  2:876079/149‑1 (MQ=255)
                                                                                                        ggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTtcttct                                                                                                          >  2:742124/1‑88 (MQ=255)
                                                                                                        ggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTtcttct                                                                                                          <  1:742124/88‑1 (MQ=255)
                                                                                                        ggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc                                                            <  1:478656/134‑1 (MQ=255)
                                                                                                         gtgtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTc                                            >  2:1073997/1‑149 (MQ=255)
                                                                                                         gtgtTTTCATGTGTCATCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGTCCATCAGCCTGTCGGTATTCACCGTGAATATCGGTTTc                                            >  4:1309518/1‑149 (MQ=255)
                                                                                                           gtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGc                                                                                 >  2:1264417/1‑110 (MQ=255)
                                                                                                               tCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                                                                        >  1:594286/1‑83 (MQ=255)
                                                                                                               tCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                                                                        <  2:594286/83‑1 (MQ=255)
                                                                                                               tCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGcc                                      <  3:388389/149‑1 (MQ=255)
                                                                                                                cATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGccc                                     <  1:785903/149‑1 (MQ=255)
                                                                                                                 aTGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCg                                    >  4:899875/1‑149 (MQ=255)
                                                                                                                 aTGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCg                                    <  4:565162/149‑1 (MQ=255)
                                                                                                                 aTGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCg                                    <  2:616321/149‑1 (MQ=255)
                                                                                                                 aTGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCg                                    >  3:582907/1‑149 (MQ=255)
                                                                                                                   gtgtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGa                                  <  4:582907/149‑1 (MQ=255)
                                                                                                                     gtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCTCTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATTCCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATa                                <  2:1448027/149‑1 (MQ=255)
                                                                                                                     gtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATa                                <  2:534082/149‑1 (MQ=255)
                                                                                                                     gtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATa                                >  4:815222/1‑149 (MQ=255)
                                                                                                                      ttCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTAt                                                                >  4:828516/1‑116 (MQ=255)
                                                                                                                      ttCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTAt                                                                <  3:828516/116‑1 (MQ=255)
                                                                                                                      ttCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAg                               >  1:201598/1‑149 (MQ=255)
                                                                                                                       tctcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTACCGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGt                              <  2:343417/149‑1 (MQ=255)
                                                                                                                         tcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATATCTTAACGGTATAGCGGCCATCAGCCTGACGGTa                                                                 >  4:1393348/1‑112 (MQ=255)
                                                                                                                         tcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTa                                                                 <  3:1393348/112‑1 (MQ=255)
                                                                                                                         tcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGa                            <  4:25119/149‑1 (MQ=255)
                                                                                                                         tcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGa                            <  1:253507/149‑1 (MQ=255)
                                                                                                                                  gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGAtt                   <  3:811186/149‑1 (MQ=255)
                                                                                                                                    aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGa                            >  1:233675/1‑138 (MQ=255)
                                                                                                                                    aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGa                            <  2:233675/138‑1 (MQ=255)
                                                                                                                                    aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTc                 >  4:142292/1‑149 (MQ=255)
                                                                                                                                    aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTc                 >  1:1421828/1‑149 (MQ=255)
                                                                                                                                     tacaaGGGAAGTGGCTTGCCACTTCCCTTTTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTca                >  2:410342/4‑149 (MQ=255)
                                                                                                                                     aaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGTGCCGATTTca                >  4:349024/1‑149 (MQ=255)
                                                                                                                                     aaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTca                >  3:738216/1‑149 (MQ=255)
                                                                                                                                          gggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAg           >  3:235492/1‑149 (MQ=255)
                                                                                                                                           ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCa                                           >  2:203979/1‑116 (MQ=255)
                                                                                                                                           ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGccc                                     >  1:981262/1‑122 (MQ=255)
                                                                                                                                            gAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAgta         <  3:529619/149‑3 (MQ=255)
                                                                                                                                             aaGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCAt        >  2:446497/1‑149 (MQ=255)
                                                                                                                                                  gCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCa   >  2:1479600/1‑149 (MQ=255)
                                                                                                                                                    aaGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGcc                                                                                <  4:494757/70‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAG  >  NC_000913/4356297‑4356592

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