Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 464,358 T→G 6.5% intergenic (+50/‑44) ybaV → / → fadM putative competence‑suppressing periplasmic helix‑hairpin‑helix DNA‑binding protein/long‑chain acyl‑CoA thioesterase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913464,3580TG 6.5% 362.4 / 13.9 105intergenic (+50/‑44)ybaV/fadMputative competence‑suppressing periplasmic helix‑hairpin‑helix DNA‑binding protein/long‑chain acyl‑CoA thioesterase III
Reads supporting (aligned to +/- strand):  ref base T (50/48);  new base G (2/5);  total (52/53)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 108 total aligned reads displayed.

GAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATGGGTTG  >  NC_000913/464213‑464506
                                                                                                                                                  |                                                                                                                                                    
gagTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTAt                                                                                                                                                   <  2:43863/149‑1 (MQ=255)
gagTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTAt                                                                                                                                                   <  2:102407/149‑1 (MQ=255)
   tACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGa                                                                                                                                                 <  3:1047034/148‑1 (MQ=255)
      ggTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctc                                                                                                                                             <  2:936659/149‑1 (MQ=255)
        tCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTg                                                                                                                                           >  2:360934/1‑149 (MQ=255)
          cGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTa                                                                                                                                         <  1:183656/149‑1 (MQ=255)
          cGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTa                                                                                                                                         <  1:629668/149‑1 (MQ=255)
          cGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAACTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttca                                                                                                                                         <  1:540161/149‑4 (MQ=255)
            tttAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACt                                                                                                                                       >  2:778257/1‑149 (MQ=255)
            tttAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACt                                                                                                                                       >  1:426632/1‑149 (MQ=255)
              tAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATg                                                                                                                                                  >  4:883191/1‑136 (MQ=255)
              tAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATg                                                                                                                                                  <  3:883191/136‑1 (MQ=255)
                 aaCTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaa                                                                                                                                  <  1:173390/149‑1 (MQ=255)
                         aggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaacaac                                                                                                                           >  3:395947/1‑148 (MQ=255)
                         aggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaacaac                                                                                                                           >  4:403046/1‑148 (MQ=255)
                          ggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACg                                                                                                                          >  4:532223/1‑148 (MQ=255)
                            aTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTaa                                                                                                                       <  2:715180/149‑1 (MQ=255)
                             tCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAg                                                                                                                      >  2:671729/1‑149 (MQ=255)
                             tCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAg                                                                                                                      >  1:916140/1‑149 (MQ=255)
                                  aGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTAt                                                                                                                 <  1:374551/149‑1 (MQ=255)
                                    cAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTAtt                                                                                                                >  3:822491/1‑148 (MQ=255)
                                       gTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgcgc                                                                                                            >  2:732799/1‑149 (MQ=255)
                                         gCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTa                                                                                                          <  1:666201/149‑1 (MQ=255)
                                          tcGGGGATGGGCAATAAGCTTGAGGAACGTAATCTGGCGGTATTATACCTGTAATTAATTTGCACAGTGGCAATTTTTTGCCAGACTGCAGAGTTCATACCAGTTATGTCCCCTGTACTTATGACTGCAGCGTAAGGTTATTGCGCTAt                                                                                                         <  2:876870/148‑1 (MQ=255)
                                               gATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                                    <  1:360934/149‑1 (MQ=255)
                                               gATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                                    <  1:403639/149‑1 (MQ=255)
                                               gATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaa                                                                                                     <  3:175770/148‑1 (MQ=255)
                                                aTGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAac                                                                                                   <  2:916140/149‑1 (MQ=255)
                                                 tGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                                    <  2:563978/147‑1 (MQ=255)
                                                    gCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGtt                                                                                                                                                     >  1:653746/1‑95 (MQ=255)
                                                    gCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGtt                                                                                                                                                     <  2:653746/95‑1 (MQ=255)
                                                    gCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAGGTTTATTGCGCTATCAAACACaa                                                                                               >  4:118688/1‑148 (MQ=255)
                                                     cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACaa                                                                                               <  3:51222/148‑1 (MQ=255)
                                                     cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACaa                                                                                               <  3:779294/148‑1 (MQ=255)
                                                     cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACaa                                                                                               <  4:287919/148‑1 (MQ=255)
                                                     cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTCATTGCGCTATGCAAACACaaa                                                                                              <  1:333710/149‑1 (MQ=255)
                                                       aTTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAAt                                                                                             >  4:169891/1‑148 (MQ=255)
                                                        ttCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCa                                                                                           >  1:607335/1‑149 (MQ=255)
                                                          cGCTGGTGGAACGTAATCTGGCGGTATTAACTCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCCCTGTACTTATACCAACACCGTAAGGTTATTGCGCTATGCAACCACAAATCaaa                                                                                         <  2:497458/149‑1 (MQ=255)
                                                          cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCaa                                                                                          <  4:621477/148‑1 (MQ=255)
                                                                 ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGt                                                                                   <  3:45818/148‑1 (MQ=255)
                                                                 ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTg                                                                                  <  1:745523/149‑1 (MQ=255)
                                                                 ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGGCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTg                                                                                  >  1:939154/1‑149 (MQ=255)
                                                                  gAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTg                                                                                  >  3:569676/1‑148 (MQ=255)
                                                                   aaCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTgg                                                                                 <  3:244848/148‑1 (MQ=255)
                                                                    aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                                    >  1:1078107/1‑128 (MQ=255)
                                                                    aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                                    <  2:1078107/128‑1 (MQ=255)
                                                                    aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGa                                                                                >  3:744/1‑148 (MQ=255)
                                                                     cGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGat                                                                               <  3:879087/148‑1 (MQ=255)
                                                                       tAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTa                                                                                                          >  3:813250/1‑119 (MQ=255)
                                                                       tAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTa                                                                                                          <  4:813250/119‑1 (MQ=255)
                                                                            tGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCAtctc                                                                       >  2:706563/1‑149 (MQ=255)
                                                                               cGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGa                                                                     >  3:967095/1‑148 (MQ=255)
                                                                                 gTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTAt                                                                                                         <  4:108096/110‑1 (MQ=255)
                                                                                 gTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTAt                                                                                                         >  3:108096/1‑110 (MQ=255)
                                                                                    ttAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTa                                                               >  2:802388/1‑149 (MQ=255)
                                                                                      aaCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAcc                                                             <  2:939154/149‑1 (MQ=255)
                                                                                       aCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCa                                                            >  1:808645/1‑149 (MQ=255)
                                                                                        cccTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCa                                                            <  4:668542/148‑1 (MQ=255)
                                                                                         ccTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGc                                                          <  2:98500/149‑1 (MQ=255)
                                                                                         ccTGTAAATTATTTGGCTTATGGCAAATTTTTGCCAGAATGGAGGGGTCCATACAGGTAAGGCCTCTTTTCCTAAAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAg                                                           <  3:118688/148‑1 (MQ=255)
                                                                                          cTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGc                                                          <  4:822491/148‑1 (MQ=255)
                                                                                           tGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCAc                                                        <  2:405891/149‑1 (MQ=255)
                                                                                            gTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCAc                                                        >  3:308428/1‑148 (MQ=255)
                                                                                                ttaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgcg                                                                                                             >  4:811915/1‑91 (MQ=255)
                                                                                                ttaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTca                                                    >  4:183526/1‑148 (MQ=255)
                                                                                                            gTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGCCCTCTGTACTTATAACAACACCGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTCCCAGCACGTCAACAACGCCCGCt                                        <  3:897229/148‑1 (MQ=255)
                                                                                                              ggCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTAc                                      >  4:181129/1‑148 (MQ=255)
                                                                                                               gCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTAcc                                     >  3:540198/1‑148 (MQ=255)
                                                                                                                cAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCt                                    >  3:630781/1‑148 (MQ=255)
                                                                                                                 aaTTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTg                                  >  2:1062204/1‑149 (MQ=255)
                                                                                                                 aaTTTTTTGCCAGAATGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTg                                  >  1:807756/1‑149 (MQ=255)
                                                                                                                     ttttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAAtt                              >  1:699363/1‑149 (MQ=255)
                                                                                                                       ttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAAttt                             <  3:403046/148‑1 (MQ=255)
                                                                                                                               gtacagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGt                                                                  >  1:732031/5‑103 (MQ=38)
                                                                                                                               gtacagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGt                                                                  <  2:732031/99‑1 (MQ=38)
                                                                                                                                tGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGaagaag                   >  1:660531/1‑149 (MQ=255)
                                                                                                                                 gAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGaagaag                   >  4:1010793/1‑148 (MQ=255)
                                                                                                                                    gagGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGAGGATATAATCTCGTCGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCg               >  2:132733/1‑149 (MQ=255)
                                                                                                                                     agGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaa                       >  1:933569/1‑140 (MQ=255)
                                                                                                                                     agGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaa                       <  2:933569/140‑1 (MQ=255)
                                                                                                                                     agGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCg               >  4:866348/1‑148 (MQ=255)
                                                                                                                                     agGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCg               >  3:48056/1‑148 (MQ=255)
                                                                                                                                        tCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTg            >  4:424353/1‑148 (MQ=255)
                                                                                                                                          cATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTgg           >  3:306339/1‑148 (MQ=255)
                                                                                                                                          cATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTgg           >  3:1030387/1‑148 (MQ=255)
                                                                                                                                           aTACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTggg          >  3:619602/1‑148 (MQ=255)
                                                                                                                                           aTACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGa         <  1:732799/149‑1 (MQ=255)
                                                                                                                                            tACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGAt        >  2:835816/1‑149 (MQ=255)
                                                                                                                                              ccAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATgg      >  1:522193/1‑149 (MQ=255)
                                                                                                                                              ccAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATgg      <  1:131169/149‑1 (MQ=255)
                                                                                                                                              ccAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATgg      >  1:61880/1‑149 (MQ=255)
                                                                                                                                              ccAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATgg      >  1:1038694/1‑149 (MQ=255)
                                                                                                                                                aGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGctgggctggg     >  4:793378/1‑148 (MQ=255)
                                                                                                                                                aGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATggg     <  3:183526/148‑1 (MQ=255)
                                                                                                                                                 gTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACaaa                                                                                              >  1:267071/1‑58 (MQ=255)
                                                                                                                                                 gTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACaaa                                                                                              <  2:267071/58‑1 (MQ=255)
                                                                                                                                                 gTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATGGGt    <  4:744/148‑1 (MQ=255)
                                                                                                                                                 gTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATGGGt    <  3:424353/148‑1 (MQ=255)
                                                                                                                                                  ttATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATGGGTTg  >  2:487783/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATGGGTTG  >  NC_000913/464213‑464506

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: