Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 570,520 T→C intergenic (+92/‑373) ybcK → / → ybcL DLP12 prophage; putative recombinase/DLP12 prophage; periplasmic protein YbcL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913570,5200TC98.8% 284.2 / ‑4.5 82intergenic (+92/‑373)ybcK/ybcLDLP12 prophage; putative recombinase/DLP12 prophage; periplasmic protein YbcL
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (42/39);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

GATAAGTGCTCTCAAATATATGGTTGATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGTTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAA  >  NC_000913/570383‑570655
                                                                                                                                         |                                                                                                                                       
gATAAGTGCTCTCAAATATATGGTTGATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACtttgctt                                                                                                                                       >  2:934720/1‑140 (MQ=255)
             aaaTATATGGTTGATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCg                                                                                                                          <  2:125658/140‑1 (MQ=255)
             aaaTATATGGTTGATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCg                                                                                                                          <  2:245483/140‑1 (MQ=255)
                     ggTTGATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGt                                                                                                                  >  2:1106517/1‑140 (MQ=255)
                     ggTTGATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGt                                                                                                                  >  2:1015829/1‑140 (MQ=255)
                      gTTGATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGTTTGCTCCCACAAGCGAGTTTTGTa                                                                                                                 <  2:243978/140‑1 (MQ=255)
                      gTTGATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTa                                                                                                                 <  2:1077569/140‑1 (MQ=255)
                       ttGATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTAc                                                                                                                >  2:139926/1‑140 (MQ=255)
                       ttGATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTAc                                                                                                                >  2:644726/1‑140 (MQ=255)
                        tGATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACg                                                                                                               >  1:365727/1‑140 (MQ=255)
                         gATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTTCCACAAGCGAGTTTTGTACgg                                                                                                              <  2:694176/140‑1 (MQ=255)
                         gATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACgg                                                                                                              <  1:1015829/140‑1 (MQ=255)
                         gATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACgg                                                                                                              <  2:222169/140‑1 (MQ=255)
                         gATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACgg                                                                                                              <  1:644726/140‑1 (MQ=255)
                         gATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACgg                                                                                                              <  2:260088/140‑1 (MQ=255)
                         gATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACgg                                                                                                              >  1:222169/1‑140 (MQ=255)
                         gATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACgg                                                                                                              <  1:10686/140‑1 (MQ=255)
                         gATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACgg                                                                                                              <  1:139926/140‑1 (MQ=255)
                         gATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACgg                                                                                                              <  2:962364/140‑1 (MQ=255)
                          aTGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGtt                                                                                                                      >  2:571252/1‑131 (MQ=255)
                          aTGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGtt                                                                                                                      <  1:571252/131‑1 (MQ=255)
                          aTGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGc                                                                                                             <  2:380044/140‑1 (MQ=255)
                          aTGGTGAGATTTATTTATAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTg                                                                                                                   >  2:151911/1‑134 (MQ=255)
                          aTGGTGAGATTTATTTATAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTg                                                                                                                   <  1:151911/134‑1 (MQ=255)
                            ggTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTg                                                                                                           <  2:494819/140‑1 (MQ=255)
                             gTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGt                                                                                                          <  2:716025/140‑1 (MQ=255)
                              tGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTa                                                                                                         >  1:173377/1‑140 (MQ=255)
                               gagaTTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTAt                                                                                                        >  1:447617/1‑140 (MQ=255)
                                agaTTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTAtt                                                                                                       <  2:980453/140‑1 (MQ=255)
                                 gaTTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTg                                                                                                      >  1:713443/1‑140 (MQ=255)
                                            aaaTAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAAt                                                                                           <  1:375725/140‑1 (MQ=255)
                                            aaaTAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAAt                                                                                           >  2:371722/1‑140 (MQ=255)
                                            aaaTAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAAt                                                                                           >  1:665775/1‑140 (MQ=255)
                                            aaaTAATGATCTCGAATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAAt                                                                                           >  1:350228/1‑140 (MQ=255)
                                               taatGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAg                                                                                        <  2:665775/140‑1 (MQ=255)
                                                 atGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGgt                                                                                                  >  1:586456/1‑128 (MQ=255)
                                                 atGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGgt                                                                                                  <  2:586456/128‑1 (MQ=255)
                                                  tGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCta                                                                                     >  2:478353/1‑140 (MQ=255)
                                                    aTCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCtata                                                                                   <  2:713443/140‑1 (MQ=255)
                                                    aTCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCtata                                                                                   <  2:525162/140‑1 (MQ=255)
                                                    aTCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCtata                                                                                   >  1:290123/1‑140 (MQ=255)
                                                       tcGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACaa                                                                                >  1:1016581/1‑140 (MQ=255)
                                                       tcGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACaa                                                                                >  1:50800/1‑140 (MQ=255)
                                                               tAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATc                                                                        >  1:922536/1‑140 (MQ=255)
                                                                aaGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCa                                                                       <  1:371722/140‑1 (MQ=255)
                                                                  gTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCAtt                                                                     <  1:1039151/140‑1 (MQ=255)
                                                                     aTGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGt                                                                  >  1:294067/1‑140 (MQ=255)
                                                                      tGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGtt                                                                 >  1:33986/1‑140 (MQ=255)
                                                                       gCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAAt                                                                         >  2:763339/1‑131 (MQ=255)
                                                                       gCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAAt                                                                         <  1:763339/131‑1 (MQ=255)
                                                                           tGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAAt                                                                         >  1:566878/1‑127 (MQ=255)
                                                                           tGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAAt                                                                         <  2:566878/127‑1 (MQ=255)
                                                                                aTAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAAt                                                                         <  1:175423/122‑1 (MQ=255)
                                                                                aTAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAAc                                                       >  1:507053/1‑140 (MQ=255)
                                                                                aTAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGCATTGGGGTAGTAAATGAGCTATACAATCTTAAt                                                                         >  2:175423/1‑122 (MQ=255)
                                                                                 tAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAAct                                                      >  2:1103438/1‑140 (MQ=255)
                                                                                     gTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTg                                                  <  2:507053/140‑1 (MQ=255)
                                                                                      tGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTgg                                                 >  2:447069/1‑140 (MQ=255)
                                                                                       gCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGt                                                >  2:192873/1‑140 (MQ=255)
                                                                                       gCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGt                                                >  1:503492/1‑140 (MQ=255)
                                                                                        cAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTc                                               >  2:534219/1‑140 (MQ=255)
                                                                                           gACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCa                                            <  2:173377/140‑1 (MQ=255)
                                                                                           gACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCa                                            >  2:310033/1‑140 (MQ=255)
                                                                                           gACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCa                                            <  1:447069/140‑1 (MQ=255)
                                                                                                     attaTCTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATa                                  <  1:192873/140‑1 (MQ=255)
                                                                                                          cTTTGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGaa                             >  2:776704/1‑140 (MQ=255)
                                                                                                             tGACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAAt                          >  2:607655/1‑140 (MQ=255)
                                                                                                              gACGCAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCCTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATa                         <  1:934720/140‑1 (MQ=255)
                                                                                                              gACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATa                         >  2:832706/1‑140 (MQ=255)
                                                                                                              gACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATa                         >  1:970461/1‑140 (MQ=255)
                                                                                                              gACGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATa                         <  2:1016581/140‑1 (MQ=255)
                                                                                                               aCGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCaaa                                    <  1:1028109/128‑1 (MQ=255)
                                                                                                               aCGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCaaa                                    >  2:1028109/1‑128 (MQ=255)
                                                                                                               aCGAAACTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATAc                        <  1:1106517/140‑1 (MQ=255)
                                                                                                                    aCTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGTGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGAc                   <  2:260912/140‑1 (MQ=255)
                                                                                                                    aCTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCaa                                     >  1:799850/1‑122 (MQ=255)
                                                                                                                    aCTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGAc                   <  1:1103438/140‑1 (MQ=255)
                                                                                                                    aCTTAATGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCCGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCaa                                     <  2:799850/122‑1 (MQ=255)
                                                                                                                        aaTGGGTAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAat               >  1:333822/1‑140 (MQ=255)
                                                                                                                              tAATTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATtattat         <  2:447617/140‑1 (MQ=255)
                                                                                                                               aaTTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATg        >  2:1076797/1‑140 (MQ=255)
                                                                                                                                aTTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGAc                   <  1:347451/128‑1 (MQ=255)
                                                                                                                                aTTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGAc                   >  2:347451/1‑128 (MQ=255)
                                                                                                                                aTTACTTTGCTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGa       >  1:369359/1‑140 (MQ=255)
                                                                                                                                     tttgcttgCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCaaa  <  2:356143/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
GATAAGTGCTCTCAAATATATGGTTGATGGTGAGATTTATTTCTAAATAATGATCTCGGATTTTAAGTTATGCTATGGTGATAAAGTGCAAGACAGAATTAATTATCTTTGACGAAACTTAATGGGTAATTACTTTGTTTGCTCCCACAAGCGAGTTTTGTACGGCTGTATTGGGGTAGTAAATGAGCTATACAATCTTAATCATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAA  >  NC_000913/570383‑570655

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