Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,744,148 A→G 6.9% intergenic (+41/+35) yfiB → / ← rplS OM lipoprotein putative positive effector of YfiN activity/50S ribosomal subunit protein L19

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,744,1480AG 6.9% 469.2 / 17.5 145intergenic (+41/+35)yfiB/rplSOM lipoprotein putative positive effector of YfiN activity/50S ribosomal subunit protein L19
Reads supporting (aligned to +/- strand):  ref base A (67/68);  new base G (5/5);  total (72/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 145 total aligned reads displayed.

CTCACCACACAGGGTTTAGGAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGTTTGACAGAAATGCTGTCAACTACCG  >  NC_000913/2744003‑2744295
                                                                                                                                                 |                                                                                                                                                   
  caccacACAGGGTTTAGGAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCct                                                                                                                                                >  3:181304/1‑149 (MQ=255)
   accacACAGGGTTTAGGAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCctt                                                                                                                                               >  4:744767/1‑149 (MQ=255)
   accacACAGGGTTTAGGAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCctt                                                                                                                                               >  1:531031/1‑149 (MQ=255)
    ccacACAGGGTTTAGGAAAAAAATATCCCATATCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTTATTACTACCCCTTAAAACCCCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCcttc                                                                                                                                              <  3:742868/149‑1 (MQ=255)
    ccacACAGGGTTTAGGAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCcttc                                                                                                                                              >  4:251728/1‑149 (MQ=255)
      acacaGGGTTTAGGAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCcttctt                                                                                                                                            <  1:469248/149‑1 (MQ=255)
        acaGGGTTTAGGAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTaa                                                                                                                                          >  1:578021/1‑149 (MQ=255)
            ggTTTAGGAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAgg                                                                                                                                      >  1:34144/1‑149 (MQ=255)
             gTTTAGGAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGa                                                                                                                                     >  4:112916/1‑149 (MQ=255)
             gTTTAGGAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGa                                                                                                                                     <  2:357283/149‑1 (MQ=255)
                  ggAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCg                                                                                                                                >  3:816289/1‑149 (MQ=255)
                   gAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCGGCCCTTCTTAACAGGATGt                                                                                                                                  <  4:86227/146‑1 (MQ=255)
                   gAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCGGCCCTTCTTAACAGGATGt                                                                                                                                  >  3:86227/1‑146 (MQ=255)
                   gAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGc                                                                                                                               >  3:25464/1‑149 (MQ=255)
                   gAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGc                                                                                                                               >  3:750905/1‑149 (MQ=255)
                   gAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGc                                                                                                                               >  3:558793/1‑149 (MQ=255)
                     aaaaaaTATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCtt                                                                                                                             >  2:624155/1‑149 (MQ=255)
                       aaaaTATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTaa                                                                                                                           <  1:820658/149‑1 (MQ=255)
                       aaaaTATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTaa                                                                                                                           <  4:708098/149‑1 (MQ=255)
                         aaTATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGc                                                                                                                         <  2:721853/149‑1 (MQ=255)
                          atatCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCg                                                                                                                        >  3:335113/1‑149 (MQ=255)
                            atCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAGCAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGaa                                                                                                                      >  4:710666/1‑149 (MQ=255)
                              cccATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAAt                                                                                                                    <  3:161927/149‑1 (MQ=255)
                               ccATAGCCAGTCACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATAGGCCAGCCCTTCTTGACAGGATGTCGCTTAAGCGAAATc                                                                                                                   >  3:581811/1‑149 (MQ=255)
                               ccATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATc                                                                                                                   >  3:232310/1‑149 (MQ=255)
                               ccATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATc                                                                                                                   >  1:864985/1‑149 (MQ=255)
                                cATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCt                                                                                                                  >  2:883180/1‑149 (MQ=255)
                                 aTAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCtt                                                                                                                 >  3:530560/1‑149 (MQ=255)
                                  tAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGGGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTa                                                                                                                >  3:600096/1‑149 (MQ=255)
                                  tAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCTCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTa                                                                                                                >  3:770984/1‑149 (MQ=255)
                                  tAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTa                                                                                                                <  4:176238/149‑1 (MQ=255)
                                  tAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTa                                                                                                                <  3:653878/149‑1 (MQ=255)
                                  tAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTa                                                                                                                <  1:352933/149‑1 (MQ=255)
                                  tAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAACACCGCCGCGTCGCAGGGGTGATTACCACTCCTTAAACCCCCAGCAAAGACAAAAAAAAGGGCCAGCTAATTGGCAAGCCCTTCTCATCAGGATTTCGCTTAAGATAAATCCTa                                                                                                                <  2:247466/149‑1 (MQ=255)
                                    gCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGt                                                                                                              <  3:698938/149‑1 (MQ=255)
                                        gTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAg                                                                                                          >  3:203457/1‑149 (MQ=255)
                                         tAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGa                                                                                                         <  4:558793/149‑1 (MQ=255)
                                             aaGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGct                                                                                                     >  3:414707/1‑149 (MQ=255)
                                             aaGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGct                                                                                                     <  2:728032/149‑1 (MQ=255)
                                             aaGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGct                                                                                                     <  3:869585/149‑1 (MQ=255)
                                                aCCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCtt                                                                                                  >  4:309691/1‑149 (MQ=255)
                                                aCCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCtt                                                                                                  >  1:456833/1‑149 (MQ=255)
                                                aCCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCtt                                                                                                  <  2:398716/149‑1 (MQ=255)
                                                   ccccAGGGCCGCGCCGAGCACCGCCGCGGCGCAGTGCTGCTTACTGCCCCTGAACACCCCCGCAAACAGATAAAAAAGGGCCAGCCAACTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAACCTTAGTTAAGCCGGTCTTTGa                                                                                               <  2:486590/148‑1 (MQ=255)
                                                     ccAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATa                                                                                             <  2:34144/149‑1 (MQ=255)
                                                      cAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATAc                                                                                            <  2:585702/149‑1 (MQ=255)
                                                      cAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATAc                                                                                            <  2:144862/149‑1 (MQ=255)
                                                       aGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACg                                                                                           >  2:532571/1‑149 (MQ=255)
                                                       aGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACg                                                                                           >  1:528464/1‑149 (MQ=255)
                                                        gggCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTAAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGa                                                                                          <  3:286859/149‑1 (MQ=255)
                                                           ccGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGagca                                                                                       >  4:186051/1‑149 (MQ=255)
                                                                       ccgccgCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAATGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTa                                                                           >  2:28915/1‑149 (MQ=255)
                                                                            gcgTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTc                                                                      <  4:395516/149‑1 (MQ=255)
                                                                                 gCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTAc                                                                                >  1:772538/1‑134 (MQ=255)
                                                                                  cAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTAAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAg                                                                <  4:428063/149‑1 (MQ=255)
                                                                                    gtggtgATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGgt                                                              >  1:193423/1‑149 (MQ=255)
                                                                                    gtggtgATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGgt                                                              <  4:25464/149‑1 (MQ=255)
                                                                                     tggtgATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGgta                                                             >  4:489711/1‑149 (MQ=255)
                                                                                     tggtgATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGgta                                                             >  1:726158/1‑149 (MQ=255)
                                                                                       gtgATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGa                                                                                                         >  4:664857/1‑103 (MQ=255)
                                                                                         gATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTAc                                                         <  3:646140/149‑1 (MQ=255)
                                                                                          aTTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACa                                                        <  1:208400/149‑1 (MQ=255)
                                                                                            tactacCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGt                                                      >  1:398849/1‑149 (MQ=255)
                                                                                               tacCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGgtag                                                            >  4:671525/1‑140 (MQ=255)
                                                                                               tacCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGgtag                                                            <  3:671525/140‑1 (MQ=255)
                                                                                               tacCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACa                                                        >  4:474699/1‑144 (MQ=255)
                                                                                               tacCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTa                                                   <  1:532571/149‑1 (MQ=255)
                                                                                               tacCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTa                                                   <  2:193423/149‑1 (MQ=255)
                                                                                                acCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAg                                                  >  1:60479/1‑149 (MQ=255)
                                                                                                   ccTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGaaa                                                                                                                     >  3:388933/1‑79 (MQ=255)
                                                                                                   ccTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGaaa                                                                                                                     <  4:388933/79‑1 (MQ=255)
                                                                                                     ttAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTa                                             <  1:883180/149‑1 (MQ=255)
                                                                                                       aaaaCACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACg                                           <  3:133779/149‑1 (MQ=255)
                                                                                                         aaCACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGaa                                         <  4:35143/149‑1 (MQ=255)
                                                                                                         aaCACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGaa                                         <  4:203457/149‑1 (MQ=255)
                                                                                                            acCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAg                                      <  3:346629/149‑1 (MQ=255)
                                                                                                            acCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAg                                      <  2:554083/149‑1 (MQ=255)
                                                                                                                gCAAACAGATAAAAAAGGGCCAGCCAATTGGACAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGGTACGAGCAGACTTACCAGTACGCTCACGCAGGTAGTACAGtt                                                     <  2:807141/130‑1 (MQ=255)
                                                                                                                  aaaCAGATAACAAAGGGCCAGCCAATTGGCCAGCTCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGcaccac                                >  2:534148/1‑149 (MQ=255)
                                                                                                                  aaaCAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACg                                           <  3:393988/138‑1 (MQ=255)
                                                                                                                  aaaCAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACg                                           >  4:393988/1‑138 (MQ=255)
                                                                                                                  aaaCAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGcaccac                                >  3:696254/1‑149 (MQ=255)
                                                                                                                   aaCAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCacg                               >  2:257246/1‑149 (MQ=255)
                                                                                                                   aaCAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCacg                               >  1:39596/1‑149 (MQ=255)
                                                                                                                   aaCAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCacg                               >  3:96109/1‑149 (MQ=255)
                                                                                                                       gATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGt                           <  3:489711/149‑1 (MQ=255)
                                                                                                                       gATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGt                           >  2:796438/1‑149 (MQ=255)
                                                                                                                        aTAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGtt                          <  2:795200/149‑1 (MQ=255)
                                                                                                                         tAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGttt                         <  1:665978/149‑1 (MQ=255)
                                                                                                                         tAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGttt                         <  1:821453/149‑1 (MQ=255)
                                                                                                                           aaaaaGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGTTTGa                       <  3:254426/149‑1 (MQ=255)
                                                                                                                           aaaaaGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGTTTGa                       <  2:60479/149‑1 (MQ=255)
                                                                                                                             gaaGGGCTGGCCAATTGGCTGGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGTTTGACa                     >  1:375716/2‑149 (MQ=255)
                                                                                                                             gaaGGGCTGGCCAATTGGCTGGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTAACAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGTTTGACa                     >  3:675330/2‑149 (MQ=255)
                                                                                                                                gggctggccaATTGGCTGGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGTTTGACAGaa                  >  2:479596/7‑149 (MQ=255)
                                                                                                                                  gccagccaATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGTTTGACAGAAAt                <  3:744767/149‑1 (MQ=255)
                                                                                                                                    cagcaaATTGGCCAGCCCTTCTTAACATGATGTCGCTTAAGCGAAATATTAGTTAAGACGCTATTTGATTCGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGACCAGCACCGCGACGTTTGACAGAAATGc              <  4:23832/149‑1 (MQ=255)
                                                                                                                                          aTTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGTTTGACAGAAATGCTGTCaa        <  4:600096/149‑1 (MQ=255)
                                                                                                                                           ttGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGTTTGACAGAAATGCTGTCAAc       >  3:448374/1‑149 (MQ=255)
                                                                                                                                                cAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGTTTGACAGAAATGCTGTCAACTACCg  <  1:772788/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTCACCACACAGGGTTTAGGAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGACGTTTGACAGAAATGCTGTCAACTACCG  >  NC_000913/2744003‑2744295

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