Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 569,757 G→A 5.4% A286T (GCA→ACA)  ybcK → DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913569,7570GA 5.4% 239.4 / 2.1 73A286T (GCA→ACA) ybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base G (36/33);  new base A (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

AATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAA  >  NC_000913/569617‑569904
                                                                                                                                            |                                                                                                                                                   
aaTTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCtttt                                                                                                                                             <  2:241620/149‑1 (MQ=255)
 attattAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTg                                                                                                                                            <  4:180878/149‑1 (MQ=255)
 attattAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTg                                                                                                                                            <  4:24016/149‑1 (MQ=255)
   tattaGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGt                                                                                                                                          <  3:166419/149‑1 (MQ=255)
        gCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAg                                                                                                                                     >  2:481376/1‑149 (MQ=255)
             tAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAAta                                                                                                                                <  1:442303/149‑1 (MQ=255)
                 gCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAAtagtag                                                                                                                            >  1:249070/1‑149 (MQ=255)
                       aTAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGaa                                                                                                                      >  4:88909/1‑149 (MQ=255)
                         aGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATc                                                                                                                    >  3:64341/1‑149 (MQ=255)
                         aGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTCAAATCGCTCGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACGGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATc                                                                                                                    <  4:350033/149‑1 (MQ=255)
                          ggTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATcc                                                                                                                   <  2:374926/149‑1 (MQ=255)
                          ggTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATcc                                                                                                                   >  1:183060/1‑149 (MQ=255)
                                     aCCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCtaataa                                                                                                        <  3:313116/149‑1 (MQ=255)
                                       cTTCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAAt                                                                                                      >  1:256828/1‑149 (MQ=255)
                                         tCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCt                                                                                                    >  2:273381/1‑149 (MQ=255)
                                           atatCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTAc                                                                                                  <  2:581407/149‑1 (MQ=255)
                                             atCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCaa                                                                                                                                  >  1:492772/1‑115 (MQ=255)
                                             atCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCaa                                                                                                                                  <  2:492772/115‑1 (MQ=255)
                                             atCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAAta                                                                                                                                <  1:89674/117‑1 (MQ=255)
                                             atCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAAta                                                                                                                                >  2:89674/1‑117 (MQ=255)
                                                 tGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGaa                                                                                            <  2:183060/149‑1 (MQ=255)
                                                  gCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAAc                                                                                           <  3:197081/149‑1 (MQ=255)
                                                    aaGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAg                                                                                         >  1:496206/1‑149 (MQ=255)
                                                    aaGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAg                                                                                         >  4:323505/1‑149 (MQ=255)
                                                          gAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGa                                                                                   >  1:13238/1‑149 (MQ=255)
                                                                 aTAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGa                                                                            >  4:318488/1‑149 (MQ=255)
                                                                  tAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAg                                                                           <  2:13238/149‑1 (MQ=255)
                                                                         aaTCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTtgtg                                                                    <  2:343527/149‑1 (MQ=255)
                                                                         aaTCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTtgtg                                                                    >  1:428063/1‑149 (MQ=255)
                                                                         aaTCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTACGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTtgtg                                                                    >  1:278091/1‑149 (MQ=255)
                                                                          aTCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTgg                                                                   >  2:362459/1‑149 (MQ=255)
                                                                              cTGGCTATTAGCACAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCGATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAAt                                                               <  2:23604/149‑1 (MQ=255)
                                                                                  cTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTACGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCa                                                           <  2:278091/149‑1 (MQ=255)
                                                                                   tattatCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCAt                                                          >  3:183910/1‑149 (MQ=255)
                                                                                    attatCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATg                                                         <  3:547561/149‑1 (MQ=255)
                                                                                    attatCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATg                                                         <  2:87282/149‑1 (MQ=255)
                                                                                        tCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTg                                                     >  4:98740/1‑149 (MQ=255)
                                                                                          ccAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGACACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGtt                                                   >  3:23997/1‑149 (MQ=255)
                                                                                            agagTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAAc                                                                                           >  4:312770/1‑107 (MQ=255)
                                                                                            agagTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAAc                                                                                           <  3:312770/107‑1 (MQ=255)
                                                                                              agTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATg                                               <  3:88909/149‑1 (MQ=255)
                                                                                               gTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGc                                              >  3:232052/1‑149 (MQ=255)
                                                                                               gTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGc                                              <  3:323505/149‑1 (MQ=255)
                                                                                                  atatCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGt                                           <  2:143194/149‑1 (MQ=255)
                                                                                                  atatCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGt                                           >  1:143194/1‑149 (MQ=255)
                                                                                                   tatCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTa                                          >  1:546184/1‑149 (MQ=255)
                                                                                                     tCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACtt                                                                                                <  2:263603/93‑1 (MQ=255)
                                                                                                     tCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACtt                                                                                                >  1:263603/1‑93 (MQ=255)
                                                                                                     tCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGACACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTtgtg                                                                    <  3:322016/121‑1 (MQ=255)
                                                                                                     tCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGACACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTtgtg                                                                    >  4:322016/1‑121 (MQ=255)
                                                                                                        gatgatTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTAt                                                                                                                <  4:364826/74‑1 (MQ=255)
                                                                                                        gatgatTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTAt                                                                                                                >  3:364826/1‑74 (MQ=255)
                                                                                                        gatgatTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAgg                                                                          >  4:68862/1‑112 (MQ=255)
                                                                                                        gatgatTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAgg                                                                          <  3:68862/112‑1 (MQ=255)
                                                                                                        gatgatTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTgg                                     >  3:289908/1‑149 (MQ=255)
                                                                                                         atgatTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGa                                    >  1:535575/1‑149 (MQ=255)
                                                                                                          tgatTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGaa                                   >  3:514380/1‑149 (MQ=255)
                                                                                                              ttGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGttt                               <  1:558355/149‑1 (MQ=255)
                                                                                                               tGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTg                              >  4:560329/1‑149 (MQ=255)
                                                                                                                 ttttACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCa                            <  1:273381/149‑1 (MQ=255)
                                                                                                                 ttttACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCa                            >  4:438746/1‑149 (MQ=255)
                                                                                                                    tACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATgg                         <  4:289908/149‑1 (MQ=255)
                                                                                                                        cTGTACAGGAAATTCGGTTGATACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACCGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCtat                     <  4:23997/149‑1 (MQ=255)
                                                                                                                             cAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGt                >  2:58774/1‑149 (MQ=255)
                                                                                                                              aGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGtt               <  2:89391/149‑1 (MQ=255)
                                                                                                                              aGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGtt               >  4:73958/1‑149 (MQ=255)
                                                                                                                               ggAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGttt              >  2:197841/1‑149 (MQ=255)
                                                                                                                               ggAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGttt              >  4:447813/1‑149 (MQ=255)
                                                                                                                                  aaTTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTc           >  4:217780/1‑149 (MQ=255)
                                                                                                                                      cGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGAt       <  1:30169/149‑1 (MQ=255)
                                                                                                                                      cGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGAt       <  1:197841/149‑1 (MQ=255)
                                                                                                                                        gTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATga     >  1:560904/1‑149 (MQ=255)
                                                                                                                                          tGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATgaga   <  1:481376/149‑1 (MQ=255)
                                                                                                                                           ggCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGagaa  <  4:514380/149‑1 (MQ=255)
                                                                                                                                           ggCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGagaa  >  1:78166/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
AATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGAAATCGCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAA  >  NC_000913/569617‑569904

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