Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,313,207 G→A 5.1% intergenic (‑29/+135) rbfA ← / ← infB 30s ribosome binding factor/translation initiation factor IF‑2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,313,2070GA 5.1% 530.4 / 12.1 158intergenic (‑29/+135)rbfA/infB30s ribosome binding factor/translation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  ref base G (69/81);  new base A (3/5);  total (72/86)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

Only 100 of 160 total aligned reads displayed.

ACCGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTA  >  NC_000913/3313059‑3313354
                                                                                                                                                    |                                                                                                                                                   
acCGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTg                                                                                                                                                     >  3:1063986/1‑149 (MQ=255)
acCGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTg                                                                                                                                                     >  2:720157/1‑149 (MQ=255)
        cGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGcccc                                                                                                                                             <  1:720157/149‑1 (MQ=255)
          gtggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAGGGGCTGTTAGCCCCttt                                                                                                                                          <  1:275191/149‑1 (MQ=255)
           tggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCttt                                                                                                                                          <  3:754555/149‑1 (MQ=255)
            ggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCtttt                                                                                                                                         <  1:158435/149‑1 (MQ=255)
            ggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTAACAGCCCCtttt                                                                                                                                         <  2:954686/149‑1 (MQ=255)
             gtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCttttt                                                                                                                                        <  3:918457/149‑1 (MQ=255)
              ccatcatGCCCAGACCTCGATCTTTAATTTCACGCTGAAGGATGAGAGCGCTCTATTTATGCATTTCCTGCGCTCAGGGCTGCCGGAGGCCGATTTCCCTCTCTATGATAAGTTCTCCTTCCATAAAAGGCGCCGTTAGACCCTTTTTa                                                                                                                                       <  1:989412/148‑1 (MQ=255)
               catcatGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTaa                                                                                                                                      <  2:464394/149‑1 (MQ=255)
                atcatGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTAACAGCCCCttt                                                                                                                                          <  4:1304192/144‑1 (MQ=255)
                atcatGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTAACAGCCCCttt                                                                                                                                          >  3:1304192/1‑144 (MQ=255)
                 tcatGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCGTTTTTaaaa                                                                                                                                    <  1:1164942/149‑1 (MQ=255)
                 tcatGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTaaaa                                                                                                                                    >  1:954959/1‑149 (MQ=255)
                 tcatGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTaaaa                                                                                                                                    <  3:717293/149‑1 (MQ=255)
                 tcatGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTaaaa                                                                                                                                    >  3:815849/1‑149 (MQ=255)
                  catGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaat                                                                                                                                   >  4:739460/1‑149 (MQ=255)
                  catGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTAACAGCCCCttttttgtca                                                                                                                                   <  4:434470/149‑6 (MQ=255)
                     gCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaattaa                                                                                                                                <  3:1330918/149‑1 (MQ=255)
                     gCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaattaa                                                                                                                                <  4:293081/149‑1 (MQ=255)
                      cccAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaattaat                                                                                                                               >  3:33189/1‑149 (MQ=255)
                       ccAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaattaatt                                                                                                                              >  3:8435/1‑149 (MQ=255)
                       ccAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaattaatt                                                                                                                              >  2:993039/1‑149 (MQ=255)
                         aGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTc                                                                                                                            >  2:1333102/1‑149 (MQ=255)
                          ggCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCa                                                                                                                           <  2:877936/149‑1 (MQ=255)
                                gATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGAATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGa                                                                                                                     >  2:463379/1‑149 (MQ=255)
                                 aTCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGaa                                                                                                                    <  4:633602/149‑1 (MQ=255)
                                 aTCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGaa                                                                                                                    >  2:1127745/1‑149 (MQ=255)
                                    tttAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAggg                                                                                                                 >  1:884298/1‑149 (MQ=255)
                                    tttAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAggg                                                                                                                 >  2:946168/1‑149 (MQ=255)
                                     ttAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGc                                                                                                                <  4:1272674/149‑1 (MQ=255)
                                       aaTTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTg                                                                                                              <  4:1206584/149‑1 (MQ=255)
                                        aTTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGt                                                                                                             <  3:403128/149‑1 (MQ=255)
                                        aTTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGt                                                                                                             >  1:1063751/1‑149 (MQ=255)
                                         tttCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGtt                                                                                                            >  2:880479/1‑149 (MQ=255)
                                              cGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGt                                                                                                       <  1:321885/149‑1 (MQ=255)
                                              cGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGt                                                                                                       <  4:1021505/149‑1 (MQ=255)
                                                 tGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGa                                                                                                    <  3:1177684/149‑1 (MQ=255)
                                                  gCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGAc                                                                                                   <  4:8435/149‑1 (MQ=255)
                                                      gATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGa                                                                                               <  4:615875/149‑1 (MQ=255)
                                                       aTGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGAt                                                                                              <  1:1299243/149‑1 (MQ=255)
                                                       aTGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGAt                                                                                              <  1:658106/149‑1 (MQ=255)
                                                        tGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATa                                                                                             >  2:580368/1‑149 (MQ=255)
                                                           gcgCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGa                                                                                          >  1:899320/3‑149 (MQ=255)
                                                            agCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGAc                                                                                         <  4:186879/149‑1 (MQ=255)
                                                                   tctTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGa                                                                                          <  1:1228727/141‑1 (MQ=255)
                                                                     tttttGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCg                                                                                >  2:910628/1‑149 (MQ=255)
                                                                     tttttGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCg                                                                                >  2:1072554/1‑149 (MQ=255)
                                                                         tGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTcac                                                                            <  1:506948/149‑1 (MQ=255)
                                                                            aTTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaattaat                                                                                                                               >  1:1121282/1‑95 (MQ=255)
                                                                            aTTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATc                                                                         >  2:1229566/1‑149 (MQ=255)
                                                                            aTTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATc                                                                         >  1:604047/1‑149 (MQ=255)
                                                                            aTTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATc                                                                         >  2:1357476/1‑149 (MQ=255)
                                                                             tttCCTGGGCTACGCGCTGCGGGCGACCAAATTCTTTCGGCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCa                                                                        >  4:290933/1‑149 (MQ=255)
                                                                             tttCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCa                                                                        >  1:313949/1‑149 (MQ=255)
                                                                              ttCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAg                                                                       <  2:324788/149‑1 (MQ=255)
                                                                              ttCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAg                                                                       <  4:1063986/149‑1 (MQ=255)
                                                                                  tGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGa                                                                                          >  2:815596/1‑126 (MQ=255)
                                                                                  tGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGa                                                                                          <  1:815596/126‑1 (MQ=255)
                                                                                     gcTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATcc                                                                <  2:1155130/149‑1 (MQ=255)
                                                                                        aCGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGtt                                                                                                      >  1:1286430/1‑108 (MQ=255)
                                                                                        aCGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGtt                                                                                                      <  2:1286430/108‑1 (MQ=255)
                                                                                         cgcgCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGc                                                            <  1:1333102/149‑1 (MQ=255)
                                                                                          gcgcTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGcc                                                           >  3:1144362/1‑149 (MQ=255)
                                                                                               gCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATg                                                      <  1:1072554/149‑1 (MQ=255)
                                                                                                   gCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGc                                                  >  2:410872/1‑149 (MQ=255)
                                                                                                   gCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGc                                                  <  2:604047/149‑1 (MQ=255)
                                                                                                      aCCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCt                                                                                                               <  2:502617/85‑1 (MQ=255)
                                                                                                       ccAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTaa                                              >  1:476845/1‑149 (MQ=255)
                                                                                                        cAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAGTTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGa                                                                                          >  1:1393307/1‑104 (MQ=255)
                                                                                                        cAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAGTTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGa                                                                                          <  2:1393307/104‑1 (MQ=255)
                                                                                                         aaaTTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAAc                                            >  1:346595/1‑149 (MQ=255)
                                                                                                         aaaTTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAAc                                            >  1:389562/1‑149 (MQ=255)
                                                                                                          aaTTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACg                                           >  4:1178381/1‑149 (MQ=255)
                                                                                                            ttCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGcc                                         <  3:937380/149‑1 (MQ=255)
                                                                                                            ttCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGc                                          >  2:615211/1‑148 (MQ=255)
                                                                                                                   gCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATccc                                  >  2:795377/1‑149 (MQ=255)
                                                                                                                        aataaATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATg                             >  4:1271466/1‑149 (MQ=255)
                                                                                                                        aataaATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATg                             <  4:33189/149‑1 (MQ=255)
                                                                                                                                 tcCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATaa                                                                                            <  1:550825/77‑1 (MQ=255)
                                                                                                                                 tcCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATaa                                                                                            >  2:550825/1‑77 (MQ=255)
                                                                                                                                   cTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACAt                                                                          <  2:1006507/93‑1 (MQ=255)
                                                                                                                                   cTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACAt                                                                          >  1:1006507/1‑93 (MQ=255)
                                                                                                                                    tGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTaaaaa                 <  2:1120543/149‑1 (MQ=255)
                                                                                                                                     gACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAAc                <  3:1090295/149‑1 (MQ=255)
                                                                                                                                       cAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCt              <  3:610563/149‑1 (MQ=255)
                                                                                                                                       cAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCt              <  3:152427/149‑1 (MQ=255)
                                                                                                                                        aaaaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATGAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCtt             >  4:1004496/1‑149 (MQ=255)
                                                                                                                                        aaaaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCtt             >  4:677340/1‑149 (MQ=255)
                                                                                                                                         aaaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTa            <  3:155947/149‑1 (MQ=255)
                                                                                                                                          aaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTaa           <  3:623484/149‑1 (MQ=255)
                                                                                                                                          aaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTaa           >  4:710895/1‑149 (MQ=255)
                                                                                                                                          aaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTaa           <  4:1324290/149‑1 (MQ=255)
                                                                                                                                           aaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAg          <  3:1131826/149‑1 (MQ=255)
                                                                                                                                           aaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAg          <  3:710895/149‑1 (MQ=255)
                                                                                                                                           aaaGGGGCTAACAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAg          >  4:1239445/1‑149 (MQ=255)
                                                                                                                                           aaaGGGGCTAACAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAg          >  1:1014980/1‑149 (MQ=255)
                                                                                                                                            aaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGc         >  1:420728/1‑149 (MQ=255)
                                                                                                                                               gggCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAAt      <  2:1440638/149‑1 (MQ=255)
                                                                                                                                                   tGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTa  >  3:544099/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACCGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTA  >  NC_000913/3313059‑3313354

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: