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.92.5% 335.1 / 23.0 80coding (126/576 nt)gmhBD,D‑heptose 1,7‑bisphosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base . (28/46);  total (31/49)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
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.

AGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGA  >  NC_000913/222819‑223092
                                                                                                                                           |                                                                                                                                      
aGGAGCTATAAACGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTa                                                                                                                                        <  1:128294/140‑1 (MQ=255)
aGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGTTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTa                                                                                                                                        <  2:222521/140‑1 (MQ=255)
     cTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATTGCACCAATAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaa                                                                                                                                   >  2:1282339/1‑140 (MQ=255)
     cTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaa                                                                                                                                   >  2:119339/1‑140 (MQ=255)
     cTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaa                                                                                                                                   >  1:879294/1‑140 (MQ=255)
     cTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaa                                                                                                                                   >  1:723262/1‑140 (MQ=255)
     cTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaa                                                                                                                                   >  1:245042/1‑140 (MQ=255)
     cTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaa                                                                                                                                   >  1:1607704/1‑140 (MQ=255)
     cTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaa                                                                                                                                   >  1:1607247/1‑140 (MQ=255)
       ataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaat                                                                                                                                 <  1:1014807/140‑2 (MQ=255)
            cgggggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGACTTTATCGACGGTGTTATTGACGCCCTGCGCGAGCT‑AAAAAAATGGGCt                                                                                                                           <  1:1482196/136‑1 (MQ=255)
             ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCtt                                                                                                                          >  1:144503/1‑140 (MQ=255)
             ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCtt                                                                                                                          >  1:403935/1‑140 (MQ=255)
               tGGCGAAGAGCGTACCCGCAATTTTTCTTGACCTTGATGGCACCATTAATGTCGATCACGGGTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGGGCGTGCT‑AATAAAATGGGCTTTg                                                                                                                        >  2:479376/1‑140 (MQ=255)
                ggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgc                                                                                                                       >  1:277838/1‑140 (MQ=255)
                 gCGAAGAGCGTCCCCGCAATTTGTCTTGACCGTGATGGCACCATTAATGTCGATCACGGATATGTCCAAGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgcg                                                                                                                      <  2:1675419/140‑1 (MQ=255)
                 gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgcg                                                                                                                      <  1:1404895/140‑1 (MQ=255)
                 gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgcg                                                                                                                      <  1:77424/140‑1 (MQ=255)
                 gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgcg                                                                                                                      <  1:77444/140‑1 (MQ=255)
                    aaGTGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGTTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtg                                                                                                                   <  2:128954/140‑1 (MQ=255)
                    aaGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTATTGTCTATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGATCT‑AAAAAAATGGGCTTTGCGCtg                                                                                                                   <  1:1203778/140‑1 (MQ=255)
                     agagCGTACCCGCAATTTTTCTTGACCGTTATGGCACATTTACTGTCGATCAAGGCTATGTCCATGTTATAGTGAAATTTGAATTTCTCGTCGGTGTTATTGACGCCTTGCGCGAGCT‑AAAAAAATGGGCTGTGCGCtgg                                                                                                                  <  1:1303687/140‑1 (MQ=255)
                     agagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACTCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtgg                                                                                                                  <  2:1522721/140‑1 (MQ=255)
                     agagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtgg                                                                                                                  <  2:1469513/140‑1 (MQ=255)
                     agagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtgg                                                                                                                  <  2:1469419/140‑1 (MQ=255)
                         cGTACCCTCAATTTTTCTTGCCCGAGATGGAACCCTAACTTTCGATCACGGCTATGTCCATGAGATCGACAACTTTGACTTTATCGACGGTGTTTTTGACGCCATGCTCGATCT‑AAAAAATTGGGCTTTGCGCTggtggt                                                                                                              <  1:1636148/140‑1 (MQ=255)
                         cGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTggtggt                                                                                                              <  2:531058/140‑1 (MQ=255)
                           tACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtag                                                                                                            >  2:1326917/1‑140 (MQ=255)
                           tACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtag                                                                                                            >  2:1291466/1‑140 (MQ=255)
                                   aTTTTTCTTGCCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCCTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  1:284147/140‑1 (MQ=255)
                                   aTTTTTCTTGACGGTGATGGCACCATTAATGTCGATCACGGCTATGTCCTTGAGTTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  2:1330630/140‑1 (MQ=255)
                                   aTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  2:70540/140‑1 (MQ=255)
                                   aTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  2:723262/140‑1 (MQ=255)
                                   aTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  1:525813/140‑1 (MQ=255)
                                   aTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGAGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  2:308061/140‑1 (MQ=255)
                                   aTTTTTCTTGACCGTGATGGCACCATTAATGGCGCTCACGGCTATGTCCATGAGATCGCCTACTTTGTATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  2:1490302/140‑1 (MQ=255)
                                   aTTTTCCTTGCCCGTGAAGGCACCACTAATCTCGAGTCCGTCAATGTCCATGAGAACCAAAACGTAGAATTAATAGACGGTGTTATTGCCGCCTTGCGCGAG‑AAAAAAAACTGTGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  2:986516/140‑1 (MQ=255)
                                      tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAaccaacca                                                                                                  <  1:1545800/140‑1 (MQ=255)
                                        tCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTc                                                                                               <  2:1306600/140‑1 (MQ=255)
                                        tCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTc                                                                                               <  2:1306114/140‑1 (MQ=255)
                                         cTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCAAGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAA‑GGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                              >  2:70364/1‑140 (MQ=255)
                                                        aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGt                                                                               <  2:755280/140‑1 (MQ=255)
                                                                 gTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTGcc                                                                      >  1:337364/1‑140 (MQ=255)
                                                                  tCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCg                                                                     <  1:1107164/140‑1 (MQ=255)
                                                                    gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGa                                                                    >  1:1145640/1‑140 (MQ=255)
                                                                    gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGaa                                                                   >  1:9199/1‑140 (MQ=255)
                                                                    gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGaa                                                                   >  1:316489/1‑140 (MQ=255)
                                                                    gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATATACCGaa                                                                   >  1:821/1‑140 (MQ=255)
                                                                     aTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAg                                                                  >  1:177151/1‑140 (MQ=255)
                                                                     aTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAg                                                                  >  1:496170/1‑140 (MQ=255)
                                                                       cACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑ACAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGca                                                                <  1:1500581/140‑1 (MQ=255)
                                                                       cACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCTTTGCTCGCGGTAAATTTACCGAAGca                                                                <  2:635377/140‑1 (MQ=255)
                                                                       cACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGca                                                                <  1:520431/140‑1 (MQ=255)
                                                                          ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCGGGCATTGCTCGCGGTAAATTTACCGAAGCACAg                                                             >  2:833018/1‑140 (MQ=255)
                                                                          ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGGTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCGGGCATTGCTCGCGGTAAATTTACCGAAGCACAg                                                             >  2:279667/1‑140 (MQ=255)
                                                                            cTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGtt                                                           >  2:1080475/1‑140 (MQ=255)
                                                                               tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGa                                                        >  1:153723/1‑140 (MQ=255)
                                                                               tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGa                                                        >  2:1505879/1‑140 (MQ=255)
                                                                                   cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGGTGGTGGTAGTTACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAATCg                                                    <  2:959658/140‑1 (MQ=255)
                                                                                   cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGGTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGTAACg                                                    <  2:959640/140‑1 (MQ=255)
                                                                                    aTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGc                                                   >  2:772628/1‑140 (MQ=255)
                                                                                    aTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGc                                                   >  2:1390646/1‑140 (MQ=255)
                                                                                      gagaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAAACAGTCTGGCATTGCTCGCGGTAAATTTACCGAATCACAGTTTGAAACGCTg                                                 >  1:561606/1‑140 (MQ=255)
                                                                                      gagaTCGACAACTTTGAAGTTATCGTCGGTGGTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTGTGCGCTGGTGGTAGTAACCCACCAGTCTGGCAGTGCTCGCGGTAAATTTACCGAAGCACAGTTGGACACGCTg                                                 >  2:388723/1‑140 (MQ=255)
                                                                                       agaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGa                                                <  2:19090/140‑1 (MQ=255)
                                                                                         aTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAc                                               >  1:519817/1‑140 (MQ=255)
                                                                                          tCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCg                                             <  1:423994/140‑1 (MQ=255)
                                                                                          tCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCg                                             <  2:1275272/140‑1 (MQ=255)
                                                                                           cGACAACTTTTAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTTACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGa                                            <  1:1522127/140‑1 (MQ=255)
                                                                                           cGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCg                                             <  1:1587522/140‑1 (MQ=255)
                                                                                           cGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGa                                            <  1:1013319/140‑1 (MQ=255)
                                                                                           cGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGa                                            <  1:1472872/140‑1 (MQ=255)
                                                                                             aCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGt                                          >  2:189268/1‑140 (MQ=255)
                                                                                             aCAACTTTGAATATATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGAAAATATACCGAAGCACAGTTTGAAACGCAGGCCGAGt                                          >  2:1597579/1‑140 (MQ=255)
                                                                                                    tGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatgga                                   <  1:810565/140‑1 (MQ=255)
                                                                                                     gAATTTATCGTCGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGAc                                  <  1:709129/140‑1 (MQ=255)
                                                                                                         ttATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTCACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGt                              <  1:1274496/140‑1 (MQ=255)
                                                                                                           aTCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGTCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCg                            >  2:439427/1‑140 (MQ=255)
                                                                                                           aTCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCg                            >  2:439538/1‑140 (MQ=255)
                                                                                                                   tgtTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGa                    >  1:144857/1‑140 (MQ=255)
                                                                                                                       aTTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAGTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAGCGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCga                <  2:350975/140‑1 (MQ=255)
                                                                                                                        ttGACGCCATTCGCGAGAT‑AAAAAATTGGGCTTTGCGCTGTTGGTAGTAACCAACCTGTCTGGCATTGCTCGCGGTAATTTTACCGAAGCACAGTTTGAAACGCTGACAGAGTGGATGGACTGGTCGCTGGCGGACCgag               <  1:1282339/140‑1 (MQ=255)
                                                                                                                             gCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTc          >  1:1292093/1‑140 (MQ=255)
                                                                                                                                aTGCGCGAGCT‑AAAATAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGTAGCACAGTTTGTTTCGCTGACCGAGTGGATGGTCTGGTCGCTGGCTGTACGATATGtgttt       <  1:119339/140‑5 (MQ=255)
                                                                                                                                aTGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt       <  2:1205141/140‑1 (MQ=255)
                                                                                                                                aTGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt       <  2:1590576/140‑1 (MQ=255)
                                                                                                                                 tGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATc      >  2:1351858/1‑140 (MQ=255)
                                                                                                                                    gcgAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTgg   <  2:177151/140‑1 (MQ=255)
                                                                                                                                     cgAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGa  <  1:1612049/140‑1 (MQ=255)
                                                                                                                                     cgAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGa  <  1:1294611/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                      
AGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGA  >  NC_000913/222819‑223092

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑