Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 464,356 A→T 5.9% intergenic (+48/‑46) 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,3560AT 5.9% 353.3 / 10.5 101intergenic (+48/‑46)ybaV/fadMputative competence‑suppressing periplasmic helix‑hairpin‑helix DNA‑binding protein/long‑chain acyl‑CoA thioesterase III
Reads supporting (aligned to +/- strand):  ref base A (49/46);  new base T (2/4);  total (51/50)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 106 total aligned reads displayed.

GCGAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATGGG  >  NC_000913/464208‑464503
                                                                                                                                                     |                                                                                                                                                   
gcgAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCa                                                                                                                                                     >  2:657050/1‑149 (MQ=255)
gcgAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCa                                                                                                                                                     >  1:405891/1‑149 (MQ=255)
 cgAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCa                                                                                                                                                     <  3:904836/148‑1 (MQ=255)
     gagTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTAt                                                                                                                                                <  2:43863/149‑1 (MQ=255)
     gagTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTAt                                                                                                                                                <  2:102407/149‑1 (MQ=255)
        tACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGa                                                                                                                                              <  3:1047034/148‑1 (MQ=255)
           ggTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATAACTGGTATGACctc                                                                                                                                          <  2:936659/149‑1 (MQ=255)
             tCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTg                                                                                                                                        >  2:360934/1‑149 (MQ=255)
               cGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTa                                                                                                                                      <  1:183656/149‑1 (MQ=255)
               cGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTa                                                                                                                                      <  1:629668/149‑1 (MQ=255)
               cGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAACTTTTTGCCAGACTGAAGAGGT‑CATAACTGGTATGACctcttca                                                                                                                                      <  1:540161/149‑4 (MQ=255)
                 tttAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACt                                                                                                                                    >  2:778257/1‑149 (MQ=255)
                 tttAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACt                                                                                                                                    >  1:426632/1‑149 (MQ=255)
                   tAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATg                                                                                                                                               >  4:883191/1‑136 (MQ=255)
                   tAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATg                                                                                                                                               <  3:883191/136‑1 (MQ=255)
                    aaaaCTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTtat                                                                                                                                 <  2:547610/149‑1 (MQ=255)
                      aaCTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATaa                                                                                                                               <  1:173390/149‑1 (MQ=255)
                              aggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATaacaacaac                                                                                                                        >  4:403046/1‑148 (MQ=255)
                              aggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATaacaacaac                                                                                                                        >  3:395947/1‑148 (MQ=255)
                               ggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACg                                                                                                                       >  4:532223/1‑148 (MQ=255)
                                 aTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTaa                                                                                                                    <  2:715180/149‑1 (MQ=255)
                                  tCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAg                                                                                                                   >  2:671729/1‑149 (MQ=255)
                                  tCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAg                                                                                                                   >  1:916140/1‑149 (MQ=255)
                                   cTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAgg                                                                                                                  <  1:892306/149‑1 (MQ=255)
                                       aGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTAt                                                                                                              <  1:374551/149‑1 (MQ=255)
                                         cAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTAtt                                                                                                             >  3:822491/1‑148 (MQ=255)
                                            gTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgcgc                                                                                                         >  2:732799/1‑149 (MQ=255)
                                              gCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTa                                                                                                       <  1:666201/149‑1 (MQ=255)
                                               tcGGGGATGGGCAATAAGCTTGAGGAACGTAATCTGGCGGTATTATACCTGTAATTAATTTGCACAGTGGCAATTTTTTGCCAGACTGCAGAGTT‑CATACCAGTTATGTCCCCTGTACTTATGACTGCAGCGTAAGGTTATTGCGCTAt                                                                                                      <  2:876870/148‑1 (MQ=255)
                                               ccGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTa                                                                                                       >  3:378935/1‑148 (MQ=255)
                                                    gATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                                 <  1:403639/149‑1 (MQ=255)
                                                    gATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                                 <  1:360934/149‑1 (MQ=255)
                                                    gATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaa                                                                                                  <  3:175770/148‑1 (MQ=255)
                                                     aTGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAac                                                                                                <  2:916140/149‑1 (MQ=255)
                                                      tGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                                 <  2:563978/147‑1 (MQ=255)
                                                      tGGGCAATTCGATGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                                 >  1:563978/1‑147 (MQ=255)
                                                         gCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGtt                                                                                                                                                  <  2:653746/95‑1 (MQ=255)
                                                          cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACaa                                                                                            <  4:287919/148‑1 (MQ=255)
                                                          cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACaa                                                                                            <  3:51222/148‑1 (MQ=255)
                                                          cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACaa                                                                                            <  3:779294/148‑1 (MQ=255)
                                                          cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTCATTGCGCTATGCAAACACaaa                                                                                           <  1:333710/149‑1 (MQ=255)
                                                            aTTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAAt                                                                                          >  4:169891/1‑148 (MQ=255)
                                                             ttCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCa                                                                                        >  1:607335/1‑149 (MQ=255)
                                                               cGCTGGTGGAACGTAATCTGGCGGTATTAACTCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCCCTGTACTTATACCAACACCGTAAGGTTATTGCGCTATGCAACCACAAATCaaa                                                                                      <  2:497458/149‑1 (MQ=255)
                                                               cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCaa                                                                                       <  4:395947/148‑1 (MQ=255)
                                                               cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCaa                                                                                       <  4:621477/148‑1 (MQ=255)
                                                                      ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGt                                                                                <  3:45818/148‑1 (MQ=255)
                                                                      ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTg                                                                               <  1:745523/149‑1 (MQ=255)
                                                                      ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGG‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTg                                                                               >  1:939154/1‑149 (MQ=255)
                                                                       gAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTg                                                                               >  3:569676/1‑148 (MQ=255)
                                                                        aaCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTgg                                                                              <  3:244848/148‑1 (MQ=255)
                                                                         aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                                 >  1:1078107/1‑128 (MQ=255)
                                                                         aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                                 <  2:1078107/128‑1 (MQ=255)
                                                                         aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGa                                                                             >  3:744/1‑148 (MQ=255)
                                                                          cGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGat                                                                            <  3:879087/148‑1 (MQ=255)
                                                                            tAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTa                                                                                                       <  4:813250/119‑1 (MQ=255)
                                                                            tAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTa                                                                                                       >  3:813250/1‑119 (MQ=255)
                                                                                 tGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCAtctc                                                                    >  2:706563/1‑149 (MQ=255)
                                                                                    cGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGa                                                                  >  3:967095/1‑148 (MQ=255)
                                                                                      gTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTAt                                                                                                      >  3:108096/1‑110 (MQ=255)
                                                                                      gTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTAt                                                                                                      <  4:108096/110‑1 (MQ=255)
                                                                                         ttAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTa                                                            >  2:802388/1‑149 (MQ=255)
                                                                                           aaCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAcc                                                          <  2:939154/149‑1 (MQ=255)
                                                                                            aCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCa                                                         >  1:808645/1‑149 (MQ=255)
                                                                                             cccTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCa                                                         <  4:668542/148‑1 (MQ=255)
                                                                                              ccTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGc                                                       <  2:98500/149‑1 (MQ=255)
                                                                                              ccTGTAAATTATTTGGCTTATGGCAAATTTTTGCCAGAATGGAGGGGTCCATA‑CAGGTAAGGCCTCTTTTCCTAAAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAg                                                        <  3:118688/148‑1 (MQ=255)
                                                                                               cTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGc                                                       <  4:822491/148‑1 (MQ=255)
                                                                                                 gTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCAc                                                     >  3:308428/1‑148 (MQ=255)
                                                                                                     ttaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgcg                                                                                                          >  4:811915/1‑91 (MQ=255)
                                                                                                     ttaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgcg                                                                                                          <  3:811915/91‑1 (MQ=255)
                                                                                                     ttaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTca                                                 >  4:183526/1‑148 (MQ=255)
                                                                                                                 gTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGCCCTCTGTACTTATAACAACACCGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTCCCAGCACGTCAACAACGCCCGCt                                     <  3:897229/148‑1 (MQ=255)
                                                                                                                   ggCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTAc                                   >  4:181129/1‑148 (MQ=255)
                                                                                                                    gCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTAcc                                  >  3:540198/1‑148 (MQ=255)
                                                                                                                     cAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCt                                 >  3:630781/1‑148 (MQ=255)
                                                                                                                      aaTTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTg                               >  2:1062204/1‑149 (MQ=255)
                                                                                                                      aaTTTTTTGCCAGAATGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTg                               >  1:807756/1‑149 (MQ=255)
                                                                                                                          ttttGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAAtt                           >  1:699363/1‑149 (MQ=255)
                                                                                                                            ttGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAAttt                          <  3:403046/148‑1 (MQ=255)
                                                                                                                                    gtacagagGT‑CATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGt                                                               >  1:732031/5‑103 (MQ=38)
                                                                                                                                    gtacagagGT‑CATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGt                                                               <  2:732031/99‑1 (MQ=38)
                                                                                                                                     tGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGaagaag                >  1:660531/1‑149 (MQ=255)
                                                                                                                                      gAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGaagaag                >  4:1010793/1‑148 (MQ=255)
                                                                                                                                         gagGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGAGGATATAATCTCGTCGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCg            >  2:132733/1‑149 (MQ=255)
                                                                                                                                          agGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaa                    >  1:933569/1‑140 (MQ=255)
                                                                                                                                          agGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaa                    <  2:933569/140‑1 (MQ=255)
                                                                                                                                          agGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCg            >  4:866348/1‑148 (MQ=255)
                                                                                                                                             t‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTg         >  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:131169/149‑1 (MQ=255)
                                                                                                                                                   ccAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATgg   >  1:522193/1‑149 (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)
                                                                                                                                                     |                                                                                                                                                   
GCGAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGT‑CATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATGGG  >  NC_000913/464208‑464503

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: ‑