Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,866,487 G→C 6.1% intergenic (+22/+72) ygbN → / ← rpoS putative transporter/RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,866,4870GC 6.1% 712.0 / 31.3 198intergenic (+22/+72)ygbN/rpoSputative transporter/RNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base G (94/91);  new base C (8/4);  total (102/96)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 203 total aligned reads displayed.

GGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTG  >  NC_000913/2866338‑2866636
                                                                                                                                                     |                                                                                                                                                      
      ggATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCg                                                                                                                                                  <  2:1706626/150‑1 (MQ=255)
        aTTGTGACCAAATCTCGGGGGTGGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGa                                                                                                                                                <  1:1174631/150‑1 (MQ=255)
         ttGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAc                                                                                                                                               >  1:2054060/1‑150 (MQ=255)
          tgtgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAcc                                                                                                                                              >  1:928148/1‑150 (MQ=255)
          tgtgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCGGAcc                                                                                                                                              <  1:1680490/150‑1 (MQ=255)
            tgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑tt                                                                                                                                           >  1:910610/1‑150 (MQ=255)
             gACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGACATCGGACC‑ttt                                                                                                                                          <  2:1947279/150‑1 (MQ=255)
              aCCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑tttt                                                                                                                                         >  1:1057348/1‑150 (MQ=255)
               ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTa                                                                                                                                        <  1:2656086/150‑1 (MQ=255)
               ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGAATCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTa                                                                                                                                        >  1:1497598/1‑150 (MQ=255)
                cAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTAt                                                                                                                                       >  2:1791394/1‑150 (MQ=255)
                cAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCGGACC‑ttttta                                                                                                                                       <  1:2200099/150‑3 (MQ=255)
                cAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGAGCGGATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCAGCGGACC‑GTTTAt                                                                                                                                       >  2:804276/1‑150 (MQ=255)
                 aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTAtt                                                                                                                                      >  1:920952/1‑150 (MQ=255)
                 aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTAtt                                                                                                                                      <  1:1827215/150‑1 (MQ=255)
                 aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCGGACC‑tttttat                                                                                                                                      <  2:1150749/150‑4 (MQ=255)
                    tatCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATtgtg                                                                                                                                   <  1:2422435/150‑1 (MQ=255)
                         ggggTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TGTTATTGTGCACAg                                                                                                                              >  1:210628/1‑150 (MQ=255)
                         ggggTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TGTTATTGTGCACAg                                                                                                                              >  1:1746851/1‑150 (MQ=255)
                            ggtGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGaaa                                                                                                                           >  1:2249681/3‑150 (MQ=255)
                            gTTGTCGTTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGaaa                                                                                                                           <  1:1706453/150‑1 (MQ=255)
                             ttGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGaaaa                                                                                                                          <  1:676062/150‑1 (MQ=255)
                             gtGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGaaaa                                                                                                                          <  1:2591802/149‑1 (MQ=255)
                               gTCGGTTGAAGAGGTTATTAATACCTGGAGAGTTTTTACGACCTTTCTAGTTTTTACATGTTTCTTAATTTCCTGGTGCTTAGTGGAGGTTTTTTTAACTTAATAAAAAGGTCCGATGTTCTTCTGACC‑TTTTATTGTGCTCTGAAAAgt                                                                                                                        <  2:2633745/150‑2 (MQ=255)
                                tCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGc                                                                                                                       >  2:247560/1‑150 (MQ=255)
                                tCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGc                                                                                                                       <  2:458244/150‑1 (MQ=255)
                                  ggTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCa                                                                                                                     <  2:1824773/150‑1 (MQ=255)
                                             gTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                          >  1:390621/1‑150 (MQ=255)
                                               cTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAACAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTgaga                                                                                                        <  2:950610/150‑1 (MQ=255)
                                               cTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTgaga                                                                                                        <  1:455776/150‑1 (MQ=255)
                                                  aaaaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGAACAGAAAAGGCCAGCCTCGCTTGAGACTg                                                                                                     >  2:63557/1‑150 (MQ=255)
                                                          gACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCt                                                                                             >  2:1914394/1‑150 (MQ=255)
                                                          gACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCt                                                                                             >  2:2455363/1‑150 (MQ=255)
                                                             aGTGTTAACGACCATTCACGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGAc                                                                                          >  2:1418277/1‑150 (MQ=255)
                                                              gtgtTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACa                                                                                         <  2:899842/150‑1 (MQ=255)
                                                               tgtTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAg                                                                                        >  2:2092388/1‑150 (MQ=255)
                                                                 ttAACGACCATTCTCTGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGTCGGTAATTTGACCTTAATAAAAAGGTCCGATTGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGAt                                                                                      <  2:336276/150‑1 (MQ=255)
                                                                 ttAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGAt                                                                                      <  2:1726510/150‑1 (MQ=255)
                                                                   aaCGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGc                                                                                    >  2:494120/1‑150 (MQ=255)
                                                                     cGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGctt                                                                                  <  2:1400700/150‑1 (MQ=255)
                                                                       aCCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttac                                                                                >  1:1143142/1‑150 (MQ=255)
                                                                        ccATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGTCCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttact                                                                               <  2:1648186/150‑1 (MQ=255)
                                                                        ccATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttact                                                                               >  1:541009/1‑150 (MQ=255)
                                                                        ccATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttact                                                                               >  1:535833/1‑150 (MQ=255)
                                                                          aTTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttactta                                                                             >  2:1602349/1‑150 (MQ=255)
                                                                           ttCTCGGTTTTACCGGATTCTTAATTACCTGGGGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGTAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCttacttac                                                                            <  2:398947/150‑1 (MQ=255)
                                                                              tcGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTcg                                                                         >  2:981979/1‑150 (MQ=255)
                                                                                         ggATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAgcgc                                                              <  2:113666/150‑1 (MQ=255)
                                                                                          gATTCTTACTTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGAGGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCt                                                             <  2:2384319/150‑1 (MQ=255)
                                                                                            ttCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTc                                                           <  2:1984423/150‑1 (MQ=255)
                                                                                               ttAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGat                                                        >  2:801905/1‑150 (MQ=255)
                                                                                               ttAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGat                                                        >  1:222011/1‑150 (MQ=255)
                                                                                               ttAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGat                                                        >  2:291663/1‑150 (MQ=255)
                                                                                                 aaTTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGatat                                                      >  1:563006/1‑150 (MQ=255)
                                                                                                   ttACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTc                                                    <  2:1676720/150‑1 (MQ=255)
                                                                                                          gTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCt                                             >  2:616789/1‑150 (MQ=255)
                                                                                                            gCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGc                                           <  1:708899/150‑1 (MQ=255)
                                                                                                            gCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGc                                           <  1:2084899/150‑1 (MQ=255)
                                                                                                              gTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGt                                         <  1:1962817/150‑1 (MQ=255)
                                                                                                              gTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGt                                         <  2:1578422/150‑1 (MQ=255)
                                                                                                               tATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGtt                                        <  2:2054060/150‑1 (MQ=255)
                                                                                                               tATGGGCGGTAATTTGACCTGAATACACAGGTCCGATGGGCATCGGACC‑TTTTAGTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACACGCGGAACAGCGCTTCGATATTCAGCCCCTGCGtt                                        <  2:1869204/150‑1 (MQ=255)
                                                                                                                   ggCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                    <  2:2622791/150‑1 (MQ=255)
                                                                                                                       gTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGAt                                >  2:1043820/1‑150 (MQ=255)
                                                                                                                       gTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGAt                                >  2:1089825/1‑150 (MQ=255)
                                                                                                                        tAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGAtt                               >  1:2523227/1‑150 (MQ=255)
                                                                                                                           tttGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTcg                            <  1:1696599/150‑1 (MQ=255)
                                                                                                                               aCCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCa                        >  1:2723895/1‑150 (MQ=255)
                                                                                                                                 cTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaaa                      <  1:2716746/150‑1 (MQ=255)
                                                                                                                                 cTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaaa                      >  1:2596397/1‑150 (MQ=255)
                                                                                                                                 cTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaaa                      >  1:2131890/1‑150 (MQ=255)
                                                                                                                                 cTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaaa                      >  1:2389266/1‑150 (MQ=255)
                                                                                                                                 cTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaaa                      >  2:2310944/1‑150 (MQ=255)
                                                                                                                                  ttAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAAc                     >  2:2371300/1‑150 (MQ=255)
                                                                                                                                   taataaAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAACCg                    <  1:407837/150‑1 (MQ=255)
                                                                                                                                     aataaaaGGTCCGATGCCCATCGGACCTTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACgg                   >  2:1261451/4‑150 (MQ=255)
                                                                                                                                       aaaaaGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCTGGATTTCGCGCAAACGgcgc                <  1:1027681/150‑1 (MQ=255)
                                                                                                                                       aaaaaGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGgcgc                <  2:2565794/150‑1 (MQ=255)
                                                                                                                                        aaaaGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCa               <  1:2390933/150‑1 (MQ=255)
                                                                                                                                        aaaaGGTCCGATGCCCATCGGACC‑TTTTAGTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCa               >  1:2065444/1‑150 (MQ=255)
                                                                                                                                         aaaGGTCCGATGGGCCTCGGACC‑TTTTATTGTGCACAGAAACGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCGTACTTACTCGCGGAACAGCGCTTCGAGATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg              <  2:2021235/150‑1 (MQ=255)
                                                                                                                                          aaGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg             >  2:860902/1‑150 (MQ=255)
                                                                                                                                           aGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACATAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTAGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc            >  1:744574/1‑150 (MQ=255)
                                                                                                                                           aGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAACCGGCGCAGGc            <  2:1498322/150‑1 (MQ=255)
                                                                                                                                           aGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc            >  1:232083/1‑150 (MQ=255)
                                                                                                                                           aGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc            <  1:1093856/150‑1 (MQ=255)
                                                                                                                                           aGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc            <  2:1369606/150‑1 (MQ=255)
                                                                                                                                            ggaccGAAGGGCATCTGACC‑TATTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc           <  2:2736999/147‑1 (MQ=255)
                                                                                                                                            ggTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc           >  2:2065842/1‑150 (MQ=255)
                                                                                                                                            ggTCCGATGCCCATCGGACC‑TTTTATGGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc           >  2:1689845/1‑150 (MQ=255)
                                                                                                                                             gTCCGATGGGCCTCGGACC‑TTTTATTGTGCCCAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCACACGGCGCAGGCCt          <  1:1482742/150‑1 (MQ=255)
                                                                                                                                             gTCCGATGCCCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt          >  1:4257/1‑150 (MQ=255)
                                                                                                                                             gTCCGATGCCCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt          >  1:2536239/1‑150 (MQ=255)
                                                                                                                                                 gATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaa      >  1:688376/1‑150 (MQ=255)
                                                                                                                                                 gATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaa      >  2:1197649/1‑150 (MQ=255)
                                                                                                                                                  aTGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCTCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc     >  2:44788/1‑150 (MQ=255)
                                                                                                                                                  aTGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc     >  2:2237515/1‑150 (MQ=255)
                                                                                                                                                   tGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCTTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaagc    <  2:743405/150‑3 (MQ=255)
                                                                                                                                                    gggCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt   >  1:1620079/1‑150 (MQ=255)
                                                                                                                                                     ggCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg  >  1:1326091/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTG  >  NC_000913/2866338‑2866636

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