Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,423,372 A→T 10.3% intergenic (+178/+51) yhdZ → / ← rrfF putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,3720AT10.3% 1410.5 / 136.9 494intergenic (+178/+51)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base A (203/233);  new base T (40/11);  total (245/249)
Fisher's exact test for biased strand distribution p-value = 1.30e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

Only 100 of 536 total aligned reads displayed.

CCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCA  >  NC_000913/3423223‑3423521
                                                                                                                                                     |                                                                                                                                                     
 cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCGCCCGGCGGGTTTGTCCTACTCGGGGGGGTGTTCACCGACAaaaaac                                                                                                                                                      >  1:2993356/1‑150 (MQ=17)
  aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca                                                                                                                                                     <  1:8169908/150‑1 (MQ=255)
    ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGa                                                                                                                                                   <  2:1059570/150‑1 (MQ=255)
      tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATa                                                                                                                                                 >  2:7960836/1‑150 (MQ=255)
      tttCGGATGGCCTTTCGCTTGATTTCATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCACCCGGCGGATTTGTCCGACTCGGGAGAGTGTTACCCGACCACCACCAGATa                                                                                                                                                 >  2:6691867/1‑150 (MQ=14)
       ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaa                                                                                                                                                <  2:6205792/150‑1 (MQ=255)
         gcggtatgataaaacagaatttgcctggcggccgtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagACATCAAATaaaa                                                                                                                                              >  1:7910373‑M1/138‑150 (MQ=255)
            aTGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTAAAAAACCCACCCGGCGGATTTGTCCGACTCGGGCGAGTGTTCACCGACAAACAACAGAaaaaaaaa                                                                                                                                           >  1:9412163/1‑150 (MQ=14)
              ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACaaaa                                                                                                                                         <  2:5324225/150‑1 (MQ=255)
                 cTTTCGCTTGATTTGCTGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACACACAACAGATACAACAAAAGGc                                                                                                                                      <  2:751376/150‑1 (MQ=255)
                 cTTTCGCTTGATTGGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGc                                                                                                                                      <  2:6778593/150‑1 (MQ=255)
                 ataaaacagaatttgcctggcggccgtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagACATCAAATAAAACAAAAGGc                                                                                                                                      <  1:826186‑M1/21‑1 (MQ=255)
                  tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGCCAAACCACAGATAAAACAAAAGGcc                                                                                                                                     >  1:6819583/1‑150 (MQ=255)
                       cagaatttgcctggcggccgtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagACATCAAATaaaa                                                                                                                                              >  1:123570‑M1/124‑136 (MQ=255)
                       cTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACGTTCAAATACCCACCCCGCGGCTTTTTGCTACTCGGGTGAGTGTTGACCCACCAAAAACAGATAAAAAAAAAAGCCCACTc                                                                                                                                >  1:3098714/1‑150 (MQ=14)
                       cTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCCAATCCCCGCCCCGCGGAGTTGTCCTACGCGGGGGAGTGTTAACCGACAAACAACACATAAAAAAAAAGGCCCCGTc                                                                                                                                >  1:8724770/1‑150 (MQ=14)
                        ttgatttgatGTCTGGCAGTTTATGGCGGGCGTCCAGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCt                                                                                                                               >  2:9850522/1‑150 (MQ=255)
                         tgatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCtt                                                                                                                              <  2:8487513/150‑1 (MQ=255)
                           atttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTcc                                                                                                                            <  2:8949651/150‑1 (MQ=39)
                            tatgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAATTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                                                                                                           >  2:291776/3‑150 (MQ=255)
                               gatGTCTGGCAGTTTAGGGCGGGCGGCCGGCCCGCCACCCGCCGGGCCGTTGCTTCGCAACTTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACt                                                                                                                        <  1:10600713/150‑1 (MQ=14)
                                   tCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACCAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGc                                                                                                                    >  2:2237800/1‑150 (MQ=38)
                                   ggcggccgtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagACATCAAATAAAACAAAAGGCCCAGTCTTCCGACTGAGc                                                                                                                    >  2:5871562‑M1/112‑150 (MQ=255)
                                    gcggcagtagcgcggtggtcccacctgaccccatgccgaactcataagtgaaacgccgtagcgccgatggtagtgtgtggtctccccatgcgagagtagggaactgccagACATCAAATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                   >  2:3227434‑M1/111‑150 (MQ=255)
                                       gcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagGCATCAAATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCttt                                                                                                                >  2:9698039‑M1/108‑150 (MQ=255)
                                       gcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagGCATCAAATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCttt                                                                                                                >  2:6274262‑M1/108‑150 (MQ=255)
                                       gcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagGCATCAAATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCttt                                                                                                                >  2:2767909‑M1/108‑150 (MQ=255)
                                        cagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagGCATCAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCctttc                                                                                                               >  1:9469911‑M1/107‑149 (MQ=255)
                                         aGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg                                                                                                              <  2:11131228/150‑1 (MQ=38)
                                           tttATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTAGCCACCTTCAAAACCGCGCCCCGCGGATTTGTCCTACGCGGGAGGGTGTTAACCCACCAACAACAGATAAAAAAAAATGCCCAGGCGTCCGGCGGAGGCGTTTggt                                                                                                            >  1:7718787/1‑148 (MQ=2)
                                           tttATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACACAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                                                                                                            <  1:11011750/150‑1 (MQ=255)
                                            ttATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGttt                                                                                                           <  1:164451/150‑1 (MQ=39)
                                              cgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagGCATCAAATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTa                                                                                                         >  1:4273411‑M1/101‑150 (MQ=255)
                                                   gggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGa                                                                                                    <  2:3396163/150‑1 (MQ=39)
                                                   gggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGa                                                                                                    <  1:5464022/150‑1 (MQ=39)
                                                    ggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                                                                                                   <  1:8073537/150‑1 (MQ=39)
                                                         ccTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg                                                                                              <  1:398853/150‑1 (MQ=255)
                                                         ccTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg                                                                                              >  2:492345/1‑150 (MQ=255)
                                                          cTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                                                                                             >  2:4414995/1‑150 (MQ=255)
                                                          cTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTGTGTTTTATTTGATGTCTgg                                                                                             >  2:8063030/1‑150 (MQ=255)
                                                           tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc                                                                                            <  1:5826419/150‑1 (MQ=255)
                                                            gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACACACAACAGATAAAACAAAAGGCCCAGTCTTCCGTCTGAGCCTTTTGTTTTATTTGATGTCTGGCa                                                                                           <  1:10392469/150‑1 (MQ=255)
                                                            gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa                                                                                           <  1:1459424/150‑1 (MQ=255)
                                                                 ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc                                                                                      <  2:9049750/150‑1 (MQ=255)
                                                                    cccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTa                                                                                   <  2:1599916/150‑1 (MQ=255)
                                                                    cccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTa                                                                                   <  1:7519524/150‑1 (MQ=255)
                                                                        ccGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActct                                                                               <  1:1112229/150‑1 (MQ=255)
                                                                            gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCCGTCTTCCGACTGAGCCTTATGTTTTATTTGATGGCTGGCAGTTCCCTACTCTCGCa                                                                           >  2:10811200/1‑150 (MQ=255)
                                                                            gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCGACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCGTCCGACTGAGCCTTTGGGTTTAGTTGATGGCTGGCAGTTCCCTACTCTCGCa                                                                           >  2:3853031/1‑150 (MQ=11)
                                                                               gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCAggg                                                                        >  2:45815/1‑150 (MQ=255)
                                                                               gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATgg                                                                        <  2:3590828/150‑1 (MQ=255)
                                                                               gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATgg                                                                        <  1:4735544/150‑1 (MQ=255)
                                                                               gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATgg                                                                        <  1:2504096/150‑1 (MQ=255)
                                                                               gTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATgg                                                                        <  1:8901252/150‑1 (MQ=255)
                                                                                ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATggg                                                                       <  2:7843680/150‑1 (MQ=255)
                                                                                ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATggg                                                                       >  2:1937405/1‑150 (MQ=255)
                                                                                  gCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGga                                                                     <  1:6104058/150‑1 (MQ=255)
                                                                                   cTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGgag                                                                    >  1:6377088/1‑150 (MQ=255)
                                                                                     tCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAc                                                                  <  2:3696358/150‑1 (MQ=255)
                                                                                     tCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAc                                                                  >  2:2370157/1‑150 (MQ=255)
                                                                                      cGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAcc                                                                 >  1:6535063/1‑150 (MQ=255)
                                                                                        cAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGCGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTGCCGACTGAGCCTTTTGTTGTATTTGATGTCTGGCAGTACCCTACGCTCGCCTGGGGAGAcccc                                                               >  2:6516794/1‑150 (MQ=14)
                                                                                        cAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAcccc                                                               >  2:7752853/1‑150 (MQ=255)
                                                                                        cAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAcccc                                                               >  2:4247927/1‑150 (MQ=255)
                                                                                         aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACACAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCca                                                              <  2:6068003/150‑1 (MQ=255)
                                                                                         aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACACGGGAGAGTGTTAACCGACAAACAACAGATAAAACAAAAGGCCCAGTCATCCGACGGAGCCTTTTGTTTTATTTGATGTATGGCAGTTCCCAACAATCGCATGGGGAGACCCca                                                              >  1:2813361/1‑150 (MQ=14)
                                                                                           cGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcaca                                                            <  1:10944549/150‑1 (MQ=255)
                                                                                             ttCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACt                                                          >  1:10005595/1‑150 (MQ=255)
                                                                                               cAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTAc                                                        <  1:2065182/150‑1 (MQ=255)
                                                                                                aaaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGa                                                                                                    <  2:2674157/105‑1 (MQ=255)
                                                                                                aaaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTAcc                                                       <  2:2993356/150‑1 (MQ=255)
                                                                                                   tcctatCCCGTCGGATTTGTCCGCATCGGGAGAGTTTTCACCGACAAACAACATATAAAAAACAAGGCCCAGTCTTCCGACTGATCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATc                                                    <  2:7175190/145‑1 (MQ=14)
                                                                                                   tCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATc                                                    <  1:535780/150‑1 (MQ=255)
                                                                                                   tCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATc                                                    <  1:2652137/150‑1 (MQ=255)
                                                                                                      gCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGgc                                                 >  2:1406991/1‑150 (MQ=255)
                                                                                                       gtCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAACACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGgcg                                                <  2:10921685/149‑1 (MQ=255)
                                                                                                           cggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAc                                            <  2:8813498/150‑1 (MQ=255)
                                                                                                             gcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCAGCGGCGCTACgg                                          >  1:7318932/1‑150 (MQ=255)
                                                                                                               ggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCTACTGAGCGTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACACCACACTACCATCGGCGCTACGGCg                                        >  1:3705707/1‑150 (MQ=255)
                                                                                                               ggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCg                                        <  1:8478964/150‑1 (MQ=255)
                                                                                                               ggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCg                                        <  1:1939354/150‑1 (MQ=255)
                                                                                                               ggATTTGTCCGACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCg                                        <  2:6465934/150‑1 (MQ=255)
                                                                                                                    tGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCa                                   <  2:8259448/150‑1 (MQ=255)
                                                                                                                      tcctccTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGTCGTTTCACt                                 <  1:10285333/150‑1 (MQ=255)
                                                                                                                        cTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTc                               >  2:1026472/1‑150 (MQ=255)
                                                                                                                         tACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                              <  1:9136913/150‑1 (MQ=255)
                                                                                                                         tACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                              <  1:1715541/150‑1 (MQ=255)
                                                                                                                              gggAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGtt                         <  2:6988781/150‑1 (MQ=255)
                                                                                                                              gggAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGtt                         >  2:5656954/1‑150 (MQ=255)
                                                                                                                                 cgagTGTTCACCGACAACCAACAGATAAAACAAAAGGCCCAGTCGTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgg                      <  1:8817878/149‑1 (MQ=255)
                                                                                                                                 agagTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgt                      <  2:10421186/150‑2 (MQ=255)
                                                                                                                                 agagTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgg                      >  1:9972469/1‑150 (MQ=255)
                                                                                                                                      gtTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg                 >  2:5747676/1‑150 (MQ=39)
                                                                                                                                       ttCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTACCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATggg                <  1:964035/150‑1 (MQ=255)
                                                                                                                                          aCCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTc             >  1:5170314/1‑150 (MQ=39)
                                                                                                                                          aCCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTc             <  1:4375829/150‑1 (MQ=39)
                                                                                                                                          aCCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTc             <  2:9365755/150‑1 (MQ=39)
                                                                                                                                          aCCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTc             <  1:3476306/150‑1 (MQ=39)
                                                                                                                                                aaacaaCAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTgg       <  1:4950540/150‑1 (MQ=39)
                                                                                                                                                  acaacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGa     >  1:6700540/1‑150 (MQ=39)
                                                                                                                                                     |                                                                                                                                                     
CCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCA  >  NC_000913/3423223‑3423521

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