Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 268,865 T→A 5.3% L32F (TTA→TTT ykfA ← CP4‑6 prophage; putative GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913268,8650TA 5.3% 303.5 / 4.2 93L32F (TTA→TTTykfACP4‑6 prophage; putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (47/41);  new base A (2/3);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCC  >  NC_000913/268725‑269008
                                                                                                                                            |                                                                                                                                               
aGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGt                                                                                                                                          >  8:140204/1‑148 (MQ=255)
  gCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTaa                                                                                                                                       >  3:197854/1‑149 (MQ=255)
     aaaGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCACCACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAAc                                                                                                                                              >  4:332552/1‑139 (MQ=255)
     aaaGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAAc                                                                                                                                              <  3:332552/139‑1 (MQ=255)
          gcAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTtgctg                                                                                                                                >  6:125637/1‑148 (MQ=255)
          gcAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTtgctg                                                                                                                                >  3:323128/1‑148 (MQ=255)
            aGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTtgctg                                                                                                                                <  1:306250/146‑1 (MQ=255)
            aGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTtgctg                                                                                                                                >  2:306250/1‑146 (MQ=255)
              tGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAg                                                                                                                             >  3:97446/1‑147 (MQ=255)
                gCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGc                                                                                                                         <  3:210216/149‑1 (MQ=255)
                  tCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGccc                                                                                                                       >  2:403605/1‑148 (MQ=255)
                  tCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCg                                                                                                                       >  2:225293/1‑149 (MQ=255)
                    aCGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGct                                                                                                                     >  7:341401/1‑149 (MQ=255)
                     cGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAAGAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGctc                                                                                                                    <  4:323128/149‑1 (MQ=255)
                     cGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGctc                                                                                                                    <  7:125033/149‑1 (MQ=255)
                     cGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGctc                                                                                                                    <  4:97446/149‑1 (MQ=255)
                      gtgtACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAg                                                                                                                          <  4:127919/142‑1 (MQ=255)
                      gtgtACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAg                                                                                                                          >  3:127919/1‑142 (MQ=255)
                      gtgtACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGctct                                                                                                                   >  4:156323/1‑149 (MQ=255)
                       tgtACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTg                                                                                                                  <  7:140204/149‑1 (MQ=255)
                          aCAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATa                                                                                                               >  7:218701/1‑149 (MQ=255)
                            aGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATAcc                                                                                                             >  1:242495/1‑149 (MQ=255)
                            aGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATAcc                                                                                                             <  5:213804/149‑1 (MQ=255)
                                 gccACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTagttagtt                                                                                                                                                 <  7:301969/108‑1 (MQ=255)
                                 gccACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTagttagtt                                                                                                                                                 >  8:301969/1‑108 (MQ=255)
                                    acaTCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACCCCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGTTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCgg                                                                                                     <  7:159080/149‑1 (MQ=255)
                                    acaTCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCgg                                                                                                     <  2:76104/149‑1 (MQ=255)
                                           tGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGAt                                                                                                >  8:320495/1‑147 (MQ=255)
                                             aCCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCa                                                                                              <  3:156323/147‑1 (MQ=255)
                                             aCCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAg                                                                                             >  3:89766/1‑148 (MQ=255)
                                             aCCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGc                                                                                            <  8:88286/149‑1 (MQ=255)
                                                ggTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAg                                                                                         <  1:326573/149‑1 (MQ=255)
                                                 gTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGc                                                                                        >  3:118394/1‑149 (MQ=255)
                                                   gATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACACATTGCGGCAGATCAGCAAGCGa                                                                                      >  4:91811/1‑149 (MQ=255)
                                                    atatatCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGAt                                                                                     <  8:87955/149‑1 (MQ=255)
                                                     tatatCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGAtt                                                                                    >  5:186768/1‑149 (MQ=255)
                                                     tatatCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTAGTAGTTAGTTAACTTGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGAtt                                                                                    <  5:40788/149‑1 (MQ=255)
                                                       tatCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTg                                                                                   <  4:89766/148‑1 (MQ=255)
                                                       tatCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGc                                                                                  <  6:186180/149‑1 (MQ=255)
                                                           tcCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCgg                                                                                                     >  8:250333/1‑126 (MQ=255)
                                                           tcCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCgg                                                                                                     <  7:250333/126‑1 (MQ=255)
                                                             cGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTtg                                                                                                                                   >  4:192816/1‑94 (MQ=255)
                                                             cGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTtg                                                                                                                                   <  3:192816/94‑1 (MQ=255)
                                                                     aGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCa                                                                     >  5:370239/1‑148 (MQ=255)
                                                                           ttGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCtt                                                               >  7:339697/1‑148 (MQ=255)
                                                                           ttGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTc                                                              <  1:225293/149‑1 (MQ=255)
                                                                            tGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTtgct                                                                                                                                 >  6:165415/1‑81 (MQ=255)
                                                                             gCACAGACTGCTCTTCCCAACCCCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTc                                                              >  6:187341/1‑147 (MQ=255)
                                                                             gCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTc                                                              <  5:187341/147‑1 (MQ=255)
                                                                              cacaGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAg                                                            <  3:91811/148‑1 (MQ=255)
                                                                              cacaGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCAATGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAg                                                            <  1:403605/148‑1 (MQ=255)
                                                                               acaGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGaa                                                          >  1:378812/1‑149 (MQ=255)
                                                                               acaGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGaa                                                          >  8:400513/1‑149 (MQ=255)
                                                                                caGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAat                                                         >  2:363773/1‑149 (MQ=255)
                                                                                 tgACTGCTCTTCCCAACACCGGTTTTACCCATACTGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTg                                                                                                        <  8:139768/100‑1 (MQ=255)
                                                                                 aGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAatt                                                        <  7:400513/149‑1 (MQ=255)
                                                                                     tGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAattatt                                                     >  2:323583/1‑148 (MQ=255)
                                                                                       ctctTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATtcat                                                  <  7:320495/149‑1 (MQ=255)
                                                                                       ctctTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATtcat                                                  <  8:218701/149‑1 (MQ=255)
                                                                                        tctTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATtcat                                                  >  6:200601/1‑148 (MQ=255)
                                                                                          ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGTCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCac                                               >  3:57212/1‑149 (MQ=255)
                                                                                          ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCac                                               >  6:45446/1‑149 (MQ=255)
                                                                                          ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCac                                               >  5:327881/1‑149 (MQ=255)
                                                                                            cccAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATtcatca                                                <  4:118394/146‑1 (MQ=255)
                                                                                            cccAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCacta                                             >  3:198863/1‑149 (MQ=255)
                                                                                               aaCACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTActc                                          >  1:299160/1‑149 (MQ=255)
                                                                                                  acCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTActcctc                                       >  6:336986/1‑149 (MQ=255)
                                                                                                  acCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCGCTGATACCCATTGAGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTACAGAATAATTCATCACTActcct                                        >  4:41629/1‑148 (MQ=255)
                                                                                                   ccGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTTCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCctct                                      <  5:125637/149‑1 (MQ=255)
                                                                                                           aCCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAaat                               <  4:197854/148‑1 (MQ=255)
                                                                                                            cccATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTActcct                                        >  1:382615/1‑138 (MQ=255)
                                                                                                            cccATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTActcct                                        <  2:382615/138‑1 (MQ=255)
                                                                                                             ccATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAatta                                                       >  6:378855/1‑122 (MQ=255)
                                                                                                             ccATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAatta                                                       <  5:378855/122‑1 (MQ=255)
                                                                                                               aTAATGCCGATCACTGGCTCGTAGTTAGTAAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAataata                           >  8:324749/1‑148 (MQ=255)
                                                                                                               aTAATGCCGATCACTGGCTCGTAGTTAGTAAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAataata                           >  8:330417/1‑148 (MQ=255)
                                                                                                               aTAATGACGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAataata                           >  5:275294/1‑148 (MQ=255)
                                                                                                                tAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCac                                               <  8:44584/127‑1 (MQ=255)
                                                                                                                tAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCac                                               >  7:44584/1‑127 (MQ=255)
                                                                                                                tAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAAtaataa                          >  8:271432/1‑148 (MQ=255)
                                                                                                                   tGCCGAGCACTGCCTCGTAGTTAGTAAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATaaaaaa                      <  8:341401/149‑1 (MQ=255)
                                                                                                                    gCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAAc                     >  3:122270/1‑149 (MQ=255)
                                                                                                                    gCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAAc                     >  8:368992/1‑149 (MQ=255)
                                                                                                                     ccGATCACTGGCTCGTAGTTAGTAAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACg                    <  7:330417/149‑1 (MQ=255)
                                                                                                                     ccGATCACTGGCTCGTAGTTAGTAAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACg                    <  7:324749/149‑1 (MQ=255)
                                                                                                                       gATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACgg                   >  4:207468/1‑148 (MQ=255)
                                                                                                                             tGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATc            <  1:363773/149‑1 (MQ=255)
                                                                                                                               gCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCctct                                      >  6:134060/1‑121 (MQ=255)
                                                                                                                               gCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCctct                                      <  5:134060/121‑1 (MQ=255)
                                                                                                                               gCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGt          >  7:318123/1‑149 (MQ=255)
                                                                                                                               gCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTACTCTGAAAAATAATAAAAAACGGTAGAATCGt          <  6:186768/149‑1 (MQ=255)
                                                                                                                                 gcGTAGTTCGTTCACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAAAAATAAAAAACGGTAGAATCGTga        <  2:242495/148‑1 (MQ=255)
                                                                                                                                   gTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTgaga      <  2:378812/149‑1 (MQ=255)
                                                                                                                                   gTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTga        <  8:339697/147‑1 (MQ=255)
                                                                                                                                    tagttagtTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGAt     <  6:327881/149‑1 (MQ=255)
                                                                                                                                     agttagttAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGAtt    >  3:40521/1‑149 (MQ=255)
                                                                                                                                       ttagttaACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTcc  <  5:336986/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                               
AGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCC  >  NC_000913/268725‑269008

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