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.98.9% 392.1 / ‑2.6 88coding (126/576 nt)gmhBD,D‑heptose 1,7‑bisphosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (38/49);  total (39/49)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCTA  >  NC_000913/222815‑223101
                                                                                                                                               |                                                                                                                                               
attaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaa                                                                                                                                            >  1:815161/1‑149 (MQ=255)
attaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTCACCGTGATGGCACCATTAATGTCGATCACGGCTTTGTCGAGGAGATCGCCAACTTTGAAGTTATCGACGGTGTTTTTGACGCCCTGCGCGAGCTaaaaaa                                                                                                                                            >  2:94990/1‑149 (MQ=255)
   aaGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatg                                                                                                                                         <  1:426826/149‑3 (MQ=255)
     ggAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatggg                                                                                                                                       >  2:486279/1‑145 (MQ=255)
        gCTATAAAAGGGGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGGCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCtt                                                                                                                                   <  2:724028/149‑1 (MQ=255)
         cTATTAAAGGTGGCGTAGAGCTTACCCGCAATTTTTCTTGACCGTGATGTCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCttt                                                                                                                                  <  2:249972/149‑1 (MQ=255)
           ataCAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTATTGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgc                                                                                                                                <  1:1281618/149‑1 (MQ=255)
           ataAAAGGTGGCGAAGAGCGTACCCGCAAGTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgc                                                                                                                                >  1:32844/1‑149 (MQ=255)
            taAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgcg                                                                                                                               >  1:952121/1‑149 (MQ=255)
            taAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgcg                                                                                                                               >  2:365345/1‑149 (MQ=255)
             aaaaGGTGGCGTAGAGCGTACCCGCAATTTTTCTTGTCCGTGATGGCACCATTTATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgcgc                                                                                                                              <  1:904298/149‑1 (MQ=255)
             aaaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTTACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgcgc                                                                                                                              <  2:1186266/149‑1 (MQ=255)
                atgTGGCGAAGAGCGTACCCGCAATTTTTGTTGACCGTGTTGGCACCATTAATGTCGTTCACGGCTATGTCCATGAGTTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtgg                                                                                                                           <  1:124022/147‑1 (MQ=255)
                 ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggt                                                                                                                          <  2:212965/149‑1 (MQ=255)
                  gTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggtg                                                                                                                         >  2:1255602/1‑149 (MQ=255)
                   tGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggtgg                                                                                                                        <  2:573863/149‑1 (MQ=255)
                   tGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggtgg                                                                                                                        <  2:363622/149‑1 (MQ=255)
                   tGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggtgg                                                                                                                        <  1:457366/149‑1 (MQ=255)
                     ccGAAGCGCGTACCCGCAATTTTTCTTCACCGTGATGGCACCATTAATGTCGATCACGGCTATGCCCATGAGAGCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                      <  1:540225/148‑1 (MQ=255)
                       gAAGAGCTTACCCTAAATTTTTCTTTACCGTGATTGCACCATTAAAGTCGATCACGGCTATGTCAATGAGATCGACAACTTTGAATTTATCGACGGTGGTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtagt                                                                                                                    <  2:526448/149‑1 (MQ=255)
                       gAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACAATTAATGTCTATCACGGCTATGTCCATGAGATCGACAACTTTGAATATATCGACGGTGTTATTGAAGCCATGCGCTAGCTAAAAAAATGGTCTTTGCGCTGGTGctagt                                                                                                                    >  2:743270/1‑149 (MQ=255)
                         agagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaa                                                                                                                  <  1:580956/149‑1 (MQ=255)
                          gagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaac                                                                                                                 <  2:300358/149‑1 (MQ=255)
                              gTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                             <  2:1200279/149‑1 (MQ=255)
                                  ccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGt                                                                                                         >  1:661243/1‑149 (MQ=255)
                                   cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTc                                                                                                        >  2:450503/1‑149 (MQ=255)
                                   cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTc                                                                                                        >  2:340098/1‑149 (MQ=255)
                                     cAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTg                                                                                                      <  2:650068/149‑1 (MQ=255)
                                     cAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGAAATGTCCATGAGATCCACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTg                                                                                                      <  1:998511/149‑1 (MQ=255)
                                        tttttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCa                                                                                                   >  1:1294769/1‑149 (MQ=255)
                                         ttttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAt                                                                                                  >  1:427098/1‑149 (MQ=255)
                                          tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAtt                                                                                                 <  1:788350/149‑1 (MQ=255)
                                                gACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTAGCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcgc                                                                                           <  1:704691/149‑1 (MQ=255)
                                                   cGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGt                                                                                        >  1:641329/1‑149 (MQ=255)
                                                   cGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACAAACAAGTCTGGCATTGCTCGCGGt                                                                                        >  2:1182550/1‑149 (MQ=255)
                                                          gCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTa                                                                                 <  2:1294769/149‑1 (MQ=255)
                                                          gCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTa                                                                                 <  1:808333/149‑1 (MQ=255)
                                                          gCACCATTAATGTCGAGCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGCCGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTa                                                                                 <  1:55061/149‑1 (MQ=255)
                                                           cACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTAATCGACGGTGGTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAc                                                                                <  2:1003869/149‑1 (MQ=255)
                                                            aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAcc                                                                               >  1:576039/1‑149 (MQ=255)
                                                               aTTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGaa                                                                            <  2:661243/149‑1 (MQ=255)
                                                                 tAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATATATCGACGGTCTAATTGACGCCATGCGCGAGCTAATAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGATCGCGGTAAATTTACCGAAGc                                                                          >  2:788832/1‑149 (MQ=255)
                                                                        gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGTTGGTGGTAGTAACCAACCAGTCTGGCATTGCTc                                                                                             <  2:617440/123‑1 (MQ=255)
                                                                        gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGTTGGTGGTAGTAACCAACCAGTCTGGCATTGCTc                                                                                             >  1:617440/1‑123 (MQ=255)
                                                                        gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGttt                                                                   <  1:120664/149‑1 (MQ=255)
                                                                        gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTAGCGGTAAATTTACCGAAGAACTGttt                                                                   >  1:554932/1‑149 (MQ=255)
                                                                        gATCACGGCTATGTCCATGAGATCGACAAATTTGAATTTATCGCCGGTGTTATTGACGCCATGCGCGATATAAAAAAATTGGCTTTGCGCTGGTGGGAGTACCCAACAAGTCTGGCATTGCTCGCGGTAAAATTAAAGAAGCCCAGgtt                                                                   >  2:425340/1‑149 (MQ=255)
                                                                           cACGGCTATGTCCATGAGATCGACAACTTTGAATTTATAGACGGTGTTATTGACGCAAAGCGCGAGCTAAAAAAATGTGCTTTGAGCTGGTGGTTGTAACCACCCAGTCTGGCATTGCTCGCGGTAAATTTAGCGAAGCAAAGTTTGaa                                                                >  1:855529/1‑149 (MQ=255)
                                                                            aCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                               <  2:1162820/149‑1 (MQ=255)
                                                                              ggatATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACg                                                             <  1:130855/146‑1 (MQ=255)
                                                                              ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACg                                                             >  1:440950/1‑149 (MQ=255)
                                                                              ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACg                                                             >  2:189579/1‑149 (MQ=255)
                                                                              ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACg                                                             >  2:173793/1‑149 (MQ=255)
                                                                                      ccATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGTGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGa                                                     >  1:390689/1‑149 (MQ=255)
                                                                                       cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTCACCCACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCGGCCCGAg                                                    >  2:928215/1‑149 (MQ=255)
                                                                                       cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAg                                                    >  2:22223/1‑149 (MQ=255)
                                                                                         tGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtg                                                  <  1:365345/149‑1 (MQ=255)
                                                                                          gagaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCACCCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtgg                                                 >  2:597855/1‑149 (MQ=255)
                                                                                            taTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggat                                               <  1:1054756/148‑1 (MQ=255)
                                                                                                taCAACTTTGAATTTTTCGACTGTGTTCTTGTCGCCATGCGCGAGCTAAAAATATGGGCTTTGCGCTGGTGGTAGTAAAAAACCAGTCTGGCATTTCTCGCTGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGAc                                           <  1:743270/148‑1 (MQ=255)
                                                                                                gACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTATCCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAAAGCTGACCGAGTGGATGGAc                                           >  1:697670/1‑149 (MQ=255)
                                                                                                 aCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGCTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACt                                          >  1:169241/1‑149 (MQ=255)
                                                                                                      tttGAATTTAACGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCg                                     >  2:541745/1‑149 (MQ=255)
                                                                                                          aaTTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGCAACGCTGCCCGAGTGGATGGCCTGGTCGCTgg                                 >  2:1044416/1‑149 (MQ=255)
                                                                                                          aaTTTATCGACGGTGTTATTGACGCCATGCGCGAGCAAAAAAAATGGGCTTTGCGCTCGTGGTAGTAACCAACCAGTATGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCtcg                                 >  1:775085/1‑147 (MQ=255)
                                                                                                           aTTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGc                                >  1:307721/1‑149 (MQ=255)
                                                                                                           aTTTATCGACGGTGTTATTGACGACATGCGCGAGCTAAAAAAATGTGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGTTCGCTGGc                                >  1:369025/1‑149 (MQ=255)
                                                                                                            tttATCGACGGTGTTATTGACGGCATGCGCGAGCTTATAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACTGTTTGAAACGCTGTCCGAGTGGATGGACTGGTCGCTGGCg                               <  1:357255/149‑1 (MQ=255)
                                                                                                              tATCGACTGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAATCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGa                             <  1:414692/149‑1 (MQ=255)
                                                                                                               aTCGACTGTGTTATTGACGCCATGCGCGAGCTAAAAAAATTGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGTTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGAc                            <  2:1255453/149‑1 (MQ=255)
                                                                                                               aTCGACGGTGTTATGGACGCCCTGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTTGCATTGCTCGCGGTAAATTTACCGAAGCACACTTTGACACGCTGCCCGAGTGGATGGACTGGTCGCTGGCGGAc                            <  2:101173/149‑1 (MQ=255)
                                                                                                               aTCGACGGTGATATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGAc                            <  1:215722/149‑1 (MQ=255)
                                                                                                                tCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAGCCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAGAGAGTCGATGGACTGGTCGCTGGCGGAcc                           >  2:1148340/1‑149 (MQ=255)
                                                                                                                tCGAAGGTGTTATTGACGCCATGCGCGAGCTAAAAAATTCGGCTTTGCGCCGGTGGTAGTAACCAAGCAGTCTGGAATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGATGACCGATTGGATGGACTGGTCGCTGGCGGAcc                           <  2:277320/149‑1 (MQ=255)
                                                                                                                 cttcGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCGCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCg                          <  2:115154/146‑1 (MQ=255)
                                                                                                                   aCGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGGGGTAGTAACCAACCAGGCTGGCATTGCTCGGGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGGGGATGGACTGGTCGCTGGCGGACCgag                        <  2:603330/149‑1 (MQ=255)
                                                                                                                     ggTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGGCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGAt                      <  2:745890/149‑1 (MQ=255)
                                                                                                                       ttttATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTTACCAACCAGTCTGGCATTGCTCGCTGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGTAAGTACTGGTCGCTGGCGGACCGAGATGt                    <  1:475337/147‑1 (MQ=255)
                                                                                                                         ttATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCg                  >  1:1122159/1‑149 (MQ=255)
                                                                                                                            ttGCCGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTACCCACCCAGTCTGGGATTGCTCGCGGTAAATTTCCCGATGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTGGATc               >  2:697380/1‑149 (MQ=255)
                                                                                                                            ttGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATc               <  2:427098/149‑1 (MQ=255)
                                                                                                                              gACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCTCTGGGGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACAGAAGCACAGTTTGAACCTGTGCCCGAGTGGATGGACTGTTCGATTGCGGACCTAGATGTCGAGCtt             >  2:1286574/1‑148 (MQ=255)
                                                                                                                              gACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTg             >  2:403849/1‑149 (MQ=255)
                                                                                                                               aCGCCATGCGCGAGCTAAAAAAATGGGCTTTCCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCTCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGAGGATGGACTGTTCGCTGGCGGACAGAGATATCGATCTgg            >  1:312946/1‑149 (MQ=255)
                                                                                                                                cGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGa           <  1:486279/149‑1 (MQ=255)
                                                                                                                                 gCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt          >  2:960190/1‑149 (MQ=255)
                                                                                                                                 gCCATGCGCGAGCTAAAAAAATGGGCGTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt          <  2:32844/149‑1 (MQ=255)
                                                                                                                                  ccATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg         <  1:649147/149‑1 (MQ=255)
                                                                                                                                    aTGCGCGAGCTAAAAAAATGGGCTTTGCGCTTGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGACGGt       <  2:377149/149‑1 (MQ=255)
                                                                                                                                     tGCGCGAGCTAAAAAAATGGGCTTTGCGCGGGTGGTAGTCCCCAACCAGGCTGGCATGGCTCGCGGTAAATTTACCGACGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTa      <  2:108561/149‑1 (MQ=255)
                                                                                                                                       cgcgAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATc    <  2:576039/149‑1 (MQ=255)
                                                                                                                                         cgAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTACCCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAACCACAGTTTGAAACGCTGACCGAGTTCATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCta  <  2:139803/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
ATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCTA  >  NC_000913/222815‑223101

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: