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)
           ctttaAAAGGTGGCGAAGAGCGTACCCGAAATTTTTCTTGACCGTGATGGCTCCATTAATGTCGATCACGGCTATGGCCATGAGTTCGACAACTTTGAATTTATCGCCGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGTCttt                                                                                                                                   <  1:316404/146‑1 (MQ=255)
           cTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCttt                                                                                                                                   >  1:300437/1‑149 (MQ=255)
            tataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGCTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCttt                                                                                                                                   <  1:295948/149‑1 (MQ=255)
              taAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgcg                                                                                                                                >  1:39176/1‑149 (MQ=255)
                aaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCt                                                                                                                              >  2:750907/1‑149 (MQ=255)
                aaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCt                                                                                                                              <  2:80192/149‑1 (MQ=255)
                 aaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCt                                                                                                                              >  2:165557/1‑149 (MQ=255)
                 aaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtg                                                                                                                             >  2:720019/1‑149 (MQ=255)
                  aGGTGGCGAAGAGCGTTCCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACTACTTTGAATTTATCTACGGTGTTATTGACGCCATGCGCGAGCTTAATAAATGGGCTTTGCGCtgg                                                                                                                            >  2:1016715/1‑149 (MQ=255)
                     gggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggtgg                                                                                                                         <  1:302382/148‑1 (MQ=255)
                      ggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTggtggt                                                                                                                        >  1:723479/1‑149 (MQ=255)
                      ggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTggtggt                                                                                                                        <  2:288537/149‑1 (MQ=255)
                       gCGAAGAGCGTACCCGCAATTTTTCTTGCCCGTGATCGCACCATTAATGTCGCTCACGGCTATGTCCATTAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGTGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                       <  2:910297/149‑1 (MQ=255)
                        cGAAGAGCGTAGCCGCAATTGTGCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGGCCATGAGATCGACAACTTTGAATTTATCGACCGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTGGCGCTGGTGggag                                                                                                                      <  2:216167/149‑1 (MQ=255)
                        cGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtag                                                                                                                      >  2:302567/1‑149 (MQ=255)
                         gAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtagt                                                                                                                     >  2:487549/1‑149 (MQ=255)
                          aaGACCGTACCCGCAATTTATCATGACCGTGAAGGCACCAATAATGTCGATCCCGGCTATGTCCATGAGATCGACAACTTAGAATTTAACGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtagta                                                                                                                    <  2:714876/149‑1 (MQ=255)
                                gTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                              <  2:300437/149‑1 (MQ=255)
                                    ccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGt                                                                                                          >  1:974344/1‑149 (MQ=255)
                                     cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACAAGTc                                                                                                         >  2:639488/1‑149 (MQ=255)
                                      gCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGAAACCAACCAGTCt                                                                                                        >  2:14056/1‑149 (MQ=255)
                                       cattttttCTTGACCGTTAGGGCACCATTAATGTCGATCACGGCTATGTCCTTGAGATCGACAACTTTGATTTTATCGACGGTGTTATTGACGCCAGGCGCGAGCTAAAAAAATTGGCTTTGCGCGGGTGGTAGTAACCAACCAGTCTg                                                                                                       <  1:581503/146‑1 (MQ=255)
                                            tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAtt                                                                                                  >  2:370670/1‑149 (MQ=255)
                                            tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGAAtt                                                                                                  >  1:289112/1‑149 (MQ=255)
                                            tttCTTGACCGTGATGGCACCATTAATGTCGATCAAGGCTATGTCCATGAGATCGACAACTTTGTATTTATCGACGGTGTTATTGACTCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAtt                                                                                                  <  2:919101/149‑1 (MQ=255)
                                                ttGACCGTGATGGCACCATTAATGTCGATCACCGCTATGTACATGAGATCGACAACATTGATTTAATCGACGGTGTAATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGCTGGTAGTAACCAACCTGTCTGGCATTGCTc                                                                                              >  2:152408/1‑149 (MQ=255)
                                                 tGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTACCCCACCAGTCTGGCATTGCTcg                                                                                             >  1:936063/1‑149 (MQ=255)
                                                        gATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaaa                                                                                      >  1:845147/1‑149 (MQ=255)
                                                             ccccATTCATCCCGATCCCGGCTATGGCCCTGAGATCGACAACTAGGAATATATCGCCGGTGGTCTTGACGCCATGCGCGAGCTAAAACACTGGGCTTTGCGCTGGGGGTAGTCACCAACCCGTCTGGCAGTGCTCGCGGTAACTTTAc                                                                                 <  2:12127/147‑1 (MQ=255)
                                                             cACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAc                                                                                 <  2:873823/149‑1 (MQ=255)
                                                             cACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTAGTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGACTGGCATTGCTCGCGGTAAATTTAc                                                                                 >  2:82396/1‑149 (MQ=255)
                                                              cccATTAATGTCGATCACGGCTATGTCCCTGATAGCGACAACTTTCTATTAATCGACGGTGTTATTGACGCCAAGCGCGAGCTAAAAACATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAcc                                                                                <  2:159270/148‑1 (MQ=255)
                                                              aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAcc                                                                                >  1:885757/1‑149 (MQ=255)
                                                                cATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGa                                                                              >  1:314830/1‑149 (MQ=255)
                                                                  ttAATGTCTATCTCGGCTATGTCCTGGAGATCGACAACTTTGTATTTTTCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTGGGGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAg                                                                            <  1:676341/149‑1 (MQ=255)
                                                                  ttAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAg                                                                            <  1:134100/149‑1 (MQ=255)
                                                                  ttAATGTCGATCAAGTCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTAGCGGTAATTTTACCGAAg                                                                            >  1:86707/1‑149 (MQ=255)
                                                                   tAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGc                                                                           <  1:975978/149‑1 (MQ=255)
                                                                   tAATGTCGATCAAGGCTATGTCCATGAGATCGACAACTTTGACATTCTCGACTGTGTTATAGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGc                                                                           <  1:979330/149‑1 (MQ=255)
                                                                       gTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAg                                                                       >  2:46805/1‑149 (MQ=255)
                                                                         cGATCACGGGTATGTCCATTAGTTCGACAACTTTGAATTTATCGTCGGTGTTATTGACGCGATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCCGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGtt                                                                     <  1:263372/149‑1 (MQ=255)
                                                                            tCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATTGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATATACAGAAGCACAGTTTGa                                                                  >  1:746913/1‑149 (MQ=255)
                                                                            tCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGa                                                                  >  1:332535/1‑149 (MQ=255)
                                                                              aCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                                <  1:880692/149‑1 (MQ=255)
                                                                                ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACTGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGGCTGGCATTGCTCGCGGTAAAGTTAGCGAAGCACAGTTTGAAACg                                                              >  2:448222/1‑149 (MQ=255)
                                                                                      gTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                        >  1:103848/1‑149 (MQ=255)
                                                                                       tCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCg                                                       >  2:1007599/1‑149 (MQ=255)
                                                                                         cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAg                                                     <  1:750907/149‑1 (MQ=255)
                                                                                              gaTCGCCAACTATGAATTCATCGACGGTGTTCATGCCGCCAGGCGCGAGCGAAAAAAATGGGCTTTGCGCTGGGGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGGAACTTTACCGACGCACAGTTTGACACGCGGACCGCGtggat                                                <  2:484707/149‑1 (MQ=255)
                                                                                              gaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGAAGTAACCAACCAGACTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTCATCCGCTGACCTAGtgga                                                 >  2:506492/1‑149 (MQ=255)
                                                                                              gaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTACCCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAACCGCTGACCGAGtggat                                                >  2:948904/1‑149 (MQ=255)
                                                                                                 cGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatgga                                             <  1:514927/149‑1 (MQ=255)
                                                                                                   aCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACt                                           <  2:437027/149‑1 (MQ=255)
                                                                                                        tttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCg                                      >  2:197062/1‑149 (MQ=255)
                                                                                                        tttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCg                                      >  1:457709/1‑149 (MQ=255)
                                                                                                         ttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGc                                     <  1:327822/149‑1 (MQ=255)
                                                                                                            aaTTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTgg                                  <  2:723479/149‑1 (MQ=255)
                                                                                                               ttATCGACGGTGTTATTGACTCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGTTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCgg                               >  1:369994/1‑149 (MQ=255)
                                                                                                                             aTTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt                 <  2:974344/149‑1 (MQ=255)
                                                                                                                             aTTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt                 >  2:699728/1‑149 (MQ=255)
                                                                                                                                 aCGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTgg             >  1:515820/1‑149 (MQ=255)
                                                                                                                                   gCCATGCGCGAGCTAAAAAACTGGGCTTTGCGCTGGTGGGAGTAACCAACCAGTCTCGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATCGACTGGGCGCTGGCGGACCGAGATGTCGATCTGGAt           <  2:234371/149‑1 (MQ=255)
                                                                                                                                   gCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt           <  1:370670/149‑1 (MQ=255)
                                                                                                                                   gCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGAGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt           >  2:204124/1‑149 (MQ=255)
                                                                                                                                          gcgAGCTAAAAAAAATGGTCTTTGCTCTGGTTGTAGTAACCAACCAGTGTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATc     <  1:166849/149‑1 (MQ=255)
                                                                                                                                            gAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCtat  <  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: