Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 222,965 (A)8→7 coding (133/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.97.7% 378.6 / 1.3 87coding (126/576 nt)gmhBD,D‑heptose 1,7‑bisphosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (36/49);  total (37/50)
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

CCTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGT  >  NC_000913/222811‑223096
                                                                                                                                                   |                                                                                                                                          
ccTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaa                                                                                                                                           >  2:339392/1‑149 (MQ=255)
 cTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGCTGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGCCAACTTTGAATATATCGACGGTGGTATTGACGCCATGCGCGAGCTaaa                                                                                                                                          >  1:1098232/1‑149 (MQ=255)
  ttattaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCTCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaa                                                                                                                                         <  1:1152632/149‑1 (MQ=255)
  ttattaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTAGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaa                                                                                                                                         >  1:802593/1‑149 (MQ=255)
    attaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACATCTTTGAATTTATCGACGGTGTTATTGAAGCCATGCGCGAGCTaaaaaa                                                                                                                                       >  2:826135/1‑149 (MQ=255)
    attaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaa                                                                                                                                       <  1:433157/149‑1 (MQ=255)
     ttaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaa                                                                                                                                      >  2:477848/1‑149 (MQ=255)
        aGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGTCTAAGTCCATGAGATCTACAACTTTGAATTTATCGACGGTGTAATTGACGCCATGCGCGACCTaaaaaaatgg                                                                                                                                   >  2:679841/1‑146 (MQ=255)
            gCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCtt                                                                                                                              <  1:479852/149‑1 (MQ=255)
               ataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgc                                                                                                                           >  1:719702/1‑149 (MQ=255)
                 aaaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgcgc                                                                                                                         >  1:113111/1‑149 (MQ=255)
                 aaaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgcgc                                                                                                                         <  2:1152544/149‑1 (MQ=255)
                 aaaaGGTGGCGAAGAGCGTACCCGCAATTTTTATTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgcgc                                                                                                                         <  2:326546/149‑1 (MQ=255)
                   aaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtg                                                                                                                       <  2:175165/149‑1 (MQ=255)
                   aaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtg                                                                                                                       >  2:1053600/1‑149 (MQ=255)
                    aGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACTGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtgg                                                                                                                      >  2:363764/1‑149 (MQ=255)
                     ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggt                                                                                                                     <  2:199825/149‑1 (MQ=255)
                     ggTGGCGAAGAGCGTACACGCAATTTTTCTTGACCGTGATGGCAGCATTAATGTCGATCTCTGCTATGTCCATGAGATCGACATCTTTGAATTTATCGCCGGTGTTATTGACGCCATGCGCGATCTATAAAAATGGGCTTTGCGCttgt                                                                                                                     <  2:300267/149‑1 (MQ=255)
                      gTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTTCGCtggtg                                                                                                                    >  2:582429/1‑149 (MQ=255)
                      gTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggtg                                                                                                                    <  1:188489/149‑1 (MQ=255)
                       gggCGAAGCGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggtgg                                                                                                                   <  2:436555/148‑1 (MQ=255)
                       gcgCGAAGACCGTACCCGCAAGTCTGCTTGACCGTCACGCCAACATAAATGTCGAACACGGCTTTGTCCAAGAGATCGACAACCTTGAATATCTCGACGGTGTTATTGCCGCCAGGCGCGAGCTAAAAAAATGGGCTTTGCGCtgttgg                                                                                                                   <  2:694854/147‑1 (MQ=255)
                         gctaatagCGTAACCGCAATTTTTCTTGACGTTGATTGCACCATTAATGTAGATCACGGTTATGTCCATGAGATCGACATCTTTGAATTTTTCGATGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                 <  2:1206221/143‑1 (MQ=255)
                         gCGAAGAGCGTACCCGGAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAGTTATCGTCGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                 <  1:732002/149‑1 (MQ=255)
                         gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                 <  2:879995/149‑1 (MQ=255)
                          cGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGGTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtag                                                                                                                <  1:885660/149‑1 (MQ=255)
                          cGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTGCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtag                                                                                                                <  1:733030/149‑1 (MQ=255)
                           gAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtagt                                                                                                               >  2:998300/1‑149 (MQ=255)
                           gAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGGGgtagt                                                                                                               <  2:473380/149‑1 (MQ=255)
                                 cGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACTACTTTGAATTTTTCGACGGTGTTATTGACGCCATGCGCGAGCTAAATAAATGGGCTTTGCGCTGGTGGTAGTaacgaa                                                                                                         <  2:54897/149‑1 (MQ=255)
                                  gTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                        <  2:753346/149‑1 (MQ=255)
                                  gTACCCGCAATTTTTCTTGACCGTCATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAATTATCGACGGTGTTATTGCCGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                        <  2:667147/149‑1 (MQ=255)
                                  gTACCCGCAATTTTTCTTGACCGTCATCGCACCAATAATGTCGATCACGGCTATGTCCATCAGAGCGACAACTTTGAAATTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATCGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                        <  2:1087367/149‑1 (MQ=255)
                                  gTACCAGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                        <  1:1083251/149‑1 (MQ=255)
                                      ccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTATCCAACCAgg                                                                                                    >  2:320821/1‑148 (MQ=255)
                                       cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATAGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTc                                                                                                   >  1:729229/1‑149 (MQ=255)
                                        gCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                  >  1:167349/1‑149 (MQ=255)
                                        gCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGGTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                  >  2:956689/1‑149 (MQ=255)
                                            tttttCTTGACCGTGATGGCACCATTAATGTCGATCACTGCTATGTCCTTGAGATCGAAAAGTTTGAATTTATCGACGGTGTTATTGACGCCATGCGGGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCa                                                                                              <  1:679841/149‑1 (MQ=255)
                                                  ttGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTc                                                                                        >  2:605439/1‑149 (MQ=255)
                                                  ttGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTc                                                                                        >  2:230495/1‑149 (MQ=255)
                                                     aCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcgcg                                                                                     <  1:606263/149‑1 (MQ=255)
                                                     aCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcgcg                                                                                     <  2:530300/149‑1 (MQ=255)
                                                      ccGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATCGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCgg                                                                                    <  1:566923/149‑1 (MQ=255)
                                                         tGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaa                                                                                 >  2:21235/1‑149 (MQ=255)
                                                          gATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaaa                                                                                >  2:873500/1‑149 (MQ=255)
                                                              gCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTa                                                                            >  1:1003116/1‑149 (MQ=255)
                                                               cACGATTAATGTCTATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTAACGCCATGCGCGAGCTAAAAAACTGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATAGCTCGCGGTAAATTTAc                                                                           <  1:760795/149‑1 (MQ=255)
                                                                aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACACCATGCGCGAGCAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAcc                                                                          <  1:34631/149‑1 (MQ=255)
                                                                     tAATGTCGATCACGGCTATGTCCATGAGATCGAAAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGc                                                                     >  2:963068/1‑149 (MQ=255)
                                                                       aTGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGcac                                                                   >  2:503943/1‑149 (MQ=255)
                                                                        tGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCCCTGGTGGTAGAAACCAACCAGACTGGCATTGATCGCGGTAAATTTACGGAAGccca                                                                  >  1:237525/1‑149 (MQ=255)
                                                                         gTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑‑AAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGt                                                                >  1:401536/1‑149 (MQ=255)
                                                                         gTCGATCACGGCTATGTCCATGAGATCGACAACTTAGAAATTATCGACGGTGTTATTGACGCCATGCGCGAGCAAAAAAAATGGTCTTTGCGCTGGTGGTAGTAACCAACCACTCTCGCCTTCCTCGCGGTAAATTTACCGAAGCACAg                                                                 >  2:712959/1‑149 (MQ=255)
                                                                          tCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGt                                                                >  2:825567/1‑149 (MQ=255)
                                                                             aTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTg                                                             >  1:201593/1‑149 (MQ=255)
                                                                             aTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACAAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTg                                                             >  1:373258/1‑149 (MQ=255)
                                                                                aCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                          <  1:413171/149‑1 (MQ=255)
                                                                                aCGGCTATGGCAATGAGATCGACAACTTTGAATTTAGCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                          <  2:808513/149‑1 (MQ=255)
                                                                                aCGGATATGTCCATTAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAATAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGTCATTGCTCGCGGTAAATTTACCGAAGCACAGTTGGaaa                                                          <  1:960379/149‑1 (MQ=255)
                                                                                   gCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTCTTTACGCCATGCGAGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTATCCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGc                                                       <  1:682282/149‑1 (MQ=255)
                                                                                   gCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGc                                                       <  1:15780/149‑1 (MQ=255)
                                                                                    cTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTAGCGCTGGTGGTATTAACCAACAAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCt                                                      >  1:1064329/1‑149 (MQ=255)
                                                                                       tGTCCATGAGTTCGTCAACTTTGAATTTATCGACGGTGTTATTGCCGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAc                                                   <  2:237525/149‑1 (MQ=255)
                                                                                        gTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                  <  1:981878/149‑1 (MQ=255)
                                                                                         tCCAAGAGAACCACAACTATCAATATATCGACGGAGTTATTGACGCCATGCGCGAGCTAAAAAAATGCGCATTGAGCTTGTGGTAGTAACCAACCAGACTCGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACCCTGAcag                                                 >  2:639904/1‑147 (MQ=255)
                                                                                          ccATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGa                                                <  1:1053600/149‑1 (MQ=255)
                                                                                           cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGCCGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAg                                               >  2:120278/1‑149 (MQ=255)
                                                                                            aTGAGATCGCCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCTCGGTAAATTTACCTAAGCACAGTTTGAAACGCTGCCCGAGt                                              >  2:1115577/1‑149 (MQ=255)
                                                                                            aTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGt                                              <  1:477848/149‑1 (MQ=255)
                                                                                             tGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATTGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtg                                             <  2:270422/149‑1 (MQ=255)
                                                                                              gagaTCGACAACTTTGAAGTTATCTACGGTGTTATTGACGACATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGTTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACAGAAGCACAGATTGAAACGCTGACCGAGtgg                                            <  2:1021631/149‑1 (MQ=255)
                                                                                                 ttCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatg                                         <  1:224447/148‑1 (MQ=255)
                                                                                                 aTCGACATCTTTGAATTTATCGACGGTGTTATTTACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatg                                         <  2:331573/149‑1 (MQ=255)
                                                                                                   cGACAACTTTGAAAGTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCAGAGTTTGAACCGCTGACCGAGtggatgga                                       <  1:1115577/149‑1 (MQ=255)
                                                                                                    gACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGGGGATGGAc                                      >  1:584735/1‑149 (MQ=255)
                                                                                                     aCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCACCCAGTCTGGCATTGGTCGCGGTAAAGTAAACGAAGCACACTTTGAAACGCTTACCGCGTGGATGGACt                                     >  2:263150/1‑149 (MQ=255)
                                                                                                     aCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACt                                     >  2:894952/1‑149 (MQ=255)
                                                                                                      cAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTg                                    <  1:854719/149‑1 (MQ=255)
                                                                                                                   aTCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGAc                       >  2:405996/1‑149 (MQ=255)
                                                                                                                     cGAAGGTGTTATTGACGCCATGCGCGAGATAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGACTCGCATTGCTCGCCGTAAATTTACCGAATCACAGATTGAAACGCTGACCTAGTGGATGGCCACGTCGATGGCGGAccc                     >  1:261635/1‑148 (MQ=255)
                                                                                                                      gACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCga                    <  2:25337/149‑1 (MQ=255)
                                                                                                                       aCGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCgag                   <  1:632353/149‑1 (MQ=255)
                                                                                                                              tATTGATGCCCTGCGCGAGCTAAAAATAGGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGa            <  2:93288/149‑1 (MQ=255)
                                                                                                                               aTTGACGCCATGCGCTAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACTGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt           <  1:320821/149‑1 (MQ=255)
                                                                                                                               aTTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt           >  1:334208/1‑149 (MQ=255)
                                                                                                                                ttGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGTGTGGAGGGACTGGTCGCTGGCGGACCGAGATGTCGATc          <  1:825567/149‑1 (MQ=255)
                                                                                                                                ttGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATc          >  1:227416/1‑149 (MQ=255)
                                                                                                                                   aCGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTgg       >  1:103652/1‑149 (MQ=255)
                                                                                                                                     gCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt     <  2:167349/149‑1 (MQ=255)
                                                                                                                                      ccATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCAGAGTTTGAAACGCTGAACGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg    <  2:729229/149‑1 (MQ=255)
                                                                                                                                      ccATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg    >  2:941947/1‑149 (MQ=255)
                                                                                                                                      ccATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg    <  1:339392/149‑1 (MQ=255)
                                                                                                                                      ccATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCCGGCGGACCGAGATGTCGATCTGGATg    <  1:363764/149‑1 (MQ=255)
                                                                                                                                       cATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGt  >  1:415567/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
CCTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGT  >  NC_000913/222811‑223096

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: