Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 585,863 T→C 12.6% intergenic (‑230/‑194) 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,8630TC12.6% 230.5 / 20.8 80intergenic (‑230/‑194)ompT/pauDDLP12 prophage; outer membrane protease VII; outer membrane protein 3b/tRNA‑OTHER
Reads supporting (aligned to +/- strand):  ref base T (37/32);  new base C (5/5);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

TAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCCGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCC  >  NC_000913/585732‑586109
                                                                                                                                   |                                                                                                                                                                                                                                                      
tAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTagag                                                                                                                                                                                                                                                   >  5:871630/1‑137 (MQ=255)
tAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTagag                                                                                                                                                                                                                                                   <  6:871630/137‑1 (MQ=255)
   tATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACTTAAGGCCTCGCAATGGCTTCCAAGGCTTTAATTGTTTTGAGGTATTGGGACGTGTGAgc                                                                                                      >  4:469657/1‑275 (MQ=255)
        aaTGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATa                                                                                                                                                                                                                                  >  4:66202/1‑146 (MQ=255)
            aTTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCAt                                                                                                                                                                                                                              <  2:66203/146‑1 (MQ=255)
             ttATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCt                                                                                                                                               >  4:120000/1‑224 (MQ=255)
                tCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAAc                                                                                                                                                                          >  4:211219/1‑194 (MQ=255)
                 cacaATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCa                                                                                                                                                                                                                                                <  8:850182/123‑1 (MQ=255)
                 cacaATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCa                                                                                                                                                                                                                                                >  7:850182/1‑123 (MQ=255)
                 cacaATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCa                                                                                                                                           <  2:120001/224‑1 (MQ=255)
                    aaTGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTgag                                                                                                                                                                      <  2:211220/194‑1 (MQ=255)
                            aTCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTCGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCACGGCTTTACATGTTTTGAGGTAGTGGGATGTGTGAGCGCAGAGTTGGCGCGGTAAGttgttg                                                                             >  1:12936/1‑275 (MQ=255)
                                         ttAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGc                                                                                                                                                                                                                                                 <  8:392444/98‑1 (MQ=255)
                                         ttAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGc                                                                                                                                                                                                                                                 >  7:392444/1‑98 (MQ=255)
                                                  tGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAAc                                                                                                                                                                          >  4:571197/1‑160 (MQ=255)
                                                     tATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATg                                                                                                                                                                                                                             >  7:843655/1‑106 (MQ=255)
                                                     tATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATg                                                                                                                                                                                                                             <  8:843655/106‑1 (MQ=255)
                                                      aTTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTgag                                                                                                                                                                      <  2:571198/160‑1 (MQ=255)
                                                      aTTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACtt                                                    >  2:276524/1‑274 (MQ=255)
                                                          tttCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAAt                                                                                                                                                                                                              >  2:540644/1‑116 (MQ=255)
                                                           ttCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGa                                                                                                                                                                                                                            >  6:664961/1‑101 (MQ=255)
                                                           ttCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGa                                                                                                                                                                                                                            <  5:664961/101‑1 (MQ=255)
                                                            tCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGCGTTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAAGGTCGTTCTAGGTACTTCTAAGTc                                             >  1:52102/1‑275 (MQ=255)
                                                              atatTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACt                                                                                                                                                                                                          <  4:540643/116‑1 (MQ=255)
                                                                  ttCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtttttt                                                                                                                                                                                                 >  8:243098/1‑121 (MQ=255)
                                                                  ttCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtttttt                                                                                                                                                                                                 <  7:243098/121‑1 (MQ=255)
                                                                   tCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCa                                                                                                                                                                                                                <  5:152141/105‑1 (MQ=255)
                                                                   tCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCa                                                                                                                                                                                                                >  6:152141/1‑105 (MQ=255)
                                                                      gTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAaga                                                                                                                                                                                    >  8:202119/1‑130 (MQ=255)
                                                                      gTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAaga                                                                                                                                                                                    <  7:202119/130‑1 (MQ=255)
                                                                        tttACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTa                                                                                                                                                                                                            <  6:854462/104‑1 (MQ=255)
                                                                        tttACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTa                                                                                                                                                                                                            >  5:854462/1‑104 (MQ=255)
                                                                           aCTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTAt                                                                                                                                                                                                                     >  5:991206/1‑92 (MQ=255)
                                                                           aCTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTAt                                                                                                                                                                                                                     <  6:991206/92‑1 (MQ=255)
                                                                                 tgGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTa                                                                                                                                                                                                                      >  4:567643/1‑85 (MQ=255)
                                                                                   gTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGaaa                                                                                                                                                                                                                                     >  5:988371/1‑68 (MQ=255)
                                                                                   gTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGaaa                                                                                                                                                                                                                                     <  6:988371/68‑1 (MQ=255)
                                                                                    tttttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCata                                                                                                                                                                                                                                              <  8:7908/58‑1 (MQ=255)
                                                                                    tttttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCata                                                                                                                                                                                                                                              >  7:7908/1‑58 (MQ=255)
                                                                                     ttttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCt                                                                                                                                                                                                                  <  2:567644/85‑1 (MQ=255)
                                                                                           aTGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATaa                                                                                                                                                                                                                                           >  3:83366/1‑54 (MQ=255)
                                                                                               ttttGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGa                                                                                                                                                                                                                                       <  1:83366/54‑1 (MQ=255)
                                                                                                       tgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATg                                                                                                                                                                                                                             >  5:39254/1‑56 (MQ=255)
                                                                                                       tgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATg                                                                                                                                                                                                                             <  6:39254/56‑1 (MQ=255)
                                                                                                       tgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtttttt                                                                                                                                                                                                 <  8:633226/84‑1 (MQ=255)
                                                                                                       tgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtttttt                                                                                                                                                                                                 >  7:633226/1‑84 (MQ=255)
                                                                                                       tgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAaga                                                                                                                                                                                    >  6:390854/1‑97 (MQ=255)
                                                                                                       tgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAaga                                                                                                                                                                                    <  5:390854/97‑1 (MQ=255)
                                                                                                       tgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCCGCAGGTTCAAATCCTGCAGGGCGCGCCATTTCTTcc  <  4:276523/275‑1 (MQ=255)
                                                                                                               tCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTg                                                                                                                                                                        <  5:234692/101‑1 (MQ=255)
                                                                                                               tCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTg                                                                                                                                                                        >  6:234692/1‑101 (MQ=255)
                                                                                                                  aTTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGa                                                                                                                                                                                >  8:401296/1‑90 (MQ=255)
                                                                                                                  aTTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGa                                                                                                                                                                                <  7:401296/90‑1 (MQ=255)
                                                                                                                   tttAAGTGCAATGCCTGTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGAAAAAACTGAGACAAATAAGAACTCGCAATGGCTTGCAAGGCTTTACATGNTTTGAGGTAGTGGGACGTTTAAGAGCACAGATGGCGCGGTAAGTTGTTGACTTCAAATGTCGTTCTAGGAACTTCTAAGTCTGGGGCCGCACGTTCGacctct                       >  3:107184/1‑237 (MQ=255)
                                                                                                                    ttAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAagag                                                                                                                                                                                   >  7:1086114/1‑85 (MQ=255)
                                                                                                                    ttAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAagag                                                                                                                                                                                   <  8:1086114/85‑1 (MQ=255)
                                                                                                                       aGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCaa                                                                                                                                                    >  7:93505/1‑113 (MQ=255)
                                                                                                                       aGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCaa                                                                                                                                                    <  8:93505/113‑1 (MQ=255)
                                                                                                                          tataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTa                                                                                                                                                                                                <  5:263492/63‑1 (MQ=39)
                                                                                                                          tataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTa                                                                                                                                                                                                >  6:263492/4‑66 (MQ=39)
                                                                                                                          tataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACaaaa                                                                                                                                                                           <  6:445266/84‑1 (MQ=255)
                                                                                                                          tataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACaaaa                                                                                                                                                                           >  5:445266/4‑87 (MQ=255)
                                                                                                                           ataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGttt                                                                                                                            >  6:106537/3‑133 (MQ=255)
                                                                                                                           ataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGttt                                                                                                                            <  5:106537/131‑1 (MQ=255)
                                                                                                                            taTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTa                                                                                                                                                                                                            >  5:775282/2‑52 (MQ=38)
                                                                                                                            taTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTa                                                                                                                                                                                                            <  6:775282/51‑1 (MQ=38)
                                                                                                                            taTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtttttt                                                                                                                                                                                                 <  5:684060/62‑1 (MQ=255)
                                                                                                                            taTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtttttt                                                                                                                                                                                                 >  6:684060/2‑63 (MQ=255)
                                                                                                                            cATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAACGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACgt                                                                                                            >  5:522464/1‑148 (MQ=255)
                                                                                                                             aTGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAg                                                                                                                    >  8:815712/1‑139 (MQ=255)
                                                                                                                             aTGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAg                                                                                                                    <  7:815712/139‑1 (MQ=255)
                                                                                                                              tGCCTTTAGAGGCATATTAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAag                                                                                                                                                                                     <  5:283440/73‑1 (MQ=255)
                                                                                                                              tGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCTTCAATTACTATGTTTTTTAGTACAAAAAag                                                                                                                                                                                     >  6:283440/1‑73 (MQ=255)
                                                                                                                              tGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATg                                                                                                                               <  5:1036958/127‑1 (MQ=255)
                                                                                                                              tGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATg                                                                                                                               >  6:1036958/1‑127 (MQ=255)
                                                                                                                              tGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGa                                                                                                                         <  7:970047/133‑1 (MQ=255)
                                                                                                                              tGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGa                                                                                                                         >  8:970047/1‑133 (MQ=255)
                                                                                                                              tGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACgtgt                                                                                                          >  5:740600/1‑148 (MQ=255)
                                                                                                                                ccTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtttttt                                                                                                                                                                                                 <  8:546864/59‑1 (MQ=255)
                                                                                                                                ccTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtttttt                                                                                                                                                                                                 >  7:546864/1‑59 (MQ=255)
                                                                                                                                  tttAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTg                                                                                                                                                                        <  8:984288/82‑1 (MQ=255)
                                                                                                                                  tttAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTg                                                                                                                                                                        >  7:984288/1‑82 (MQ=255)
                                                                                                                                   ttAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCaa                                                                                                                                                    <  8:934520/101‑1 (MQ=255)
                                                                                                                                   ttAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCaa                                                                                                                                                    >  7:934520/1‑101 (MQ=255)
                                                                                                                                   ttAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGtgtg                                                                                                         >  6:248135/1‑144 (MQ=255)
                                                                                                                                   ttAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGtgtg                                                                                                         <  5:248135/144‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                                                                                                                                      
TAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCCGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCC  >  NC_000913/585732‑586109

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