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)
            cgggggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGACTTTATCGACGGTGTTATTGACGCCCTGCGCGAGCTAAAAAAATGGGCt                                                                                                                           <  1:1482196/136‑1 (MQ=255)
             ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCtt                                                                                                                          >  1:144503/1‑140 (MQ=255)
             ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCtt                                                                                                                          >  1:403935/1‑140 (MQ=255)
               tGGCGAAGAGCGTACCCGCAATTTTTCTTGACCTTGATGGCACCATTAATGTCGATCACGGGTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGGGCGTGCTAATAAAATGGGCTTTg                                                                                                                        >  2:479376/1‑140 (MQ=255)
                ggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgc                                                                                                                       >  1:277838/1‑140 (MQ=255)
                 gCGAAGAGCGTCCCCGCAATTTGTCTTGACCGTGATGGCACCATTAATGTCGATCACGGATATGTCCAAGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgcg                                                                                                                      <  2:1675419/140‑1 (MQ=255)
                 gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgcg                                                                                                                      <  1:1404895/140‑1 (MQ=255)
                 gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgcg                                                                                                                      <  1:77424/140‑1 (MQ=255)
                 gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgcg                                                                                                                      <  1:77444/140‑1 (MQ=255)
                    aaGTGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGTTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtg                                                                                                                   <  2:128954/140‑1 (MQ=255)
                    aaGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTATTGTCTATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGATCTAAAAAAATGGGCTTTGCGCtg                                                                                                                   <  1:1203778/140‑1 (MQ=255)
                     agagCGTACCCGCAATTTTTCTTGACCGTTATGGCACATTTACTGTCGATCAAGGCTATGTCCATGTTATAGTGAAATTTGAATTTCTCGTCGGTGTTATTGACGCCTTGCGCGAGCTAAAAAAATGGGCTGTGCGCtgg                                                                                                                  <  1:1303687/140‑1 (MQ=255)
                     agagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACTCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtgg                                                                                                                  <  2:1522721/140‑1 (MQ=255)
                     agagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtgg                                                                                                                  <  2:1469513/140‑1 (MQ=255)
                     agagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtgg                                                                                                                  <  2:1469419/140‑1 (MQ=255)
                         cGTACCCTCAATTTTTCTTGCCCGAGATGGAACCCTAACTTTCGATCACGGCTATGTCCATGAGATCGACAACTTTGACTTTATCGACGGTGTTTTTGACGCCATGCTCGATCTAAAAAATTGGGCTTTGCGCTggtggt                                                                                                              <  1:1636148/140‑1 (MQ=255)
                         cGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTggtggt                                                                                                              <  2:531058/140‑1 (MQ=255)
                           tACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtag                                                                                                            >  2:1326917/1‑140 (MQ=255)
                           tACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtag                                                                                                            >  2:1291466/1‑140 (MQ=255)
                                   aTTTTTCTTGCCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCCTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  1:284147/140‑1 (MQ=255)
                                   aTTTTTCTTGACGGTGATGGCACCATTAATGTCGATCACGGCTATGTCCTTGAGTTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  2:1330630/140‑1 (MQ=255)
                                   aTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  2:70540/140‑1 (MQ=255)
                                   aTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  2:723262/140‑1 (MQ=255)
                                   aTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  1:525813/140‑1 (MQ=255)
                                   aTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGAGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  2:308061/140‑1 (MQ=255)
                                   aTTTTTCTTGACCGTGATGGCACCATTAATGGCGCTCACGGCTATGTCCATGAGATCGCCTACTTTGTATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  2:1490302/140‑1 (MQ=255)
                                   aTTTTCCTTGCCCGTGAAGGCACCACTAATCTCGAGTCCGTCAATGTCCATGAGAACCAAAACGTAGAATTAATAGACGGTGTTATTGCCGCCTTGCGCGAGAAAAAAAACTGTGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                    <  2:986516/140‑1 (MQ=255)
                                      tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAaccaacca                                                                                                  <  1:1545800/140‑1 (MQ=255)
                                        tCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTc                                                                                               <  2:1306600/140‑1 (MQ=255)
                                        tCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTc                                                                                               <  2:1306114/140‑1 (MQ=255)
                                         cTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCAAGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAAGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                              >  2:70364/1‑140 (MQ=255)
                                                        aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGt                                                                               <  2:755280/140‑1 (MQ=255)
                                                                 gTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTGcc                                                                      >  1:337364/1‑140 (MQ=255)
                                                                  tCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCg                                                                     <  1:1107164/140‑1 (MQ=255)
                                                                    gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGa                                                                    >  1:1145640/1‑140 (MQ=255)
                                                                    gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGaa                                                                   >  1:9199/1‑140 (MQ=255)
                                                                    gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGaa                                                                   >  1:316489/1‑140 (MQ=255)
                                                                    gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATATACCGaa                                                                   >  1:821/1‑140 (MQ=255)
                                                                     aTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAg                                                                  >  1:177151/1‑140 (MQ=255)
                                                                     aTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAg                                                                  >  1:496170/1‑140 (MQ=255)
                                                                       cACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTACAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGca                                                                <  1:1500581/140‑1 (MQ=255)
                                                                       cACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCTTTGCTCGCGGTAAATTTACCGAAGca                                                                <  2:635377/140‑1 (MQ=255)
                                                                       cACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGca                                                                <  1:520431/140‑1 (MQ=255)
                                                                          ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCGGGCATTGCTCGCGGTAAATTTACCGAAGCACAg                                                             >  2:833018/1‑140 (MQ=255)
                                                                          ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGGTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCGGGCATTGCTCGCGGTAAATTTACCGAAGCACAg                                                             >  2:279667/1‑140 (MQ=255)
                                                                            cTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGtt                                                           >  2:1080475/1‑140 (MQ=255)
                                                                               tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGa                                                        >  1:153723/1‑140 (MQ=255)
                                                                               tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGa                                                        >  2:1505879/1‑140 (MQ=255)
                                                                                   cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGGTGGTGGTAGTTACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAATCg                                                    <  2:959658/140‑1 (MQ=255)
                                                                                   cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGGTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGTAACg                                                    <  2:959640/140‑1 (MQ=255)
                                                                                    aTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGc                                                   >  2:772628/1‑140 (MQ=255)
                                                                                    aTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGc                                                   >  2:1390646/1‑140 (MQ=255)
                                                                                      gagaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAAACAGTCTGGCATTGCTCGCGGTAAATTTACCGAATCACAGTTTGAAACGCTg                                                 >  1:561606/1‑140 (MQ=255)
                                                                                      gagaTCGACAACTTTGAAGTTATCGTCGGTGGTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTGTGCGCTGGTGGTAGTAACCCACCAGTCTGGCAGTGCTCGCGGTAAATTTACCGAAGCACAGTTGGACACGCTg                                                 >  2:388723/1‑140 (MQ=255)
                                                                                       agaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGa                                                <  2:19090/140‑1 (MQ=255)
                                                                                         aTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAc                                               >  1:519817/1‑140 (MQ=255)
                                                                                          tCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCg                                             <  1:423994/140‑1 (MQ=255)
                                                                                          tCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCg                                             <  2:1275272/140‑1 (MQ=255)
                                                                                           cGACAACTTTTAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTTACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGa                                            <  1:1522127/140‑1 (MQ=255)
                                                                                           cGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCg                                             <  1:1587522/140‑1 (MQ=255)
                                                                                           cGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGa                                            <  1:1013319/140‑1 (MQ=255)
                                                                                           cGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGa                                            <  1:1472872/140‑1 (MQ=255)
                                                                                             aCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGt                                          >  2:189268/1‑140 (MQ=255)
                                                                                             aCAACTTTGAATATATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGAAAATATACCGAAGCACAGTTTGAAACGCAGGCCGAGt                                          >  2:1597579/1‑140 (MQ=255)
                                                                                                    tGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatgga                                   <  1:810565/140‑1 (MQ=255)
                                                                                                     gAATTTATCGTCGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGAc                                  <  1:709129/140‑1 (MQ=255)
                                                                                                         ttATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTCACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGt                              <  1:1274496/140‑1 (MQ=255)
                                                                                                           aTCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGTCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCg                            >  2:439427/1‑140 (MQ=255)
                                                                                                           aTCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCg                            >  2:439538/1‑140 (MQ=255)
                                                                                                                   tgtTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGa                    >  1:144857/1‑140 (MQ=255)
                                                                                                                       aTTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAGTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAGCGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCga                <  2:350975/140‑1 (MQ=255)
                                                                                                                        ttGACGCCATTCGCGAGATAAAAAATTGGGCTTTGCGCTGTTGGTAGTAACCAACCTGTCTGGCATTGCTCGCGGTAATTTTACCGAAGCACAGTTTGAAACGCTGACAGAGTGGATGGACTGGTCGCTGGCGGACCgag               <  1:1282339/140‑1 (MQ=255)
                                                                                                                             gCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTc          >  1:1292093/1‑140 (MQ=255)
                                                                                                                                aTGCGCGAGCTAAAATAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGTAGCACAGTTTGTTTCGCTGACCGAGTGGATGGTCTGGTCGCTGGCTGTACGATATGtgttt       <  1:119339/140‑5 (MQ=255)
                                                                                                                                aTGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt       <  2:1205141/140‑1 (MQ=255)
                                                                                                                                aTGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt       <  2:1590576/140‑1 (MQ=255)
                                                                                                                                 tGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATc      >  2:1351858/1‑140 (MQ=255)
                                                                                                                                    gcgAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTgg   <  2:177151/140‑1 (MQ=255)
                                                                                                                                     cgAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGa  <  1:1612049/140‑1 (MQ=255)
                                                                                                                                     cgAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGa  <  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: