Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 399,810 A→C 5.7% intergenic (+5/+19) yaiZ → / ← ddlA DUF2754 family putative inner membrane protein/D‑alanine‑D‑alanine ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913399,8100AC 5.7% 374.1 / 12.0 124intergenic (+5/+19)yaiZ/ddlADUF2754 family putative inner membrane protein/D‑alanine‑D‑alanine ligase A
Reads supporting (aligned to +/- strand):  ref base A (47/70);  new base C (4/3);  total (51/73)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01

Only 100 of 132 total aligned reads displayed.

CTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTG  >  NC_000913/399662‑399958
                                                                                                                                                    |                                                                                                                                                    
cTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAata                                                                                                                                                      >  1:680150/1‑149 (MQ=255)
 ttAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAatat                                                                                                                                                     >  2:1460777/1‑149 (MQ=255)
  tAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAAtatt                                                                                                                                                    >  2:784770/1‑149 (MQ=255)
           ggtggGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGt                                                                                                                                           <  1:1216746/149‑1 (MQ=255)
            gtggGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTAGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGtt                                                                                                                                          >  1:108400/1‑149 (MQ=255)
            gtggGGGTACTGGGATTTGAAGAAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGTGTTATTAGNTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGAGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGtt                                                                                                                                          <  1:1103855/149‑1 (MQ=255)
                ggTTACTTGGATTTGAAGCAAAAGGTTGGCAGACCTATGGCGTCGGTCTGGTGACGTGTGTGATTAGGTTCTGGCTGGCGGTGTTGATTATTCGTCGTCGCGATGAAGAAACTGTAAAAGCCCGATAATATTCTTATTGGGCGttttca                                                                                                                                      <  2:898261/149‑5 (MQ=255)
                ggTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTaaaa                                                                                                                                      <  1:486581/149‑1 (MQ=255)
                 ggtACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCGATAATATTCTTATTGGGCGttttcag                                                                                                                                     <  2:1327031/147‑6 (MQ=255)
                  ttACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATTCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAata                                                                                                                                    <  2:330036/149‑1 (MQ=255)
                     cTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGGTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTa                                                                                                                                 <  1:1409681/149‑1 (MQ=255)
                      tGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTAc                                                                                                                                >  2:771037/1‑149 (MQ=255)
                      tGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTAc                                                                                                                                >  2:776329/1‑149 (MQ=255)
                      ggggATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGGTGATTAATCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTAc                                                                                                                                <  1:44582/148‑1 (MQ=255)
                        ggATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATTAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACAt                                                                                                                              >  2:128430/1‑149 (MQ=255)
                           tttGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATtgt                                                                                                                           >  2:452811/1‑149 (MQ=255)
                            ttGAAGCAAAAGGTTGGCAGACCTATGCCGTCTGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATtgtg                                                                                                                          <  1:1190272/149‑1 (MQ=255)
                            ttGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATtgtg                                                                                                                          <  1:54935/149‑1 (MQ=255)
                            ttGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATtgtg                                                                                                                          >  1:402/1‑149 (MQ=255)
                             tGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTgg                                                                                                                         <  1:1460777/149‑1 (MQ=255)
                             tGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTgg                                                                                                                         <  2:62168/149‑1 (MQ=255)
                               aaGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGtt                                                                                                                       >  2:1136711/1‑149 (MQ=255)
                                 gCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGtttt                                                                                                                     >  2:511063/1‑149 (MQ=255)
                                 gCAAAAGGTTGGCAGACATATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCCCCATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGtttt                                                                                                                     >  2:1048202/1‑149 (MQ=255)
                                  cAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTc                                                                                                                    >  1:226317/1‑149 (MQ=255)
                                  cAAAAGGTTGGCAGACCTATGCCATCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTc                                                                                                                    >  1:459132/1‑149 (MQ=255)
                                  cAAAAGGGTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTc                                                                                                                    >  1:1555405/1‑149 (MQ=255)
                                    aaaGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCaa                                                                                                                  <  2:944515/149‑1 (MQ=255)
                                    aaaGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCaa                                                                                                                  >  1:923868/1‑149 (MQ=255)
                                    aaaGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCaa                                                                                                                  >  1:694723/1‑149 (MQ=255)
                                    aaaGGTTGGCAGACCTATGCCGTCGGTCTGGGGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCaa                                                                                                                  <  1:515793/149‑1 (MQ=255)
                                       ggTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGc                                                                                                               >  2:1228192/1‑149 (MQ=255)
                                       ggTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGc                                                                                                               >  1:1084541/1‑149 (MQ=255)
                                             cAGACCTATGCCTGCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTAATCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATc                                                                                                         <  1:809063/149‑1 (MQ=255)
                                                aCCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGc                                                                                                      >  2:1595141/1‑149 (MQ=255)
                                                     tGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGt                                                                                                 >  2:1561465/1‑149 (MQ=255)
                                                        cGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGGCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGAc                                                                                              <  2:1589882/149‑1 (MQ=255)
                                                          tCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGc                                                                                            <  2:921784/149‑1 (MQ=255)
                                                                tGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAgc                                                                                      >  1:377687/1‑149 (MQ=255)
                                                                 ggTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTGGATTATTCGGCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAgcg                                                                                     <  2:625181/149‑1 (MQ=255)
                                                                  gTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAgcgc                                                                                    <  1:784770/149‑1 (MQ=255)
                                                                  gTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAgcgc                                                                                    >  1:268783/1‑149 (MQ=255)
                                                                    gACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCa                                                                                  <  1:1136711/149‑1 (MQ=255)
                                                                       gTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGtt                                                                               <  1:378297/149‑1 (MQ=255)
                                                                          ggTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTtcaa                                                                            <  1:432650/149‑1 (MQ=255)
                                                                          ggTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTtcaa                                                                            <  2:680150/149‑1 (MQ=255)
                                                                           tggATTAGTTTCTGGCTGGCTGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAATACGCCCAATAAGAATATTATCGGGCGTTAAAATATAACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTtcaat                                                                           <  1:239852/147‑1 (MQ=255)
                                                                               ttAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCaaa                                                                       >  1:505658/1‑149 (MQ=255)
                                                                               ttAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCaaa                                                                       <  1:1157948/149‑1 (MQ=255)
                                                                                tAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTTGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAAc                                                                      <  1:776329/149‑1 (MQ=255)
                                                                                tAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAAc                                                                      <  2:1181211/149‑1 (MQ=255)
                                                                                tAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAAc                                                                      >  1:478559/1‑149 (MQ=255)
                                                                                tAGTTTCTGGCTGGCGGGGTTGATTATTCGGCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAAc                                                                      <  1:771037/149‑1 (MQ=255)
                                                                                 aGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACg                                                                     >  1:734475/1‑149 (MQ=255)
                                                                                   tttCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACgtg                                                                   <  2:1140763/149‑1 (MQ=255)
                                                                                     tCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGtgtg                                                                 >  1:381912/1‑149 (MQ=255)
                                                                                       tggctggcGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGAt                                                               <  2:734475/149‑1 (MQ=255)
                                                                                             gcGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATc                                                         >  2:1555977/1‑149 (MQ=255)
                                                                                             gcGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATc                                                         <  2:986602/149‑1 (MQ=255)
                                                                                                gtgTTTATTATTCGTCGTAGAGATGAAGAAACTGAAATCGCACAATAAGAATATTTTCGGGCGTTAAAATATTTCATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGgt                                                      <  1:198328/147‑1 (MQ=255)
                                                                                                gggTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGgt                                                      <  2:407118/149‑1 (MQ=255)
                                                                                                gggTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGgt                                                      >  2:1423295/1‑149 (MQ=255)
                                                                                                 ggTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATTCGTTATCCGCAGCGTGACGCTCCATCTCCAGTTCAATCAAACGTGTGATCAGATCGgtt                                                     >  2:1130519/1‑148 (MQ=255)
                                                                                                   ttGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTa                                                   >  2:713649/1‑149 (MQ=255)
                                                                                                   ttGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTa                                                   >  2:1207725/1‑149 (MQ=255)
                                                                                                    tGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTaa                                                  <  1:1087122/149‑1 (MQ=255)
                                                                                                      attattCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAAcc                                                <  2:459132/149‑1 (MQ=255)
                                                                                                      attattCGTCGTCGCGATGAAGAAACTGAAAACGCCAAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAAcc                                                <  2:108400/149‑1 (MQ=255)
                                                                                                         attCTTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAAt                                                                                                                 <  1:993866/81‑1 (MQ=255)
                                                                                                              tcgtcgCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGc                                        <  2:135442/149‑1 (MQ=255)
                                                                                                                   gcgATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCt                                   <  1:128430/149‑1 (MQ=255)
                                                                                                                    cgATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCtt                                  >  1:1588712/1‑149 (MQ=255)
                                                                                                                    cgATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCtt                                  <  2:1097949/149‑1 (MQ=255)
                                                                                                                    cgATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCtt                                  >  1:159924/1‑149 (MQ=255)
                                                                                                                     gATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTtcaa                                                                            <  2:804786/106‑1 (MQ=255)
                                                                                                                     gATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTtcaa                                                                            >  1:804786/1‑106 (MQ=255)
                                                                                                                     gATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTg                                 >  2:791819/1‑149 (MQ=255)
                                                                                                                     gATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTg                                 >  1:219121/1‑149 (MQ=255)
                                                                                                                            aaaCTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCGATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGc                          >  2:948530/1‑149 (MQ=255)
                                                                                                                              tattttaaCGCCCGATAATATTCTTATTGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGTTGTAACCCAGACCGCTGGCTTGCCACAGCtt                        >  2:709033/7‑149 (MQ=38)
                                                                                                                               cTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTc                       <  2:694723/149‑1 (MQ=255)
                                                                                                                               cTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTc                       <  2:923868/149‑1 (MQ=255)
                                                                                                                                tGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCg                      <  2:325061/149‑1 (MQ=255)
                                                                                                                                  aaaaCGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTgg                                     <  1:799131/132‑1 (MQ=255)
                                                                                                                                  aaaaCGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTgg                                     >  2:799131/1‑132 (MQ=255)
                                                                                                                                  aaaaCGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGa                    <  2:203082/149‑1 (MQ=255)
                                                                                                                                  aaaaCGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGa                    <  2:346797/149‑1 (MQ=255)
                                                                                                                                    aaCGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGGTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCGTCGGata                  <  1:64564/149‑1 (MQ=255)
                                                                                                                                     aCGCCCGATAATATTCTTATTGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGatac                 >  1:973294/1‑149 (MQ=255)
                                                                                                                                        cccGATAATATTCTTATTGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGatacata              >  1:135772/1‑149 (MQ=255)
                                                                                                                                         ccAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGatacatac             <  1:1228192/149‑1 (MQ=255)
                                                                                                                                          cAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACt            <  2:687298/149‑1 (MQ=255)
                                                                                                                                          aaaTAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACt            <  2:1417238/148‑1 (MQ=255)
                                                                                                                                           aaTAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTg           >  2:1560328/1‑149 (MQ=255)
                                                                                                                                            aTAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACAGACTGa          >  1:96063/1‑149 (MQ=255)
                                                                                                                                              aaGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGata        >  2:1249981/1‑149 (MQ=255)
                                                                                                                                              aaGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGata        <  2:1277411/149‑1 (MQ=255)
                                                                                                                                                gAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATAtt      >  2:1321225/1‑149 (MQ=255)
                                                                                                                                                 aaTATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTa     <  2:402/149‑1 (MQ=255)
                                                                                                                                                   tattatCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGt   <  1:1236006/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTG  >  NC_000913/399662‑399958

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: