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)
            gCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCtt                                                                                                                              <  1:479852/149‑1 (MQ=255)
               ataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgc                                                                                                                           >  1:719702/1‑149 (MQ=255)
                 aaaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgcgc                                                                                                                         >  1:113111/1‑149 (MQ=255)
                 aaaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgcgc                                                                                                                         <  2:1152544/149‑1 (MQ=255)
                 aaaaGGTGGCGAAGAGCGTACCCGCAATTTTTATTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgcgc                                                                                                                         <  2:326546/149‑1 (MQ=255)
                   aaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtg                                                                                                                       <  2:175165/149‑1 (MQ=255)
                   aaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtg                                                                                                                       >  2:1053600/1‑149 (MQ=255)
                    aGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACTGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtgg                                                                                                                      >  2:363764/1‑149 (MQ=255)
                     ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggt                                                                                                                     <  2:199825/149‑1 (MQ=255)
                     ggTGGCGAAGAGCGTACACGCAATTTTTCTTGACCGTGATGGCAGCATTAATGTCGATCTCTGCTATGTCCATGAGATCGACATCTTTGAATTTATCGCCGGTGTTATTGACGCCATGCGCGATCT‑ATAAAAATGGGCTTTGCGCttgt                                                                                                                     <  2:300267/149‑1 (MQ=255)
                      gTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTTCGCtggtg                                                                                                                    >  2:582429/1‑149 (MQ=255)
                      gTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggtg                                                                                                                    <  1:188489/149‑1 (MQ=255)
                       gggCGAAGCGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggtgg                                                                                                                   <  2:436555/148‑1 (MQ=255)
                       gcgCGAAGACCGTACCCGCAAGTCTGCTTGACCGTCACGCCAACATAAATGTCGAACACGGCTTTGTCCAAGAGATCGACAACCTTGAATATCTCGACGGTGTTATTGCCGCCAGGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtgttgg                                                                                                                   <  2:694854/147‑1 (MQ=255)
                         gctaatagCGTAACCGCAATTTTTCTTGACGTTGATTGCACCATTAATGTAGATCACGGTTATGTCCATGAGATCGACATCTTTGAATTTTTCGATGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                 <  2:1206221/143‑1 (MQ=255)
                         gCGAAGAGCGTACCCGGAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAGTTATCGTCGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                 <  1:732002/149‑1 (MQ=255)
                         gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                 <  2:879995/149‑1 (MQ=255)
                          cGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGGTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtag                                                                                                                <  1:885660/149‑1 (MQ=255)
                          cGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTGCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtag                                                                                                                <  1:733030/149‑1 (MQ=255)
                           gAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtagt                                                                                                               >  2:998300/1‑149 (MQ=255)
                           gAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGGGgtagt                                                                                                               <  2:473380/149‑1 (MQ=255)
                                 cGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACTACTTTGAATTTTTCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAATAAATGGGCTTTGCGCTGGTGGTAGTaacgaa                                                                                                         <  2:54897/149‑1 (MQ=255)
                                  gTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                        <  2:753346/149‑1 (MQ=255)
                                  gTACCCGCAATTTTTCTTGACCGTCATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAATTATCGACGGTGTTATTGCCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                        <  2:667147/149‑1 (MQ=255)
                                  gTACCCGCAATTTTTCTTGACCGTCATCGCACCAATAATGTCGATCACGGCTATGTCCATCAGAGCGACAACTTTGAAATTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATCGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                        <  2:1087367/149‑1 (MQ=255)
                                  gTACCAGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                        <  1:1083251/149‑1 (MQ=255)
                                      ccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTATCCAACCAgg                                                                                                    >  2:320821/1‑148 (MQ=255)
                                       cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATAGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTc                                                                                                   >  1:729229/1‑149 (MQ=255)
                                        gCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                  >  1:167349/1‑149 (MQ=255)
                                        gCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGGTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                  >  2:956689/1‑149 (MQ=255)
                                            tttttCTTGACCGTGATGGCACCATTAATGTCGATCACTGCTATGTCCTTGAGATCGAAAAGTTTGAATTTATCGACGGTGTTATTGACGCCATGCGGGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCa                                                                                              <  1:679841/149‑1 (MQ=255)
                                                  ttGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTc                                                                                        >  2:605439/1‑149 (MQ=255)
                                                  ttGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTc                                                                                        >  2:230495/1‑149 (MQ=255)
                                                     aCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcgcg                                                                                     <  1:606263/149‑1 (MQ=255)
                                                     aCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcgcg                                                                                     <  2:530300/149‑1 (MQ=255)
                                                      ccGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATCGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCgg                                                                                    <  1:566923/149‑1 (MQ=255)
                                                         tGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaa                                                                                 >  2:21235/1‑149 (MQ=255)
                                                          gATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaaa                                                                                >  2:873500/1‑149 (MQ=255)
                                                              gCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTa                                                                            >  1:1003116/1‑149 (MQ=255)
                                                               cACGATTAATGTCTATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTAACGCCATGCGCGAGCT‑AAAAAACTGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATAGCTCGCGGTAAATTTAc                                                                           <  1:760795/149‑1 (MQ=255)
                                                                aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACACCATGCGCGAGC‑AAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAcc                                                                          <  1:34631/149‑1 (MQ=255)
                                                                     tAATGTCGATCACGGCTATGTCCATGAGATCGAAAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGc                                                                     >  2:963068/1‑149 (MQ=255)
                                                                       aTGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGcac                                                                   >  2:503943/1‑149 (MQ=255)
                                                                        tGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCCCTGGTGGTAGAAACCAACCAGACTGGCATTGATCGCGGTAAATTTACGGAAGccca                                                                  >  1:237525/1‑149 (MQ=255)
                                                                         gTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑‑AAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGt                                                                >  1:401536/1‑149 (MQ=255)
                                                                         gTCGATCACGGCTATGTCCATGAGATCGACAACTTAGAAATTATCGACGGTGTTATTGACGCCATGCGCGAGC‑AAAAAAAATGGTCTTTGCGCTGGTGGTAGTAACCAACCACTCTCGCCTTCCTCGCGGTAAATTTACCGAAGCACAg                                                                 >  2:712959/1‑149 (MQ=255)
                                                                          tCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGt                                                                >  2:825567/1‑149 (MQ=255)
                                                                             aTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTg                                                             >  1:201593/1‑149 (MQ=255)
                                                                             aTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACAAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTg                                                             >  1:373258/1‑149 (MQ=255)
                                                                                aCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                          <  1:413171/149‑1 (MQ=255)
                                                                                aCGGCTATGGCAATGAGATCGACAACTTTGAATTTAGCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                          <  2:808513/149‑1 (MQ=255)
                                                                                aCGGATATGTCCATTAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAATAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGTCATTGCTCGCGGTAAATTTACCGAAGCACAGTTGGaaa                                                          <  1:960379/149‑1 (MQ=255)
                                                                                   gCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTCTTTACGCCATGCGAGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTATCCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGc                                                       <  1:682282/149‑1 (MQ=255)
                                                                                   gCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGc                                                       <  1:15780/149‑1 (MQ=255)
                                                                                    cTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTAGCGCTGGTGGTATTAACCAACAAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCt                                                      >  1:1064329/1‑149 (MQ=255)
                                                                                       tGTCCATGAGTTCGTCAACTTTGAATTTATCGACGGTGTTATTGCCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAc                                                   <  2:237525/149‑1 (MQ=255)
                                                                                        gTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                  <  1:981878/149‑1 (MQ=255)
                                                                                         tCCAAGAGAACCACAACTATCAATATATCGACGGAGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGCGCATTGAGCTTGTGGTAGTAACCAACCAGACTCGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACCCTGAcag                                                 >  2:639904/1‑147 (MQ=255)
                                                                                          ccATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGa                                                <  1:1053600/149‑1 (MQ=255)
                                                                                           cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGCCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAg                                               >  2:120278/1‑149 (MQ=255)
                                                                                            aTGAGATCGCCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCTCGGTAAATTTACCTAAGCACAGTTTGAAACGCTGCCCGAGt                                              >  2:1115577/1‑149 (MQ=255)
                                                                                            aTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGt                                              <  1:477848/149‑1 (MQ=255)
                                                                                             tGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATTGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtg                                             <  2:270422/149‑1 (MQ=255)
                                                                                              gagaTCGACAACTTTGAAGTTATCTACGGTGTTATTGACGACATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGTTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACAGAAGCACAGATTGAAACGCTGACCGAGtgg                                            <  2:1021631/149‑1 (MQ=255)
                                                                                                 ttCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatg                                         <  1:224447/148‑1 (MQ=255)
                                                                                                 aTCGACATCTTTGAATTTATCGACGGTGTTATTTACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatg                                         <  2:331573/149‑1 (MQ=255)
                                                                                                   cGACAACTTTGAAAGTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCAGAGTTTGAACCGCTGACCGAGtggatgga                                       <  1:1115577/149‑1 (MQ=255)
                                                                                                    gACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGGGGATGGAc                                      >  1:584735/1‑149 (MQ=255)
                                                                                                     aCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCACCCAGTCTGGCATTGGTCGCGGTAAAGTAAACGAAGCACACTTTGAAACGCTTACCGCGTGGATGGACt                                     >  2:263150/1‑149 (MQ=255)
                                                                                                     aCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACt                                     >  2:894952/1‑149 (MQ=255)
                                                                                                      cAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTg                                    <  1:854719/149‑1 (MQ=255)
                                                                                                                   aTCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGAc                       >  2:405996/1‑149 (MQ=255)
                                                                                                                     cGAAGGTGTTATTGACGCCATGCGCGAGAT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGACTCGCATTGCTCGCCGTAAATTTACCGAATCACAGATTGAAACGCTGACCTAGTGGATGGCCACGTCGATGGCGGAccc                     >  1:261635/1‑148 (MQ=255)
                                                                                                                      gACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCga                    <  2:25337/149‑1 (MQ=255)
                                                                                                                       aCGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCgag                   <  1:632353/149‑1 (MQ=255)
                                                                                                                              tATTGATGCCCTGCGCGAGCT‑AAAAATAGGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGa            <  2:93288/149‑1 (MQ=255)
                                                                                                                               aTTGACGCCATGCGCTAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACTGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt           <  1:320821/149‑1 (MQ=255)
                                                                                                                               aTTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt           >  1:334208/1‑149 (MQ=255)
                                                                                                                                ttGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGTGTGGAGGGACTGGTCGCTGGCGGACCGAGATGTCGATc          <  1:825567/149‑1 (MQ=255)
                                                                                                                                ttGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATc          >  1:227416/1‑149 (MQ=255)
                                                                                                                                   aCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTgg       >  1:103652/1‑149 (MQ=255)
                                                                                                                                     gCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt     <  2:167349/149‑1 (MQ=255)
                                                                                                                                      ccATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCAGAGTTTGAAACGCTGAACGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg    <  2:729229/149‑1 (MQ=255)
                                                                                                                                      ccATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg    >  2:941947/1‑149 (MQ=255)
                                                                                                                                      ccATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg    <  1:339392/149‑1 (MQ=255)
                                                                                                                                      ccATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCCGGCGGACCGAGATGTCGATCTGGATg    <  1:363764/149‑1 (MQ=255)
                                                                                                                                       cATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAA‑CAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGt  >  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: ‑