Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 222,958 Δ1 bp 100% coding (126/576 nt) gmhB → D,D‑heptose 1,7‑bisphosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913222,9580A.88.4% 233.8 / 27.9 69coding (126/576 nt)gmhBD,D‑heptose 1,7‑bisphosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base . (33/28);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

TTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCTAT  >  NC_000913/222813‑223102
                                                                                                                                                 |                                                                                                                                                
ttattaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaa                                                                                                                                               >  1:31767/1‑149 (MQ=255)
  attaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaa                                                                                                                                             >  2:902331/1‑149 (MQ=255)
    taAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaat                                                                                                                                           <  1:751112/149‑2 (MQ=255)
    taAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGGGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaat                                                                                                                                           <  1:469359/149‑2 (MQ=255)
      aGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatgg                                                                                                                                         >  2:753659/1‑146 (MQ=255)
      aGGAGCTATAAAAGGCGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatgg                                                                                                                                         >  1:873823/1‑146 (MQ=255)
       ggAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGATATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatggg                                                                                                                                        <  1:673245/149‑5 (MQ=255)
          gCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGGGTTTTTGACGCCATGCGCGAGCTAAAAAAAATGGGCt                                                                                                                                     >  2:185694/1‑149 (MQ=255)
           ctttaAAAGGTGGCGAAGAGCGTACCCGAAATTTTTCTTGACCGTGATGGCTCCATTAATGTCGATCACGGCTATGGCCATGAGTTCGACAACTTTGAATTTATCGCCGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGTCttt                                                                                                                                   <  1:316404/146‑1 (MQ=255)
           cTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCttt                                                                                                                                   >  1:300437/1‑149 (MQ=255)
            tataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGCTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCttt                                                                                                                                   <  1:295948/149‑1 (MQ=255)
              taAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgcg                                                                                                                                >  1:39176/1‑149 (MQ=255)
                aaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCt                                                                                                                              >  2:750907/1‑149 (MQ=255)
                aaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCt                                                                                                                              <  2:80192/149‑1 (MQ=255)
                 aaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCt                                                                                                                              >  2:165557/1‑149 (MQ=255)
                 aaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtg                                                                                                                             >  2:720019/1‑149 (MQ=255)
                  aGGTGGCGAAGAGCGTTCCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACTACTTTGAATTTATCTACGGTGTTATTGACGCCATGCGCGAGCT‑TAATAAATGGGCTTTGCGCtgg                                                                                                                            >  2:1016715/1‑149 (MQ=255)
                     gggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggtgg                                                                                                                         <  1:302382/148‑1 (MQ=255)
                      ggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTggtggt                                                                                                                        >  1:723479/1‑149 (MQ=255)
                      ggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTggtggt                                                                                                                        <  2:288537/149‑1 (MQ=255)
                       gCGAAGAGCGTACCCGCAATTTTTCTTGCCCGTGATCGCACCATTAATGTCGCTCACGGCTATGTCCATTAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGTGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                       <  2:910297/149‑1 (MQ=255)
                        cGAAGAGCGTAGCCGCAATTGTGCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGGCCATGAGATCGACAACTTTGAATTTATCGACCGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTGGCGCTGGTGggag                                                                                                                      <  2:216167/149‑1 (MQ=255)
                        cGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtag                                                                                                                      >  2:302567/1‑149 (MQ=255)
                         gAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtagt                                                                                                                     >  2:487549/1‑149 (MQ=255)
                          aaGACCGTACCCGCAATTTATCATGACCGTGAAGGCACCAATAATGTCGATCCCGGCTATGTCCATGAGATCGACAACTTAGAATTTAACGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtagta                                                                                                                    <  2:714876/149‑1 (MQ=255)
                                gTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                              <  2:300437/149‑1 (MQ=255)
                                    ccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGt                                                                                                          >  1:974344/1‑149 (MQ=255)
                                     cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACAAGTc                                                                                                         >  2:639488/1‑149 (MQ=255)
                                      gCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGAAACCAACCAGTCt                                                                                                        >  2:14056/1‑149 (MQ=255)
                                       cattttttCTTGACCGTTAGGGCACCATTAATGTCGATCACGGCTATGTCCTTGAGATCGACAACTTTGATTTTATCGACGGTGTTATTGACGCCAGGCGCGAGCT‑AAAAAAATTGGCTTTGCGCGGGTGGTAGTAACCAACCAGTCTg                                                                                                       <  1:581503/146‑1 (MQ=255)
                                            tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAtt                                                                                                  >  2:370670/1‑149 (MQ=255)
                                            tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGAAtt                                                                                                  >  1:289112/1‑149 (MQ=255)
                                            tttCTTGACCGTGATGGCACCATTAATGTCGATCAAGGCTATGTCCATGAGATCGACAACTTTGTATTTATCGACGGTGTTATTGACTCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAtt                                                                                                  <  2:919101/149‑1 (MQ=255)
                                                ttGACCGTGATGGCACCATTAATGTCGATCACCGCTATGTACATGAGATCGACAACATTGATTTAATCGACGGTGTAATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGCTGGTAGTAACCAACCTGTCTGGCATTGCTc                                                                                              >  2:152408/1‑149 (MQ=255)
                                                 tGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTACCCCACCAGTCTGGCATTGCTcg                                                                                             >  1:936063/1‑149 (MQ=255)
                                                        gATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaaa                                                                                      >  1:845147/1‑149 (MQ=255)
                                                             ccccATTCATCCCGATCCCGGCTATGGCCCTGAGATCGACAACTAGGAATATATCGCCGGTGGTCTTGACGCCATGCGCGAGCT‑AAAACACTGGGCTTTGCGCTGGGGGTAGTCACCAACCCGTCTGGCAGTGCTCGCGGTAACTTTAc                                                                                 <  2:12127/147‑1 (MQ=255)
                                                             cACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAc                                                                                 <  2:873823/149‑1 (MQ=255)
                                                             cACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTAGTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGACTGGCATTGCTCGCGGTAAATTTAc                                                                                 >  2:82396/1‑149 (MQ=255)
                                                              cccATTAATGTCGATCACGGCTATGTCCCTGATAGCGACAACTTTCTATTAATCGACGGTGTTATTGACGCCAAGCGCGAGCT‑AAAAACATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAcc                                                                                <  2:159270/148‑1 (MQ=255)
                                                              aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAcc                                                                                >  1:885757/1‑149 (MQ=255)
                                                                cATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGa                                                                              >  1:314830/1‑149 (MQ=255)
                                                                  ttAATGTCTATCTCGGCTATGTCCTGGAGATCGACAACTTTGTATTTTTCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTGGGGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAg                                                                            <  1:676341/149‑1 (MQ=255)
                                                                  ttAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAg                                                                            <  1:134100/149‑1 (MQ=255)
                                                                  ttAATGTCGATCAAGTCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTAGCGGTAATTTTACCGAAg                                                                            >  1:86707/1‑149 (MQ=255)
                                                                   tAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGc                                                                           <  1:975978/149‑1 (MQ=255)
                                                                   tAATGTCGATCAAGGCTATGTCCATGAGATCGACAACTTTGACATTCTCGACTGTGTTATAGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGc                                                                           <  1:979330/149‑1 (MQ=255)
                                                                       gTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAg                                                                       >  2:46805/1‑149 (MQ=255)
                                                                         cGATCACGGGTATGTCCATTAGTTCGACAACTTTGAATTTATCGTCGGTGTTATTGACGCGATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCCGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGtt                                                                     <  1:263372/149‑1 (MQ=255)
                                                                            tCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATTGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATATACAGAAGCACAGTTTGa                                                                  >  1:746913/1‑149 (MQ=255)
                                                                            tCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGa                                                                  >  1:332535/1‑149 (MQ=255)
                                                                              aCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                                <  1:880692/149‑1 (MQ=255)
                                                                                ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACTGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGGCTGGCATTGCTCGCGGTAAAGTTAGCGAAGCACAGTTTGAAACg                                                              >  2:448222/1‑149 (MQ=255)
                                                                                      gTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                        >  1:103848/1‑149 (MQ=255)
                                                                                       tCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCg                                                       >  2:1007599/1‑149 (MQ=255)
                                                                                         cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAg                                                     <  1:750907/149‑1 (MQ=255)
                                                                                              gaTCGCCAACTATGAATTCATCGACGGTGTTCATGCCGCCAGGCGCGAGC‑GAAAAAAATGGGCTTTGCGCTGGGGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGGAACTTTACCGACGCACAGTTTGACACGCGGACCGCGtggat                                                <  2:484707/149‑1 (MQ=255)
                                                                                              gaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGAAGTAACCAACCAGACTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTCATCCGCTGACCTAGtgga                                                 >  2:506492/1‑149 (MQ=255)
                                                                                              gaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTACCCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAACCGCTGACCGAGtggat                                                >  2:948904/1‑149 (MQ=255)
                                                                                                 cGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatgga                                             <  1:514927/149‑1 (MQ=255)
                                                                                                   aCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACt                                           <  2:437027/149‑1 (MQ=255)
                                                                                                        tttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCg                                      >  2:197062/1‑149 (MQ=255)
                                                                                                        tttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCg                                      >  1:457709/1‑149 (MQ=255)
                                                                                                         ttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGc                                     <  1:327822/149‑1 (MQ=255)
                                                                                                            aaTTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTgg                                  <  2:723479/149‑1 (MQ=255)
                                                                                                               ttATCGACGGTGTTATTGACTCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGTTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCgg                               >  1:369994/1‑149 (MQ=255)
                                                                                                                             aTTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt                 <  2:974344/149‑1 (MQ=255)
                                                                                                                             aTTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt                 >  2:699728/1‑149 (MQ=255)
                                                                                                                                 aCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTgg             >  1:515820/1‑149 (MQ=255)
                                                                                                                                   gCCATGCGCGAGCT‑AAAAAACTGGGCTTTGCGCTGGTGGGAGTAACCAACCAGTCTCGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATCGACTGGGCGCTGGCGGACCGAGATGTCGATCTGGAt           <  2:234371/149‑1 (MQ=255)
                                                                                                                                   gCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt           <  1:370670/149‑1 (MQ=255)
                                                                                                                                   gCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGAGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt           >  2:204124/1‑149 (MQ=255)
                                                                                                                                          gcgAGCTAAAAAAAATGGTCTTTGCTCTGGTTGTAGTAACCAACCAGTGTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATc     <  1:166849/149‑1 (MQ=255)
                                                                                                                                            gAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCtat  <  1:46805/149‑1 (MQ=255)
                                                                                                                                                tAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGt        >  1:524820/1‑140 (MQ=255)
                                                                                                                                                tAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGt        <  2:524820/140‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCTAT  >  NC_000913/222813‑223102

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑