Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 518,730 A→G R140R (CGT→CGC)  ybbO ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913518,7300AG99.3% 575.0 / ‑6.1 152R140R (CGT→CGC) ybbOputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (81/70);  total (81/71)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 153 total aligned reads displayed.

GGGTTCGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAG  >  NC_000913/518583‑518877
                                                                                                                                                   |                                                                                                                                                   
 ggTTCGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGcc                                                                                                                                                 >  4:551247/1‑151 (MQ=255)
 ggTTCGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGcc                                                                                                                                                 >  2:520055/1‑151 (MQ=255)
    tCGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCt                                                                                                                                               <  2:1711161/150‑1 (MQ=255)
     cGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCtt                                                                                                                                              >  2:316183/1‑150 (MQ=255)
      gATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTc                                                                                                                                             >  1:677196/1‑150 (MQ=255)
      gATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCa                                                                                                                                            <  3:7459/151‑1 (MQ=255)
      gATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCCCCATGCGCAGTGCATCTGACCCCGCCTCCAGCGCATATTTACTGGCCGCGGAAGCGCCACGACCCGGCGCGGAGATTAATCCCATCACCGCTGATGTCATCACAATGCGCTCTTc                                                                                                                                             >  1:1528904/1‑150 (MQ=255)
       aTCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCAc                                                                                                                                           >  1:1025710/1‑151 (MQ=255)
          aGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGt                                                                                                                                        >  1:548847/1‑151 (MQ=255)
            gCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGc                                                                                                                                      <  4:44343/151‑1 (MQ=255)
                 cTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTa                                                                                                                                  >  2:122646/1‑150 (MQ=255)
                  tttAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAAc                                                                                                                                >  3:482880/1‑151 (MQ=255)
                   ttAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAAc                                                                                                                                <  1:499173/150‑1 (MQ=255)
                    tAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACAt                                                                                                                              >  2:150701/1‑151 (MQ=255)
                    tAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACAt                                                                                                                              >  2:81883/1‑151 (MQ=255)
                     aaTTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACAt                                                                                                                              <  3:124095/150‑1 (MQ=255)
                       ttCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATcgc                                                                                                                           >  2:489066/1‑151 (MQ=255)
                       ttCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATcgc                                                                                                                           >  1:1273178/1‑151 (MQ=255)
                          cGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCggg                                                                                                                        >  3:283203/1‑151 (MQ=255)
                             tgtgGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCggg                                                                                                                        >  2:823240/1‑148 (MQ=255)
                              gtgGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTa                                                                                                                      <  4:451321/149‑1 (MQ=255)
                              gtgGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTa                                                                                                                      <  2:1140128/149‑1 (MQ=255)
                                ggCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAg                                                                                                                  <  1:896615/151‑1 (MQ=255)
                                  cgcAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGgc                                                                                                                >  2:1135477/1‑151 (MQ=255)
                                  cgcAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGgc                                                                                                                <  2:508416/151‑1 (MQ=255)
                                      gCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCAt                                                                                                            <  2:1103455/151‑1 (MQ=255)
                                         ccATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATgg                                                                                                          >  2:747111/1‑150 (MQ=255)
                                         ccATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGt                                                                                                         <  2:304470/151‑1 (MQ=255)
                                          cATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTg                                                                                                        >  1:1131075/1‑151 (MQ=255)
                                           aTGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGa                                                                                                       <  2:95268/151‑1 (MQ=255)
                                            tGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTg                                                                                                        >  2:1405441/1‑149 (MQ=255)
                                               gcAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTg                                                                                                   <  2:37124/151‑1 (MQ=255)
                                               gcAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTg                                                                                                   <  2:83647/151‑1 (MQ=255)
                                                 aGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGgt                                                                                                 <  2:198209/151‑1 (MQ=255)
                                                  gTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGgtg                                                                                                <  2:1276568/151‑1 (MQ=255)
                                                    gCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGgtgt                                                                                               <  2:1105775/150‑1 (MQ=255)
                                                    gCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGtgtg                                                                                              >  1:1705255/1‑151 (MQ=255)
                                                      aTCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTgcg                                                                                            >  2:1540844/1‑151 (MQ=255)
                                                          gACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCg                                                                                         <  1:828477/150‑1 (MQ=255)
                                                             cACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGaaaa                                                                                     <  3:578742/151‑1 (MQ=255)
                                                             cACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGaaa                                                                                      >  2:1558532/1‑150 (MQ=255)
                                                               cGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGaaaaa                                                                                    <  2:1322372/150‑1 (MQ=255)
                                                               cGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAg                                                                                   <  2:1252765/151‑1 (MQ=255)
                                                                 ccTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGtt                                                                                 >  2:1014901/1‑151 (MQ=255)
                                                                 ccTCCAACGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCTATGATTTCATTTTAATACGCCCTGCGCCCCCCGGTAACATCTCGGGTAACCGTCCCTTGTTGCCCTAGCGTCCGCCGTAAAAGtc                                                                                 <  1:174012/151‑2 (MQ=255)
                                                                  cTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTg                                                                                >  2:604297/1‑151 (MQ=255)
                                                                  cTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTg                                                                                >  1:1181638/1‑151 (MQ=255)
                                                                    ccAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTg                                                                                <  1:321311/149‑1 (MQ=255)
                                                                    ccAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGc                                                                              >  3:192564/1‑151 (MQ=255)
                                                                     cAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCg                                                                             >  1:1641068/1‑151 (MQ=255)
                                                                           atatTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGaaaa                                                                        <  2:1390302/150‑1 (MQ=255)
                                                                             atTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAActg                                                                     <  1:1211950/151‑1 (MQ=255)
                                                                             atTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAct                                                                      <  2:160983/150‑1 (MQ=255)
                                                                              tttACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAActgc                                                                    >  1:465308/1‑151 (MQ=255)
                                                                                tACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAActgctg                                                                  >  1:1410668/1‑151 (MQ=255)
                                                                                tACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGAGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAActgctg                                                                  >  2:340755/1‑151 (MQ=255)
                                                                                   tGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTc                                                               >  4:419874/1‑151 (MQ=255)
                                                                                   tGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTc                                                               >  2:252800/1‑151 (MQ=255)
                                                                                    ggCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTcc                                                              <  1:1095202/151‑1 (MQ=255)
                                                                                    ggCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTc                                                               <  2:498764/150‑1 (MQ=255)
                                                                                    ggCAGCGTGAGCGGCACGACCCGGCGTGGAGATTGATCCCATCTCCGATGATGTCATCGCAATGCGCCCTTCACCGAGCGGTAACATCGCGGGTAACAGGCGCGTGGTGAGCTGGTGTGCGCCGAAAAAGTTAGCGGAAAACTGCAGTTc                                                               <  1:890683/150‑1 (MQ=255)
                                                                                      ccGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCAt                                                            >  1:1450671/1‑151 (MQ=255)
                                                                                          gTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTgc                                                        >  1:734473/1‑151 (MQ=255)
                                                                                              gcgcCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCAc                                                    >  1:677589/1‑151 (MQ=255)
                                                                                                 ccACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGc                                                 <  2:611/151‑1 (MQ=255)
                                                                                                 ccACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGc                                                 <  3:98697/151‑1 (MQ=255)
                                                                                                 ccACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGc                                                 >  2:1345820/1‑151 (MQ=255)
                                                                                                  cACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCt                                                <  2:1217162/151‑1 (MQ=255)
                                                                                                   aCGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTg                                               <  1:489843/151‑1 (MQ=255)
                                                                                                     gACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAt                                             >  1:918169/1‑151 (MQ=255)
                                                                                                      aCCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAtg                                            >  1:1339884/1‑151 (MQ=255)
                                                                                                      aCCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAtg                                            >  1:408693/1‑151 (MQ=255)
                                                                                                      aCCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAtg                                            <  1:913709/151‑1 (MQ=255)
                                                                                                       cccGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAtgg                                           <  1:30578/151‑1 (MQ=255)
                                                                                                       cccGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAtg                                            <  1:1037274/150‑1 (MQ=255)
                                                                                                        ccGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAtgg                                           >  2:749247/1‑150 (MQ=255)
                                                                                                         cGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAtggtg                                         >  2:1199222/1‑151 (MQ=255)
                                                                                                            cGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGaa                                      <  1:996659/151‑1 (MQ=255)
                                                                                                            cGTGGAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGaa                                      <  2:1605023/151‑1 (MQ=255)
                                                                                                               ggAGATTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAg                                    <  2:681246/150‑1 (MQ=255)
                                                                                                                gagaTTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAggg                                  >  4:4987/1‑151 (MQ=255)
                                                                                                                  gaTTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAggggg                                >  1:1546594/1‑151 (MQ=255)
                                                                                                                   aTTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAggggg                                >  2:684238/1‑150 (MQ=255)
                                                                                                                   aTTAATCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGc                               <  1:1273808/151‑1 (MQ=255)
                                                                                                                      aaTCCCATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCAt                            <  1:50300/151‑1 (MQ=255)
                                                                                                                         cccATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACa                         <  1:393439/151‑1 (MQ=255)
                                                                                                                           cATCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATg                       >  4:169787/1‑151 (MQ=255)
                                                                                                                            aTCACCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGc                      >  2:1487997/1‑151 (MQ=255)
                                                                                                                               aCCGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCg                    >  1:580427/1‑150 (MQ=255)
                                                                                                                                ccGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGa                   <  2:402988/150‑1 (MQ=255)
                                                                                                                                 cGATGATGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAAt                 >  3:459736/1‑151 (MQ=255)
                                                                                                                                  gatgatGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATc                <  2:1446381/151‑1 (MQ=255)
                                                                                                                                  gatgatGTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATc                >  1:351113/1‑151 (MQ=255)
                                                                                                                                   atgatgTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATcc               >  1:1331908/1‑151 (MQ=255)
                                                                                                                                   atgatgTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATcc               >  1:426108/1‑151 (MQ=255)
                                                                                                                                   atgatgTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATcc               >  2:321978/1‑151 (MQ=255)
                                                                                                                                    tgatgTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGTGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATcc               <  2:1278712/150‑1 (MQ=255)
                                                                                                                                        gTCATCACAATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCAt          <  2:448014/151‑1 (MQ=255)
                                                                                                                                             cacaATGCGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGtt      >  1:728481/1‑150 (MQ=255)
                                                                                                                                                aaTGCGCCCTTCACCGTGCGGAAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTaa    >  2:357128/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGGTTCGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAG  >  NC_000913/518583‑518877

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