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.96.9% 412.6 / 5.9 97coding (126/576 nt)gmhBD,D‑heptose 1,7‑bisphosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base . (48/46);  total (50/47)
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
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 106 total aligned reads displayed.

CCTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCT  >  NC_000913/222811‑223100
                                                                                                                                                   |                                                                                                                                              
ccTTATTAAGCAGCAAAAAAACCTGCCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGCCGGTGTTATTGACGCCATGCGCGAGCTaa                                                                                                                                               <  2:779152/149‑1 (MQ=255)
   tattaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATTGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTCTTGACGCCATGCGCGAGCTaaaaa                                                                                                                                            <  1:429108/149‑1 (MQ=255)
   tattaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaa                                                                                                                                            >  2:781406/1‑149 (MQ=255)
   tattaAGGAGCGATAAAAGGGGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATCGCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaa                                                                                                                                            <  2:1007547/149‑1 (MQ=255)
    attaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGAGAGCTaaaaaa                                                                                                                                           >  1:1218638/1‑149 (MQ=255)
      taAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaat                                                                                                                                         >  2:79035/1‑148 (MQ=255)
        aGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatgg                                                                                                                                       >  1:1421690/1‑146 (MQ=255)
         ggAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatggg                                                                                                                                      >  1:191192/1‑145 (MQ=255)
               ataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgc                                                                                                                               >  1:508015/1‑149 (MQ=255)
               ataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgc                                                                                                                               >  1:1146839/1‑149 (MQ=255)
                taAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgcg                                                                                                                              >  1:1368700/1‑149 (MQ=255)
                taAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgcc                                                                                                                              >  1:1086851/1‑148 (MQ=255)
                 aaaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgcgc                                                                                                                             <  1:800463/149‑1 (MQ=255)
                      gTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggtg                                                                                                                        <  1:786144/149‑1 (MQ=255)
                       tGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGATATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggtgg                                                                                                                       >  2:1103090/1‑149 (MQ=255)
                        ggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGCCACCTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAACAATGGGCTTTGCGCTggtggt                                                                                                                      >  1:1534028/1‑149 (MQ=255)
                        ggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTggtggt                                                                                                                      <  2:191192/149‑1 (MQ=255)
                         gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTggtggtt                                                                                                                     >  1:151470/1‑148 (MQ=255)
                         gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                     >  1:1174687/1‑149 (MQ=255)
                          cGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtag                                                                                                                    >  2:341461/1‑149 (MQ=255)
                           gAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtagt                                                                                                                   >  1:693043/1‑149 (MQ=255)
                            aaGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtagta                                                                                                                  <  1:803911/149‑1 (MQ=255)
                            aaGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtagta                                                                                                                  <  1:396335/149‑1 (MQ=255)
                             agagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaa                                                                                                                 <  2:1212843/149‑1 (MQ=255)
                              gagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaac                                                                                                                >  1:802346/1‑149 (MQ=255)
                              gagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaac                                                                                                                >  1:596291/1‑149 (MQ=255)
                              gagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaac                                                                                                                >  1:1167193/1‑149 (MQ=255)
                                    aCCCGCAATTTTTCTTGACCGGGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAaccaacca                                                                                                          <  2:828598/149‑1 (MQ=255)
                                     cccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATAGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAATAATGGGCTTTGCGCTGGTGGTAGTACCCATCCAg                                                                                                         >  2:1030668/1‑149 (MQ=255)
                                        gCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                      >  2:225828/1‑149 (MQ=255)
                                          aaTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTCGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTgg                                                                                                    <  1:405347/149‑1 (MQ=255)
                                          aaTTTGTCTTGACCGTGACGTCAACATTAATGTCGCTCCCGGCTATGTCCATGAGATCGACAACTTTCAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTgg                                                                                                    <  2:1230806/149‑1 (MQ=255)
                                           aTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGc                                                                                                   >  2:334788/1‑149 (MQ=255)
                                           aTTTTTCGTGCCCGGGATGGCACCCTTAACGGCGATCCCGGCTATGTCCATGAGCACGCCAACTTTGACTTAATCGACGGTGTTATTCACGCCATCCGCGAGCT‑AACAAAATCGGCAATTCGCTCGAGGTAGTAACCAACCAGTCTGGc                                                                                                   <  2:1181378/149‑1 (MQ=255)
                                              tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAtt                                                                                                >  2:40418/1‑149 (MQ=255)
                                                   tGACCGTGATTGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGATATTGACGCCATGCGCGAGCT‑AAAATAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcg                                                                                           <  2:60019/149‑1 (MQ=255)
                                                   tGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcg                                                                                           >  2:734310/1‑149 (MQ=255)
                                                    gACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGC‑AAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcgc                                                                                          >  2:699850/1‑149 (MQ=255)
                                                     aCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGAGCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATTGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcgcg                                                                                         <  2:1471490/149‑1 (MQ=255)
                                                       cGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGCCGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGt                                                                                       >  1:1497719/1‑149 (MQ=255)
                                                         tGATGGCCCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCCGTCTGGCATTTCTCGCGGTaa                                                                                     >  2:390951/1‑149 (MQ=255)
                                                          gATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaaa                                                                                    <  2:1402344/149‑1 (MQ=255)
                                                          gATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaaa                                                                                    <  2:1477329/149‑1 (MQ=255)
                                                           aTGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGCCAACTTTGAATTTATCGACGGTGTTATTGACGACAAGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTACTAACCAACCAGTCTGGCATTGCTCCCGGTaact                                                                                   >  2:1244708/1‑147 (MQ=255)
                                                               cACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAc                                                                               <  1:24080/149‑1 (MQ=255)
                                                               cACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAc                                                                               <  1:1459977/149‑1 (MQ=255)
                                                                 ccATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTTCGCTGGTGGTTGTAACCCACCAGTCTGGCATTGCTCGCTGTAAATTTACCg                                                                             >  1:117878/1‑149 (MQ=255)
                                                                   aTTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTTAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTATCCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGaa                                                                           >  1:993524/1‑149 (MQ=255)
                                                                   aTTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGaa                                                                           >  2:1107198/1‑149 (MQ=255)
                                                                      aatggagATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGca                                                                        <  2:300429/143‑1 (MQ=255)
                                                                       aTGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGcac                                                                       >  2:511716/1‑149 (MQ=255)
                                                                        tGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGcaca                                                                      <  2:1069636/149‑1 (MQ=255)
                                                                         gTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGAAGCAATGCGCGAGCT‑AAAAAAATGTGCTTTGAGATGGTGGTAGTAACAAACCAGACTGGCAATGCTAGCGGTAAATATACAGAAGCACAg                                                                     >  2:470264/1‑149 (MQ=255)
                                                                            gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGttt                                                                  <  1:837890/149‑1 (MQ=255)
                                                                             aTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTg                                                                 <  1:79631/149‑1 (MQ=255)
                                                                                aCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                              <  1:341461/149‑1 (MQ=255)
                                                                                  ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACCGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACg                                                            >  2:1243380/1‑149 (MQ=255)
                                                                                     tATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTGATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTACCCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTg                                                         <  1:719265/149‑1 (MQ=255)
                                                                                        gTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                      >  1:34569/1‑149 (MQ=255)
                                                                                        gTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGC‑AAAAAAAATGGGCTATGCGCTGGTGGTAGTAACCAACCAGTCTGGAATTGCTCGAGGTAAATTAACCGAAGCACACATTGAAACGCTGAcc                                                      >  1:1139953/1‑149 (MQ=255)
                                                                                          ccATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGa                                                    >  1:416410/1‑149 (MQ=255)
                                                                                            aTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGt                                                  >  1:414431/1‑149 (MQ=255)
                                                                                            aTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGt                                                  >  2:727337/1‑149 (MQ=255)
                                                                                            aTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAGTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGt                                                  >  1:858436/1‑149 (MQ=255)
                                                                                             tGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTCGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtg                                                 >  2:1184432/1‑149 (MQ=255)
                                                                                              gagaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtgg                                                >  2:300470/1‑149 (MQ=255)
                                                                                                gaTCGACAACTTTGATTTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggat                                              <  1:880105/149‑1 (MQ=255)
                                                                                                 aTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGCGtggatg                                             <  2:1421690/149‑1 (MQ=255)
                                                                                                   cGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatgga                                           <  1:781406/149‑1 (MQ=255)
                                                                                                    gACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGAc                                          >  2:227131/1‑149 (MQ=255)
                                                                                                     aCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCTACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACt                                         <  2:209500/149‑1 (MQ=255)
                                                                                                       aaCTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGTCTTTGCTCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTgg                                       <  2:596291/149‑1 (MQ=255)
                                                                                                          tttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAACCGCTGACCGAGTGGATGGACTGGTCg                                    >  2:1462610/1‑149 (MQ=255)
                                                                                                          tttGAATTTATCGACGGGGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTCCGCTGGTGGTAGAAACCAACCAGTCTGGCATTGCACGCGGTAAAATTACCGAACCACAGTAGTAAACCCTGACCAAGTGCATGGACTGGACg                                    >  2:1301732/1‑149 (MQ=255)
                                                                                                           ttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGc                                   <  2:887275/149‑1 (MQ=255)
                                                                                                              aaTTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTgg                                <  1:1287872/149‑1 (MQ=255)
                                                                                                                tttATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCg                              <  1:551259/149‑1 (MQ=255)
                                                                                                                  tATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGa                            <  2:1146839/149‑1 (MQ=255)
                                                                                                                  tATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGGGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGa                            <  2:508015/149‑1 (MQ=255)
                                                                                                                      gACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAAACAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCCAGTGGATGGACTGCTCGCTGGCGGACAga                        >  2:1187233/1‑149 (MQ=255)
                                                                                                                         ggTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCgagct                     >  2:766487/1‑147 (MQ=255)
                                                                                                                         ggTGTTAATGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGGATTGCTCGCGGTAAATTAACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGAt                     <  1:281272/149‑1 (MQ=255)
                                                                                                                           tgtTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGt                   <  1:79035/149‑1 (MQ=255)
                                                                                                                           tgtTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTCGTGGTAGTAACCAACCAGTCTGGCATTGCACGCGCTAAATTTACCCAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGACATGt                   >  1:1086848/1‑149 (MQ=255)
                                                                                                                           tgtTATTGACGCCATGCGAGAGCT‑AAAAAAATTGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGt                   <  2:1218638/149‑1 (MQ=255)
                                                                                                                               aTTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt               >  2:267520/1‑149 (MQ=255)
                                                                                                                               aTTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt               >  1:825580/1‑149 (MQ=255)
                                                                                                                                 tgccGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGAGGGACTGGTCGCTGGCGGACCGAGATGTCGATCt             <  1:1030668/146‑1 (MQ=255)
                                                                                                                                 tGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCt             >  2:353308/1‑149 (MQ=255)
                                                                                                                                 tGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCt             <  2:802346/149‑1 (MQ=255)
                                                                                                                                 tGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAGCt             <  2:686683/149‑1 (MQ=255)
                                                                                                                                  gACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTg            >  1:749798/1‑149 (MQ=255)
                                                                                                                                   aCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCTAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTgg           >  2:185400/1‑149 (MQ=255)
                                                                                                                                    cGCCATGCGCGAGCT‑AAAAACATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGGTTGATACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGa          <  2:157854/149‑1 (MQ=255)
                                                                                                                                    cGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGa          <  1:300470/149‑1 (MQ=255)
                                                                                                                                    cGCCATGCGCGAGC‑AAAATACATGGGCTTTGCGCTGGGGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTTAATTTCCCGAAGCACAGTTTGAACCGCGGACCGAGGGGATGGACTGGGCGCTGGCGGACCGAGATCTCGATCTGGa          <  1:63805/149‑1 (MQ=255)
                                                                                                                                       cATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATgg       <  1:334788/149‑1 (MQ=255)
                                                                                                                                          gcgcGAGCT‑AAAAAAATGGGCTTTGCGCTGGGGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTAt    <  2:1174687/149‑1 (MQ=255)
                                                                                                                                           cgcgAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAACCGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATCTCGATCTGGATGGTATc   >  1:1412996/1‑149 (MQ=255)
                                                                                                                                            gcgAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCt  <  1:225828/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CCTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCT  >  NC_000913/222811‑223100

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: ‑