Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,429,128 A→G 11.4% K27R (AAA→AGA)  stfR → Rac prophage; putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,1280AG11.4% 239.7 / 19.4 88K27R (AAA→AGA) stfRRac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base A (38/40);  new base G (4/6);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

GATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACC  >  NC_000913/1428983‑1429276
                                                                                                                                                 |                                                                                                                                                    
gATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACAGGTACAGAACTGCACAATCCAGCAGAAAGCAAAACg                                                                                                                                                   >  1:36532/1‑149 (MQ=255)
 aTTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACg                                                                                                                                                   >  6:132791/1‑148 (MQ=255)
 aTTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACg                                                                                                                                                   <  2:65613/148‑1 (MQ=255)
     gCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTaa                                                                                                                                                >  1:176766/1‑147 (MQ=255)
        aGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGca                                                                                                                                           >  7:342172/1‑149 (MQ=255)
          gTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGcac                                                                                                                                          >  3:17006/1‑148 (MQ=255)
              cacaGAGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGAACCACgg                                                                                                                                     >  3:21314/1‑149 (MQ=255)
                   tgtTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACggtggtg                                                                                                                                <  1:259696/149‑1 (MQ=255)
                     ttttGACAGCCCGGCGGTTGAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggtggt                                                                                                                              <  4:21314/146‑1 (MQ=255)
                     ttATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggtgg                                                                                                                               <  1:277857/148‑1 (MQ=255)
                     ttATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggtgg                                                                                                                               <  6:131867/148‑1 (MQ=255)
                     ttATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggtgg                                                                                                                               <  4:155731/148‑1 (MQ=255)
                        tGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGaa                                                                                                                           >  4:108215/1‑149 (MQ=255)
                         gACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAac                                                                                                                          >  5:289743/1‑149 (MQ=255)
                          aCAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAac                                                                                                                          <  8:413364/148‑1 (MQ=255)
                           cAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAacac                                                                                                                        >  8:103525/1‑149 (MQ=255)
                             gCCCGCCGGAACCGGTGGGCTTTTTTGTGTGGTGAATATGGCAGTAAAGATTTCAGGAGTCATGAAAGACGGCACAGTAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTTGTGAACACgg                                                                                                                      >  6:57356/1‑148 (MQ=255)
                                                 ttttttGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGAt                                                                                                  <  4:64991/149‑1 (MQ=255)
                                                   ttttGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATg                                                                                                 >  7:324970/1‑148 (MQ=255)
                                                       gtgGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGcc                                                                                            >  5:139984/1‑149 (MQ=255)
                                                        tgGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGcc                                                                                            <  3:322371/148‑1 (MQ=255)
                                                        tgGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCg                                                                                           <  2:176766/149‑1 (MQ=255)
                                                          gggTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTAAAGAACTGAACAATCCAGCTGAAAGCAAAAAGTAACAGCACCACGGTGGTGGTGACCACGCTGCCCTCAGAAAATCCGGATGAAGCCggg                                                                                         >  5:207159/1‑149 (MQ=255)
                                                           ggTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGc                                                                                        <  5:114439/149‑1 (MQ=255)
                                                             tGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGt                                                                                      <  5:132791/149‑1 (MQ=255)
                                                               aaTATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGtt                                                                                     <  3:378080/148‑1 (MQ=255)
                                                                atatGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTAc                                                                                   <  3:242169/149‑1 (MQ=255)
                                                                   tGGCAGTAAAGATTTCAGGTGTTCTGAAAGACGTCACAGGAAAACCGGTACTGAACTGGACAATCCAGCTGAAACGAAAACGTAACAGCACCACGGTGGTTGTGAACACGCTGTCCTCAGAAATTCCGGATGAAGCCGGGCGTTACAg                                                                                 <  3:287775/148‑1 (MQ=255)
                                                                     gCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCa                                                                               <  1:357244/148‑1 (MQ=255)
                                                                      cAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATg                                                                             <  4:50403/149‑1 (MQ=255)
                                                                       aGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAAc                                                                                                                                               <  5:253893/82‑1 (MQ=255)
                                                                       aGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAAc                                                                                                                                               >  6:253893/1‑82 (MQ=255)
                                                                       aGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCAAAATCAAGCTGAAAGCAAAACGTAACAGCACCAAGGTGGTGGTGAAAAAGCTGGCCTCATAAAATCCTTATGAAGACGGTCGTTACAGCATg                                                                             >  1:403860/1‑148 (MQ=255)
                                                                        gTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGa                                                                           >  4:131771/1‑149 (MQ=255)
                                                                        gTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGa                                                                           >  4:211108/1‑149 (MQ=255)
                                                                         tAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGAc                                                                          >  2:77119/1‑149 (MQ=255)
                                                                           aaGTTTTCAGGTGTACTGAAAGACGTCACAGGAAAACCGGTACATAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGATCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGt                                                                        <  6:207159/149‑1 (MQ=255)
                                                                            aGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTTCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTTGCCTCAGAAAATCCGGATTAAGCCGGGCGTTACAGCATGGACGtt                                                                       <  7:242040/149‑1 (MQ=255)
                                                                            aGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCgg                                                                                          <  1:309154/130‑1 (MQ=255)
                                                                            aGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCgg                                                                                          >  2:309154/1‑130 (MQ=255)
                                                                            aGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGtt                                                                       <  7:36967/149‑1 (MQ=255)
                                                                            aGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGtt                                                                       <  6:139984/149‑1 (MQ=255)
                                                                            aGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGtt                                                                       <  3:211108/149‑1 (MQ=255)
                                                                            aGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGtt                                                                       <  1:196759/149‑1 (MQ=255)
                                                                                 tCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCTTTACAGCATGGACGTTGAGt                                                                   <  5:42562/148‑1 (MQ=255)
                                                                                 tCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACAGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTa                                                                  >  7:245071/1‑149 (MQ=255)
                                                                                   aGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACg                                                                <  5:328771/149‑1 (MQ=255)
                                                                                      tgtACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGt                                                              >  3:311530/1‑148 (MQ=255)
                                                                                      agtccTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTc                                                             <  5:42843/145‑1 (MQ=255)
                                                                                       gtccTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGGGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCa                                                            <  8:342172/146‑1 (MQ=255)
                                                                                        tACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCa                                                            <  3:131771/148‑1 (MQ=255)
                                                                                        tACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAACCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGTGGGTTAGAGCATGGACGTTGAGTACGGTCAg                                                           >  4:264369/1‑149 (MQ=255)
                                                                                         aCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCCCGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGt                                                              >  5:301044/1‑145 (MQ=255)
                                                                                         aCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGt                                                              <  6:301044/145‑1 (MQ=255)
                                                                                          cTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTa                                                         <  3:108215/149‑1 (MQ=255)
                                                                                           tGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTa                                                         >  4:76680/1‑148 (MQ=255)
                                                                                           tGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTa                                                         >  6:204354/1‑148 (MQ=255)
                                                                                           tGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTa                                                         >  4:87618/1‑148 (MQ=255)
                                                                                           tGAAAGACGGCACAGGAAAACAGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTAc                                                        >  6:287203/1‑149 (MQ=255)
                                                                                             aaaGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACa                                                       <  8:324970/148‑1 (MQ=255)
                                                                                              aaGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGtg                                                                       >  7:378086/1‑130 (MQ=255)
                                                                                              aaGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGtg                                                                       <  8:378086/131‑2 (MQ=255)
                                                                                               aGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtg                                                    <  4:17006/149‑3 (MQ=255)
                                                                                                      acaGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCt                                             <  3:360390/149‑1 (MQ=255)
                                                                                                      acaGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTTTCTg                                            >  2:187373/1‑149 (MQ=255)
                                                                                                      acaGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGCTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGGTATTCt                                             >  4:296502/1‑149 (MQ=255)
                                                                                                        aGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTAt                                                >  5:342721/1‑144 (MQ=255)
                                                                                                        aGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTAt                                                <  6:342721/144‑1 (MQ=255)
                                                                                                         ggAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACGACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGtt                                          <  3:87618/149‑1 (MQ=255)
                                                                                                         ggAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGtt                                          >  1:105654/1‑149 (MQ=255)
                                                                                                                cGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTtggtgg                                     >  6:39030/1‑147 (MQ=255)
                                                                                                                cGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGaa                                   <  4:234807/149‑1 (MQ=255)
                                                                                                                      aGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTc                             <  4:311530/149‑1 (MQ=255)
                                                                                                                        aaCTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTccc                           >  4:181589/1‑149 (MQ=255)
                                                                                                                         aCTGAACAATCCAGCTGAAAGCAAAACGTAACATCACCACGGTGGTGTTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACTTTGAGTACGGTCAGTACAGAGTTATTCTGTTGGTTGAAGTATTccc                           <  2:403860/148‑1 (MQ=255)
                                                                                                                                aaTCCAGCTGAAAGCAAAACGTAACAGCAGCACGGTGGTTGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACTGTCATTAAAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTcac                   <  2:36532/149‑1 (MQ=255)
                                                                                                                                aaTCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACa                                                                                  <  7:230629/86‑1 (MQ=255)
                                                                                                                                aaTCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACa                                                                                  >  8:230629/1‑86 (MQ=255)
                                                                                                                                aaTCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTca                    >  7:29058/1‑148 (MQ=255)
                                                                                                                                    cAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGc               >  5:324201/1‑149 (MQ=255)
                                                                                                                                    cAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGc               >  1:92829/1‑149 (MQ=255)
                                                                                                                                     aGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGcc              <  3:76680/149‑1 (MQ=255)
                                                                                                                                       cTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGGCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCgg            <  5:287203/149‑1 (MQ=255)
                                                                                                                                          aaaGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGAc         >  2:19276/1‑149 (MQ=255)
                                                                                                                                          aaaGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAAACCGGATGCAGCCGGGCGTTACAGCATGGACGTAGAGTACGGTCAGTACAGCGTTAATCTGTGGGTGGAAGGATTCCCGCCGTCACATGCCgggcc         >  2:330687/1‑147 (MQ=255)
                                                                                                                                           aaGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGAcc        <  6:289743/149‑1 (MQ=255)
                                                                                                                                           aaGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGAc         >  6:380930/1‑148 (MQ=255)
                                                                                                                                              cAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCAt      >  8:262502/1‑148 (MQ=255)
                                                                                                                                               aaaaCGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTAc                                                                 <  4:126322/88‑1 (MQ=255)
                                                                                                                                               aaaaCGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTAc                                                                 >  3:126322/1‑88 (MQ=255)
                                                                                                                                               aaaaCGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTAAGTACAGAGTTATTCTGTTGGTGGAAGGATTCCCGCCGTACCATGCCGGGACCATTa    >  6:249657/1‑149 (MQ=255)
                                                                                                                                                aaaCGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTAc   >  8:406268/1‑149 (MQ=255)
                                                                                                                                                 aaCGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATACGGATGACGCCGGGCGTGACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTAGTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTAcc  >  5:5955/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACC  >  NC_000913/1428983‑1429276

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: