Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,079,004 G→C 7.7% intergenic (+73/+28) yoeF → / ← yeeX pseudogene, CP4‑44 putative prophage remnant;Phage or Prophage Related/UPF0265 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,079,0040GC 7.7% 273.8 / 10.5 87intergenic (+73/+28)yoeF/yeeXpseudogene, CP4‑44 putative prophage remnant;Phage or Prophage Related/UPF0265 family protein
Reads supporting (aligned to +/- strand):  ref base G (33/47);  new base C (5/2);  total (38/49)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

GTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTC  >  NC_000913/2078857‑2079151
                                                                                                                                                   |                                                                                                                                                   
gTAACCAGTGCAGTTAATAACAATGGAGTTACAGCACTTTTTGCCTTGCTGCATTCAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCgg                                                                                                                                                    >  4:51912/1‑149 (MQ=255)
gTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCgg                                                                                                                                                    <  4:233163/149‑1 (MQ=255)
gTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCgg                                                                                                                                                    <  3:167450/149‑1 (MQ=255)
 tAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCggg                                                                                                                                                   >  2:24856/1‑149 (MQ=255)
       gTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATTAGCAATGAGGGTTGCCGGGCAAccc                                                                                                                                             >  3:195763/1‑149 (MQ=255)
          cAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCa                                                                                                                                          >  8:229834/1‑149 (MQ=255)
          cAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTGAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTTCCCGGCAACCCTCa                                                                                                                                          >  7:401510/1‑149 (MQ=255)
            gTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCAtt                                                                                                                                        <  6:131890/149‑1 (MQ=255)
               acaaaCACTGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGGGCTCGGATGAGAATGATGACGGTGATCCACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGa                                                                                                                                      <  6:87174/145‑1 (MQ=255)
                ataaCAATGGAGTTACAGAACTTTTTTCATTGCTGCATACAGAAGAAGGATGTCGTTAATTTCGATTGGGTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGaa                                                                                                                                     <  5:107989/148‑1 (MQ=255)
                 taaCAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAAt                                                                                                                                    >  4:349252/1‑148 (MQ=255)
                        ggAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAAcc                                                                                                                                              <  5:229655/131‑1 (MQ=255)
                        ggAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAAcc                                                                                                                                              >  6:229655/1‑131 (MQ=255)
                         gAGTGACAGAACTGTTTGCCTTGCTGCATACAGAAGAAGGATGTTGTTAATTTCTAGTGGTTCTCGGAGGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGa                                                                                                                           <  6:287088/149‑1 (MQ=255)
                           gTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAAt                                                                                                                         <  1:212943/149‑1 (MQ=255)
                              aCAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTAc                                                                                                                      >  8:336699/1‑149 (MQ=255)
                                 gAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTActtc                                                                                                                   <  6:140504/149‑1 (MQ=255)
                                  aaCTTTTTGCCTTGCTGCATACAGAAGAGGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTActtcg                                                                                                                  >  3:90861/1‑149 (MQ=255)
                                  aaCTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGTAATTActtcg                                                                                                                  >  5:271258/1‑149 (MQ=255)
                                  aaCTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGGTGCCGGGCAACCCTCATTGAATAAAACGGGAATTActtcg                                                                                                                  <  2:344313/149‑1 (MQ=255)
                                     tttttGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGTGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTActtcgctt                                                                                                               >  1:338465/1‑149 (MQ=255)
                                     tttttGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTActtcgctt                                                                                                               >  4:109272/1‑149 (MQ=255)
                                     tttttGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTActtcgctt                                                                                                               >  1:218316/1‑149 (MQ=255)
                                      ttttGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTActtcgcttc                                                                                                              >  7:124742/1‑149 (MQ=255)
                                       tttGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTActtcgcttc                                                                                                              <  8:81183/148‑1 (MQ=255)
                                        ttGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTActtcgcttcg                                                                                                             <  8:124742/148‑1 (MQ=255)
                                         tGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACttcgcttcgc                                                                                                            <  4:97766/148‑1 (MQ=255)
                                            cTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGtt                                                                                                        <  2:218316/149‑1 (MQ=255)
                                             ttGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGttt                                                                                                       >  5:367588/1‑149 (MQ=255)
                                               gctgcATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGttt                                                                                                       <  1:172702/147‑1 (MQ=255)
                                                  gcATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACttcgcttcgc                                                                                                            <  7:384277/139‑1 (MQ=255)
                                                  gcATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACttcgcttcgc                                                                                                            >  8:384277/1‑139 (MQ=255)
                                                     tACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTc                                                                                               >  6:200461/1‑149 (MQ=255)
                                                     tACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTc                                                                                               <  1:160987/149‑1 (MQ=255)
                                                      aCAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGaattaatt                                                                                                                    >  5:318773/1‑127 (MQ=255)
                                                      aCAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTActt                                                                                                                    <  6:318773/127‑1 (MQ=255)
                                                      aCAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCg                                                                                              >  8:287429/1‑149 (MQ=255)
                                                       cAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGc                                                                                             <  7:229834/149‑1 (MQ=255)
                                                       cAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGc                                                                                             <  6:367588/149‑1 (MQ=255)
                                                        agaagaAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGcc                                                                                            <  3:349252/149‑1 (MQ=255)
                                                        agaagaAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGcc                                                                                            >  7:91057/1‑149 (MQ=255)
                                                                    ttgttAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTc                                                                                  <  8:323139/147‑1 (MQ=255)
                                                                       ttAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTTATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGttt                                                                                                       <  1:342484/123‑1 (MQ=255)
                                                                       ttAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGttt                                                                                                       >  2:342484/1‑123 (MQ=255)
                                                                       ttAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCtt                                                                             >  1:331827/1‑149 (MQ=255)
                                                                         aaTTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCtttt                                                                           >  3:88589/1‑149 (MQ=255)
                                                                          aTTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCtttt                                                                           <  1:350878/148‑1 (MQ=255)
                                                                           tttCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTg                                                                         >  5:204026/1‑149 (MQ=255)
                                                                           tttCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTg                                                                         <  6:361415/149‑1 (MQ=255)
                                                                           tttCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTg                                                                         <  4:90861/149‑1 (MQ=255)
                                                                              cGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGga                                                                       <  7:287429/148‑1 (MQ=255)
                                                                               gATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGgag                                                                      <  3:109272/148‑1 (MQ=255)
                                                                                          ggATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCgcggc                                                            <  3:96533/147‑1 (MQ=255)
                                                                                          ggATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCgcggc                                                            <  4:88589/147‑1 (MQ=255)
                                                                                               agaATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCtt                                                      <  8:91057/148‑1 (MQ=255)
                                                                                                  atgatgACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATaa                                                                                                                                  <  8:155081/69‑1 (MQ=255)
                                                                                                  atgatgACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGga                                                  <  1:24856/149‑1 (MQ=255)
                                                                                                  atgatgACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGga                                                  <  7:322646/149‑1 (MQ=255)
                                                                                                  atgatgAAGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGTCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGga                                                  <  3:246448/149‑1 (MQ=255)
                                                                                                   tgatgaCGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTg                                                    >  1:259538/1‑146 (MQ=255)
                                                                                                   tgatgaCGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTg                                                    >  5:300496/1‑146 (MQ=255)
                                                                                                      tgccGGCGATACACCGATAAGCAGCAATGAGCAAGGAGGGATGCCGGGCAACCCTCATTGAATACAACCGGAAGTACTTCGCTTCGCCGTTTTTCATTTCCCCCATCGCTTGCAGCTTTTTGGAGATATCGCGGCGCTCGTTGGAGAGc                                              <  4:21491/146‑1 (MQ=255)
                                                                                                         cGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCTCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCg                                           <  6:271258/149‑1 (MQ=255)
                                                                                                               tACAAAGATAAGAATCAATCAGCAATGAGGGTTGCCGGGCAACCGTCATTGAATAACTCGGGTATTAGTTCGCTTCGCCGTTTTTCATTTCGCCGATAGCTTTCTGCTTTTTGGAGATATTGCGGCGTTCTTTGGAGAGCTCGGCAtt                                      <  7:349495/148‑1 (MQ=255)
                                                                                                                acacAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCAtttt                                    >  6:238921/1‑149 (MQ=255)
                                                                                                                 cacaGATAATAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTg                                                    <  6:79332/132‑1 (MQ=255)
                                                                                                                 cacaGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCAtttt                                    <  6:204026/148‑1 (MQ=255)
                                                                                                                 cacaGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGATTTTTGGAGATATCGCGGCGTTCTTTg                                                    >  5:79332/1‑132 (MQ=255)
                                                                                                                  acaGCTAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAACCGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCAttttt                                   >  3:81418/1‑148 (MQ=255)
                                                                                                                  acaGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCACCCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCAttttt                                   >  5:190188/1‑148 (MQ=255)
                                                                                                                    agctAAGAATCAATGAGCAAGGAGGGTTGCCGGGCAACCCTCATTGAATAAAACTGTAATTCATTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTtgat                                <  6:32835/146‑1 (MQ=255)
                                                                                                                      aTAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTtgatga                              <  3:51912/149‑1 (MQ=255)
                                                                                                                       tAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTgatgat                             <  5:200461/149‑1 (MQ=255)
                                                                                                                          gAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGt                           >  5:178827/1‑148 (MQ=255)
                                                                                                                          gAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCACTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGc                                                                                             >  8:232435/1‑82 (MQ=255)
                                                                                                                          gAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCACTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGc                                                                                             <  7:232435/82‑1 (MQ=255)
                                                                                                                           aaTCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTaa                         >  1:205234/1‑149 (MQ=255)
                                                                                                                             tCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATc                       <  5:238921/149‑1 (MQ=255)
                                                                                                                               aatgagcaatgagGGTTGCCGGGCACCCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGt                     >  7:60461/1‑149 (MQ=255)
                                                                                                                               aatgagcaatgagGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCg                      <  2:338465/148‑1 (MQ=255)
                                                                                                                                ttattcaatgagGGTTGCCCGGCAACCCTCATTGAATAAAACGGGATTTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGt                     >  7:369523/6‑148 (MQ=255)
                                                                                                                                ttattcaatgagGGTTGCCCGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGt                     >  3:180029/6‑148 (MQ=255)
                                                                                                                                    gcaatgagGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAAcg                >  5:103358/1‑149 (MQ=255)
                                                                                                                                     caatgagGGTTGCCCGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTc                                                                                  >  7:37911/1‑82 (MQ=255)
                                                                                                                                     caatgagGGTTGCCCGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTc                                                                                  <  8:37911/82‑1 (MQ=255)
                                                                                                                                       atgagGGTTGCCCGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGa             <  5:377293/149‑1 (MQ=255)
                                                                                                                                        tgagGGTTGCCCGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGAt            >  8:251688/1‑149 (MQ=255)
                                                                                                                                           gggTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCtt         >  5:275091/1‑149 (MQ=255)
                                                                                                                                               tGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATa     <  6:178827/149‑1 (MQ=255)
                                                                                                                                                gCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATa     <  7:336699/148‑1 (MQ=255)
                                                                                                                                                  cGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTc  <  6:300496/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTC  >  NC_000913/2078857‑2079151

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