Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 585,866 G→A 10.6% intergenic (‑233/‑191) ompT ← / → pauD DLP12 prophage; outer membrane protease VII; outer membrane protein 3b/tRNA‑OTHER

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913585,8660GA10.6% 366.9 / 27.9 123intergenic (‑233/‑191)ompT/pauDDLP12 prophage; outer membrane protease VII; outer membrane protein 3b/tRNA‑OTHER
Reads supporting (aligned to +/- strand):  ref base G (57/53);  new base A (7/6);  total (64/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

Only 100 of 124 total aligned reads displayed.

TTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCA  >  NC_000913/585719‑586012
                                                                                                                                                   |                                                                                                                                                  
tttATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTaga                                                                                                                                                   <  4:321814/149‑1 (MQ=255)
           tttAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGc                                                                                                                                        >  2:914586/1‑149 (MQ=255)
           tttAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTCTAAAGGCatg                                                                                                                                             >  2:579895/1‑143 (MQ=255)
                  ttGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATa                                                                                                                                 >  2:262482/1‑149 (MQ=255)
                  ttGAATGATTATCACAATGCGCCATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATa                                                                                                                                 >  1:521005/1‑149 (MQ=255)
                         aTTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAg                                                                                                                          <  1:914586/149‑1 (MQ=255)
                          ttATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGc                                                                                                                               <  1:570652/143‑1 (MQ=255)
                          ttATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGc                                                                                                                               >  2:570652/1‑143 (MQ=255)
                            aTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATa                                                                                                                                 <  1:1061199/139‑1 (MQ=255)
                            aTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATa                                                                                                                                 >  2:1061199/1‑139 (MQ=255)
                             tCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGc                                                                                                                               <  2:899564/140‑1 (MQ=255)
                             tCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGc                                                                                                                               >  1:899564/1‑140 (MQ=255)
                                      gcTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATaa                                                                                                                                          >  2:831208/1‑120 (MQ=255)
                                      gcTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATaa                                                                                                                                          <  1:831208/120‑1 (MQ=255)
                                      gcTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAAt                                                                                                             >  2:1278133/1‑149 (MQ=255)
                                          tCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAg                                                                                                                          >  3:616387/1‑132 (MQ=255)
                                          tCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAg                                                                                                                          <  4:616387/132‑1 (MQ=255)
                                                    gTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTa                                                                                               <  1:1278133/149‑1 (MQ=255)
                                                        aaTTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtt                                                                                                    <  2:1240366/140‑1 (MQ=255)
                                                        aaTTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtt                                                                                                    <  2:458190/140‑1 (MQ=255)
                                                        aaTTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtt                                                                                                    >  1:458190/1‑140 (MQ=255)
                                                        aaTTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtt                                                                                                    >  1:1240366/1‑140 (MQ=255)
                                                          ttATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGc                                                                                                                                        >  3:1146753/1‑102 (MQ=255)
                                                          ttATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACaa                                                                                         >  1:484910/1‑149 (MQ=255)
                                                             tCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCaa                                                                                                              >  2:225004/1‑125 (MQ=255)
                                                             tCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtt                                                                                                    >  2:655695/1‑135 (MQ=255)
                                                             tCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtt                                                                                                    <  1:655695/135‑1 (MQ=255)
                                                              cTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACaaaaaa                                                                                     <  2:1070758/149‑1 (MQ=255)
                                                                gTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAaga                                                                                   <  1:262482/149‑1 (MQ=255)
                                                                 ttATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTCTAAAGGCatgg                                                                                                                                            <  1:1120126/91‑3 (MQ=255)
                                                                 ttATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTCTAAAGGCatgg                                                                                                                                            >  2:1120126/1‑89 (MQ=255)
                                                                 ttATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTCTAAAGGCat                                                                                                                                              >  1:1770140/1‑89 (MQ=255)
                                                                 ttATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTCTAAAGGCat                                                                                                                                              <  2:1770140/89‑1 (MQ=255)
                                                                    ttGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGa                                                                               <  2:484910/149‑1 (MQ=255)
                                                                       tttCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCa                                                                                                                              <  1:1656987/99‑1 (MQ=255)
                                                                                tCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGac                                                                   >  3:512120/1‑149 (MQ=255)
                                                                                    ttttACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCaa                                                                                                              >  3:361553/1‑102 (MQ=255)
                                                                                    ttttACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCaa                                                                                                              <  4:361553/102‑1 (MQ=255)
                                                                                    ttttACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAAt                                                                                                             >  1:1818966/1‑103 (MQ=255)
                                                                                    ttttACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAAt                                                                                                             <  2:1818966/103‑1 (MQ=255)
                                                                                    ttttACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAagag                                                                                  <  3:895/130‑1 (MQ=255)
                                                                                    ttttACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAagag                                                                                  >  4:895/1‑130 (MQ=255)
                                                                                         cTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGc                                                          <  4:512120/149‑1 (MQ=255)
                                                                                            tgtgGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGc                                                                                                                                        <  3:816165/68‑1 (MQ=255)
                                                                                            tgtgGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGc                                                                                                                               >  2:1283822/1‑77 (MQ=255)
                                                                                            tgtgGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGc                                                                                                                               <  1:1283822/77‑1 (MQ=255)
                                                                                            tgtgGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTc                                                       >  1:329479/1‑149 (MQ=255)
                                                                                                 tttttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTg                                                                       <  3:539633/128‑1 (MQ=255)
                                                                                                 tttttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAAt                                                  >  3:1343331/1‑149 (MQ=255)
                                                                                                    ttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGc                                               <  2:521005/149‑1 (MQ=255)
                                                                                                            ttttGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACaaaaaa                                                                                     >  1:1544720/1‑103 (MQ=255)
                                                                                                            ttttGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACaaaaaa                                                                                     <  2:1544720/103‑1 (MQ=255)
                                                                                                            ttttGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAgg                                       >  3:1012803/1‑149 (MQ=255)
                                                                                                            ttttGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAgg                                       >  2:957242/1‑149 (MQ=255)
                                                                                                                   gtgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTga                                                                      >  2:1230518/1‑111 (MQ=255)
                                                                                                                    tgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACaaaaaa                                                                                     <  2:219473/95‑1 (MQ=255)
                                                                                                                    tgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACaaaaaa                                                                                     >  1:219473/1‑95 (MQ=255)
                                                                                                                     gCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAggg                                                                                >  4:1022843/1‑99 (MQ=255)
                                                                                                                     gCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAggg                                                                                <  3:1022843/99‑1 (MQ=255)
                                                                                                                       ttttATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGGGAAAATAGCATGAGGt                                                                                                                        >  4:1104282/1‑57 (MQ=255)
                                                                                                                       ttttATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGt                                                                                                                        <  3:1104282/57‑1 (MQ=255)
                                                                                                                       ttttATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCaa                                                                                                              >  3:835496/1‑67 (MQ=255)
                                                                                                                       ttttATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCaa                                                                                                              <  4:835496/67‑1 (MQ=255)
                                                                                                                       ttttATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGttttt                                                                                                 <  1:203284/80‑1 (MQ=255)
                                                                                                                        tttATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtttttt                                                                                                <  4:368373/80‑1 (MQ=255)
                                                                                                                           aTCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGa                        <  4:1343331/149‑1 (MQ=255)
                                                                                                                           aTCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGa                        <  1:585504/149‑1 (MQ=255)
                                                                                                                           aTCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGa                        >  1:1676573/1‑149 (MQ=255)
                                                                                                                              tATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGc                                      >  2:1373303/1‑132 (MQ=255)
                                                                                                                               aTTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAagag                                                                                  >  1:1035844/1‑87 (MQ=255)
                                                                                                                               aTTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAagag                                                                                  <  1:1119613/87‑1 (MQ=255)
                                                                                                                               aTTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAagag                                                                                  <  2:1035844/87‑1 (MQ=255)
                                                                                                                               aTTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTa                    >  3:1716941/1‑149 (MQ=255)
                                                                                                                                tttAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATc                                                                                                                   >  3:565058/1‑53 (MQ=255)
                                                                                                                                tttAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATc                                                                                                                   <  4:565058/53‑1 (MQ=255)
                                                                                                                                tttAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAg                   >  4:1380998/1‑149 (MQ=255)
                                                                                                                                   aaGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTgg                <  2:1676573/149‑1 (MQ=255)
                                                                                                                                      ttataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCaa                                                   >  3:1883007/5‑111 (MQ=255)
                                                                                                                                      ttataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCaa                                                   <  4:1883007/107‑1 (MQ=255)
                                                                                                                                       tataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAaga                                                                                   >  1:350616/4‑78 (MQ=255)
                                                                                                                                       tataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAaga                                                                                   <  2:350616/75‑1 (MQ=255)
                                                                                                                                        ataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACgt           >  1:1476675/3‑149 (MQ=255)
                                                                                                                                          aTGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAagag                                                                                  <  3:1259491/76‑1 (MQ=255)
                                                                                                                                          aTGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAagag                                                                                  >  4:1259491/1‑76 (MQ=255)
                                                                                                                                          aTGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAaga                                                                                   >  4:647166/1‑75 (MQ=255)
                                                                                                                                          aTGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCt                                                        >  3:735332/1‑102 (MQ=255)
                                                                                                                                          aTGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCt                                                        <  4:735332/102‑1 (MQ=255)
                                                                                                                                          aTGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCt                                                        <  2:1696514/102‑1 (MQ=255)
                                                                                                                                          aTGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGa                        <  1:1251692/134‑1 (MQ=255)
                                                                                                                                          aTGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACgtgt         <  1:41609/149‑1 (MQ=255)
                                                                                                                                          aTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTACAAGGCTTTACATGTTTTGAGGTAGTGGGACgtgt         >  1:741248/1‑149 (MQ=255)
                                                                                                                                           tGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCACTTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTga                                                                      <  4:1101379/87‑1 (MQ=255)
                                                                                                                                           tGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGtgtg        >  4:859558/1‑149 (MQ=255)
                                                                                                                                            gCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCaa                                         <  3:991211/115‑1 (MQ=255)
                                                                                                                                            gCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCaa                                         >  4:991211/1‑115 (MQ=255)
                                                                                                                                            gCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTACAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGa       <  2:741248/149‑1 (MQ=255)
                                                                                                                                              cTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGttt                           <  1:1203593/127‑1 (MQ=255)
                                                                                                                                              cTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGttt                           >  2:1203593/1‑127 (MQ=255)
                                                                                                                                               tttAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAgcg    >  3:22767/1‑149 (MQ=255)
                                                                                                                                                 tAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCa  <  4:1012803/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCA  >  NC_000913/585719‑586012

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