Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 292,302 C→G 7.1% intergenic (+71/+20) yagJ → / ← yagK CP4‑6 prophage; uncharacterized protein;Phage or Prophage Related/CP4‑6 prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913292,3020CG 7.1% 345.1 / 14.8 99intergenic (+71/+20)yagJ/yagKCP4‑6 prophage; uncharacterized protein;Phage or Prophage Related/CP4‑6 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base C (51/41);  new base G (3/4);  total (54/45)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

ATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGT  >  NC_000913/292155‑292450
                                                                                                                                                   |                                                                                                                                                    
aTGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACt                                                                                                                                                     <  2:62614/149‑1 (MQ=255)
  gtgCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGCAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTga                                                                                                                                                   <  1:605984/149‑1 (MQ=255)
   tgCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTgag                                                                                                                                                  >  2:4210/1‑149 (MQ=255)
      aGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAg                                                                                                                                               <  1:525151/149‑1 (MQ=255)
       ggATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAgg                                                                                                                                              <  3:458716/149‑1 (MQ=255)
        gATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGc                                                                                                                                             <  1:181707/149‑1 (MQ=255)
        gATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGc                                                                                                                                             >  4:628124/1‑149 (MQ=255)
         aTGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcc                                                                                                                                            <  4:664465/149‑1 (MQ=255)
         aTGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcc                                                                                                                                            <  3:280302/149‑1 (MQ=255)
         aTGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGcc                                                                                                                                            <  3:308551/149‑1 (MQ=255)
                 gCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGa                                                                                                                                    <  1:257332/149‑1 (MQ=255)
                  cTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGAc                                                                                                                                   <  1:293739/149‑1 (MQ=255)
                  cTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGAc                                                                                                                                   >  3:370215/1‑149 (MQ=255)
                       gCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACa                                                                                                                              <  1:773235/149‑1 (MQ=255)
                       gCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACa                                                                                                                              <  4:499843/149‑1 (MQ=255)
                         tGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCCGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGt                                                                                                                            >  1:570929/1‑149 (MQ=255)
                            aaaTGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCg                                                                                                                         >  4:674909/1‑149 (MQ=255)
                                  aGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGc                                                                                                                                             >  1:411713/1‑123 (MQ=255)
                                  aGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGc                                                                                                                                             <  2:411713/123‑1 (MQ=255)
                                  aGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTgg                                                                                                                   <  4:155757/149‑1 (MQ=255)
                                  aGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTgg                                                                                                                   <  4:47395/149‑1 (MQ=255)
                                  aGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTgg                                                                                                                   >  2:198758/1‑149 (MQ=255)
                                   gAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGc                                                                                                                  >  4:433122/1‑149 (MQ=255)
                                     aaCGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTa                                                                                                                >  2:79585/1‑149 (MQ=255)
                                       cGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACa                                                                                                              >  1:583170/1‑149 (MQ=255)
                                         ttATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGc                                                                                                            >  4:394661/1‑149 (MQ=255)
                                           aTTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCg                                                                                                          >  2:564726/1‑149 (MQ=255)
                                            tttAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGa                                                                                                         >  3:263314/1‑149 (MQ=255)
                                             ttAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGaa                                                                                                        <  4:259394/149‑1 (MQ=255)
                                             ttAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGaa                                                                                                        <  3:507158/149‑1 (MQ=255)
                                                gtagCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGtt                                                                                                     >  2:69028/1‑149 (MQ=255)
                                                  agCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTAc                                                                                                   >  3:747817/1‑149 (MQ=255)
                                                   gCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACg                                                                                                  <  3:6272/149‑1 (MQ=255)
                                                   gCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACg                                                                                                  <  1:141573/149‑1 (MQ=255)
                                                    cGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACgg                                                                                                 <  4:481934/149‑1 (MQ=255)
                                                      gACTTGCTGATGACTATCAGTAATCATTAGAAAATTCCGTAACGAAACTGCATTTTAGTCATTTTTTGACACAAAACAAAGAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACATATCGACTTCGCTACCGCCGAAGTTACGGtc                                                                                               <  4:469854/149‑1 (MQ=255)
                                                         ttGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTcc                                                                                            >  2:235360/1‑149 (MQ=255)
                                                         ttGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTcc                                                                                            >  4:40110/1‑149 (MQ=255)
                                                           gCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGCCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCtt                                                                                          >  4:77883/1‑149 (MQ=255)
                                                               atgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCAcc                                                                                      >  3:212554/1‑149 (MQ=255)
                                                                    tATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAgg                                                                                                                                              >  2:179558/1‑88 (MQ=255)
                                                                    tATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAgg                                                                                                                                              <  1:179558/88‑1 (MQ=255)
                                                                    tATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAaca                                                                                 >  3:754806/1‑149 (MQ=255)
                                                                    tATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAaca                                                                                 >  1:492900/1‑149 (MQ=255)
                                                                        aGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACttt                                                                             <  2:179317/149‑1 (MQ=255)
                                                                              cTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTc                                                                       <  4:370215/149‑1 (MQ=255)
                                                                               tataGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCa                                                                      >  2:636178/1‑149 (MQ=255)
                                                                                 taGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAAc                                                                    >  3:272396/1‑149 (MQ=255)
                                                                                   gAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACtt                                                                  <  3:394661/149‑1 (MQ=255)
                                                                                    aaaaTTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACttt                                                                 <  3:116851/149‑1 (MQ=255)
                                                                                     aaaTTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTg                                                                <  2:590155/149‑1 (MQ=255)
                                                                                      aaTTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTgg                                                               >  4:397741/1‑149 (MQ=255)
                                                                                            gTAACGAAACTGCATTTTAGTACTTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGt                                                         >  2:200053/1‑149 (MQ=255)
                                                                                            gTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGt                                                         >  1:172761/1‑149 (MQ=255)
                                                                                             tAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTa                                                        <  1:749415/149‑1 (MQ=255)
                                                                                              aaCGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAgg                                                                                                                                              >  3:155965/1‑62 (MQ=255)
                                                                                              aaCGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAgg                                                                                                                                              <  4:155965/62‑1 (MQ=255)
                                                                                              aaCGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAg                                                       >  3:11950/1‑149 (MQ=255)
                                                                                               aCGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGt                                                      <  4:89656/149‑1 (MQ=255)
                                                                                                     cTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTc                                                >  2:676797/1‑149 (MQ=255)
                                                                                                       gCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCgg                                              <  1:564726/149‑1 (MQ=255)
                                                                                                       gCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCgg                                              <  3:628124/149‑1 (MQ=255)
                                                                                                         aTTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGt                                            <  2:582082/149‑1 (MQ=255)
                                                                                                         aTTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGt                                            >  4:700161/1‑149 (MQ=255)
                                                                                                          ttttAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGtt                                           >  1:167070/1‑149 (MQ=255)
                                                                                                                tAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgt                                     >  3:67118/1‑148 (MQ=255)
                                                                                                                tAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgc                                     >  2:766382/1‑149 (MQ=255)
                                                                                                                  aTTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgcgc                                   >  3:9541/1‑149 (MQ=255)
                                                                                                                  aTTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgcgc                                   >  2:316775/1‑149 (MQ=255)
                                                                                                                    tttttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTAAGGTCTCCTTCACCAAACACCTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCt                                 <  2:105391/149‑1 (MQ=255)
                                                                                                                      tttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTgg                               <  4:754806/149‑1 (MQ=255)
                                                                                                                              aaTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctg                       >  3:711912/1‑149 (MQ=255)
                                                                                                                              aaTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCCTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctg                       <  1:69028/149‑1 (MQ=255)
                                                                                                                               aTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgc                      >  2:278589/1‑149 (MQ=255)
                                                                                                                               aTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgc                      <  2:172761/149‑1 (MQ=255)
                                                                                                                                tACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgct                     >  1:309774/1‑149 (MQ=255)
                                                                                                                                tACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgct                     >  1:607984/1‑149 (MQ=255)
                                                                                                                                tACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgct                     >  2:469771/1‑149 (MQ=255)
                                                                                                                                tACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgct                     >  4:167492/1‑149 (MQ=255)
                                                                                                                                  cAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTgg                   <  2:492900/149‑1 (MQ=255)
                                                                                                                                   aaaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGa                  >  1:528621/1‑149 (MQ=255)
                                                                                                                                   aaaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGa                  >  2:93296/1‑149 (MQ=255)
                                                                                                                                    aaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGaa                 <  3:361821/149‑1 (MQ=255)
                                                                                                                                    aaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGaa                 <  2:607984/149‑1 (MQ=255)
                                                                                                                                    aaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGaa                 <  1:64604/149‑1 (MQ=255)
                                                                                                                                    aaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGaa                 <  1:636178/149‑1 (MQ=255)
                                                                                                                                      tAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGt               >  4:327605/1‑149 (MQ=255)
                                                                                                                                        aGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCa             >  2:713964/1‑149 (MQ=255)
                                                                                                                                        aGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTAACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCa             >  4:619808/1‑149 (MQ=255)
                                                                                                                                          gCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGt           >  3:479093/1‑149 (MQ=255)
                                                                                                                                           ccTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTa          >  1:69952/1‑149 (MQ=255)
                                                                                                                                            cTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTAc         >  1:715527/1‑149 (MQ=255)
                                                                                                                                             ttCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACt        <  4:272396/149‑1 (MQ=255)
                                                                                                                                               ctcACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTAt      >  3:200465/1‑149 (MQ=255)
                                                                                                                                                  aCTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTCACAGCGCCTGGTAGATctgc                      >  3:665046/1‑130 (MQ=255)
                                                                                                                                                  aCTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTg                                <  2:579203/120‑1 (MQ=255)
                                                                                                                                                  aCTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTg                                >  1:579203/1‑120 (MQ=255)
                                                                                                                                                  aCTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgc                      <  4:665046/130‑1 (MQ=255)
                                                                                                                                                  aCTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTgg   >  1:639238/1‑149 (MQ=255)
                                                                                                                                                   cTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGgt  <  4:365570/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGT  >  NC_000913/292155‑292450

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