Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,364,998 C→A 6.4% intergenic (+19/+20) nudI → / ← ais nucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,364,9980CA 6.4% 711.3 / 25.9 203intergenic (+19/+20)nudI/aisnucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (108/82);  new base A (8/5);  total (116/87)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-01

Only 100 of 203 total aligned reads displayed.

CAGCTGCTTTTGACAGAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGATTTACCTGATTTTTGACTGCGTTTCTGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACGATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACC  >  NC_000913/2364734‑2365271
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                 
cAGCTGCTTTTGACAGAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGATTTACCTGATTTTTGACTGCGTTTCTGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTCCTGGCTGTCAg                                                                                                                                                                                                                                                                         >  1:191893/1‑275 (MQ=255)
             cAGAAAGCACGCCGTGGAACTTCAGCGATGATATTCGAAACAAGACGTAGGCGGATGGTCGCAAGGAAGAGATTTATATGATTTACCTGATTTTTGACTGCGGTGCTGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAAttt                                                                                                                                                                                                                                                              <  1:22685/273‑1 (MQ=255)
                        ccGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGATTTACCTGATTTTTGACTGCGTTTCTGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGaa                                                                                                                                                                                                                                                 <  2:315453/275‑1 (MQ=255)
                                                     cAAGACGTATGCAGATTGTGGCAAGGANNAGAGTTATATGATTTACCTGATTTTTTACGNCGTTTCTTCCAAACGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGGGGGTAAAACCTGAAGATCTGGNGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCttt                                                                                                                                                                                                                      <  3:80985/273‑1 (MQ=255)
                                                                                agagATTTATATGATTTACCTGATTTTTGACTGCGTTTCTGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGAATGACTAAACCATCCAGAt                                                                                                                                                                                         >  1:100374/1‑275 (MQ=255)
                                                                                         tatGATTTACCTGATTTTTGACTGCGTTTCTGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGt                                                                                                                                                                                 >  1:77337/1‑274 (MQ=255)
                                                                                              tttACCTGATTTTTGACTGCGTTTCTGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAAt                                                                                                                                                                           >  2:626650/1‑275 (MQ=255)
                                                                                                    tggtttttGACTGCGTTTCTGCCAAACGANNAGTTAAAATAAACGAAGAGTTTCAGGACTCNGGGTGGGTAAAAACTGAAGAGCTTGTGGATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGNCTTCTGTAACAATTATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcgcg                                                                                                                                                                     <  3:100374/272‑1 (MQ=255)
                                                                                                                        gCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                                                                                                                                                                                                               >  7:820450/1‑149 (MQ=255)
                                                                                                                           aaCCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGt                                                                                                                                                                                                                                                                            >  6:1024057/1‑149 (MQ=255)
                                                                                                                           aaCCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATCATCAGGTTTAAATGTCGCGTCCCGCTTAACTTTAGCAatata                                                                                                                                              >  2:264793/1‑275 (MQ=255)
                                                                                                                            aCCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTc                                                                                                                                                                                                                                                                           <  7:569430/149‑1 (MQ=255)
                                                                                                                            aCCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTAGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTc                                                                                                                                                                                                                                                                           <  7:168915/149‑1 (MQ=255)
                                                                                                                             ccGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCa                                                                                                                                                                                                                                                                          <  7:1014329/149‑1 (MQ=255)
                                                                                                                              cGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAg                                                                                                                                                                                                                                                                         >  7:94182/1‑149 (MQ=255)
                                                                                                                                agaAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAt                                                                                                                                                                                                                                                                       <  5:752726/149‑1 (MQ=255)
                                                                                                                                agaAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCa                                                                                                                                         >  2:533446/1‑275 (MQ=255)
                                                                                                                                     tGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtttt                                                                                                                                                                                                                     >  4:456639/1‑194 (MQ=255)
                                                                                                                                          aTAAACGAAGAGTTTCAGGACTACGCGTGGGTATAACCTGATGATCTGGTGGATTATGATTTGAATGTCGCCAACCGAAAAACGTTACAGTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCGTCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGAtt                                                                                                                               <  2:537688/275‑1 (MQ=255)
                                                                                                                                          aTAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTa                                                                                                                                                                                                                                                             <  6:608401/149‑1 (MQ=255)
                                                                                                                                          aTAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTa                                                                                                                                                                                                                                                             <  5:995670/149‑1 (MQ=255)
                                                                                                                                              aCGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaaa                                                                                                                                                                                                                                                                 >  8:785504/1‑141 (MQ=255)
                                                                                                                                              aCGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGt                                                                                                                                                                                                                                                                            <  7:540365/130‑1 (MQ=255)
                                                                                                                                              aCGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATCATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaaa                                                                                                                                                                                                                                                                 <  7:785504/141‑1 (MQ=255)
                                                                                                                                               cGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaa                                                                                                                                                                                                                                                                   >  4:340294/1‑138 (MQ=255)
                                                                                                                                                     gTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGa                                                                                                                                                                                                                                                  >  5:935378/1‑149 (MQ=255)
                                                                                                                                                      tttCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCa                                                                                                                                     >  2:437468/1‑257 (MQ=255)
                                                                                                                                                          aGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACg                                                                                                               <  4:569515/275‑1 (MQ=255)
                                                                                                                                                          aGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACg                                                                                                               <  4:626650/275‑1 (MQ=255)
                                                                                                                                                          aGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAAAGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACg                                                                                                               <  4:135528/275‑1 (MQ=255)
                                                                                                                                                                aCGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACGATAttt                                                                                                         >  2:243675/1‑275 (MQ=255)
                                                                                                                                                                 cgcgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACGATAttt                                                                                                         >  4:149249/1‑274 (MQ=255)
                                                                                                                                                                  gcgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAg                                                                                                                                                                                                                                     >  8:558380/1‑149 (MQ=255)
                                                                                                                                                                  gcgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAg                                                                                                                                                                                                                                     >  6:113192/1‑149 (MQ=255)
                                                                                                                                                                       ggTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATaaa                                                                                                                                                                                                                                <  6:476040/149‑1 (MQ=255)
                                                                                                                                                                        gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAAc                                                                                                                                                                                                                               >  7:84947/1‑149 (MQ=255)
                                                                                                                                                                        gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAAc                                                                                                                                                                                                                               >  8:226283/1‑149 (MQ=255)
                                                                                                                                                                        gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCca                                                                                                                                                                                                                 >  4:560285/1‑163 (MQ=255)
                                                                                                                                                                         tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACt                                                                                                                                                                                                                              >  6:1076728/1‑149 (MQ=255)
                                                                                                                                                                         tAAAACCTGAAGATCTGGTGCATTATGAGTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACt                                                                                                                                                                                                                              >  7:691981/1‑149 (MQ=255)
                                                                                                                                                                           aaaCCTGAAGATCTGGTGAATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGt                                                                                                                                                                    <  3:19612/205‑1 (MQ=255)
                                                                                                                                                                            aaCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACGATATTTTTATCAGGCgct                                                                                             <  4:243675/275‑1 (MQ=255)
                                                                                                                                                                               cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCt                                                                                                                                                                                                                        >  5:485357/1‑149 (MQ=255)
                                                                                                                                                                                tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtt                                                                                                                                                                                                                       <  7:416461/149‑1 (MQ=255)
                                                                                                                                                                                  aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACGATATTTTTTCCGGGcctttgcttt                                                                                      >  2:50728/1‑275 (MQ=255)
                                                                                                                                                                                        tGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCcaca                                                                                                                                                                                                               <  5:1019841/149‑1 (MQ=255)
                                                                                                                                                                                             cATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgcat                                                                                                                                                                                                          <  8:74996/149‑1 (MQ=255)
                                                                                                                                                                                               ttATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGa                                                                                                                                                                                                        <  5:113192/149‑1 (MQ=255)
                                                                                                                                                                                               ttATGATTTGAATGTCGCAACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAAATTTGCATTTTTCCACATGCATAACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGAGTATCTTTAGCAATATATGTCAGGCAATGATTAGGGGGGAAGATAACGATATTTTTATAAGGCGCGTTTATTTGAAAGTCCtt                                                                            >  1:75266/1‑273 (MQ=255)
                                                                                                                                                                                                    aTTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGAc                                                                                                                                                                                                       >  2:56185/1‑145 (MQ=255)
                                                                                                                                                                                                     tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTaa                                                                                                                                                                                                                                                     >  5:490886/1‑98 (MQ=255)
                                                                                                                                                                                                     tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                                                                                                                                                                                                  >  7:394994/1‑149 (MQ=255)
                                                                                                                                                                                                       tGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCca                                                                                                                                                                                                                 >  2:676892/1‑132 (MQ=255)
                                                                                                                                                                                                          aTGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACGATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACtata                                                                <  2:149249/274‑1 (MQ=255)
                                                                                                                                                                                                           tGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatg                                                                                                                                                                                                             <  4:676892/132‑1 (MQ=255)
                                                                                                                                                                                                              cGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGAt                                                                                                                                                                                         >  7:324551/1‑149 (MQ=255)
                                                                                                                                                                                                               gCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGATTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATa                                                                                                                                                                                        >  8:611979/1‑149 (MQ=255)
                                                                                                                                                                                                                ccACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATaa                                                                                                                                                                                       >  7:464626/1‑149 (MQ=255)
                                                                                                                                                                                                                 cACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAAt                                                                                                                                                                                      <  5:371577/149‑1 (MQ=255)
                                                                                                                                                                                                                  aCCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATc                                                                                                                                                                                     >  6:183228/1‑149 (MQ=255)
                                                                                                                                                                                                                    ccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAg                                                                                                                                                                                   >  8:630375/1‑149 (MQ=255)
                                                                                                                                                                                                                       aaaaaCGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGtt                                                                                                                                                                                <  6:388568/149‑1 (MQ=255)
                                                                                                                                                                                                                           acgttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTaaa                                                                                                                                                                            <  8:324551/149‑1 (MQ=255)
                                                                                                                                                                                                                              ttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGt                                                                                                                                                                         >  7:407078/1‑149 (MQ=255)
                                                                                                                                                                                                                              ttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACGATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTgaag                                           >  4:88946/1‑275 (MQ=255)
                                                                                                                                                                                                                                acgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcg                                                                                                                                                                       >  8:253873/1‑149 (MQ=255)
                                                                                                                                                                                                                                acgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcg                                                                                                                                                                       <  7:670720/149‑1 (MQ=255)
                                                                                                                                                                                                                                 cgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAatatat                                                                                                                                             >  2:459423/1‑174 (MQ=255)
                                                                                                                                                                                                                                  gttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTg                                                                                                                                                                                                                                                                                 >  5:150122/1‑41 (MQ=255)
                                                                                                                                                                                                                                  gttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTg                                                                                                                                                                                                                                                                                 <  6:150122/41‑1 (MQ=255)
                                                                                                                                                                                                                                       aaaGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGGCTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACg                                                                                                                                                                >  8:825459/1‑149 (MQ=255)
                                                                                                                                                                                                                                       aaaGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACg                                                                                                                                                                >  7:435818/1‑149 (MQ=255)
                                                                                                                                                                                                                                           gTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTc                                                                                                                                                                                                                                              <  8:513213/67‑1 (MQ=255)
                                                                                                                                                                                                                                           gTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAAt                                                                                                                                                                           >  8:402972/1‑134 (MQ=255)
                                                                                                                                                                                                                                              ttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCt                                                                                                                                                         >  5:642822/1‑149 (MQ=255)
                                                                                                                                                                                                                                              ttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCt                                                                                                                                                         >  5:467285/1‑149 (MQ=255)
                                                                                                                                                                                                                                                ctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCttt                                                                                                                                                       <  7:253873/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTcc                                                                                                                                                                                                                  >  6:505524/1‑89 (MQ=255)
                                                                                                                                                                                                                                                 tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTcc                                                                                                                                                                                                                  <  5:505524/89‑1 (MQ=255)
                                                                                                                                                                                                                                                 tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTa                                                                                                                                                      <  7:207988/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTTCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAg                                                                                                                                                     >  5:332328/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  gTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAg                                                                                                                                                     <  5:823850/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     aCAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCaa                                                                                                                                                  >  5:893285/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      cAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAat                                                                                                                                                 <  8:445328/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaTGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCa                                                                                                                                                                                                                                          <  7:494943/59‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaTGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAata                                                                                                                                                >  5:850765/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          tctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAatatat                                                                                                                                             >  7:1065146/3‑149 (MQ=255)
                                                                                                                                                                                                                                                          tctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAatatat                                                                                                                                             >  8:585677/3‑149 (MQ=255)
                                                                                                                                                                                                                                                            tgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAatatat                                                                                                                                             <  6:806436/147‑1 (MQ=255)
                                                                                                                                                                                                                                                            tgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGt                                                                                                                                           >  8:71815/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            taaCAGCCATNACTTGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCAACCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGt                                                                                                                                           <  3:77250/147‑1 (MQ=255)
                                                                                                                                                                                                                                                             gaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTc                                                                                                                                          >  5:288657/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             gaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCACGTCACGCTTATCTTTAGCAATATATGTc                                                                                                                                          >  5:262989/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              aCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTTCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCa                                                                                                                                         <  6:332328/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              aCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAAt                                                                                                                                                                                      >  6:728261/1‑104 (MQ=255)
                                                                                                                                                                                                                                                                 gCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAAt                                                                                                                                                                           <  5:435869/112‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCa                                                                                                                                     >  7:563913/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   cATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCaa                                                                                                                                    <  6:829325/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                      tACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATGTGTCAGGCAATGa                                                                                                                                 <  8:182695/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aCTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACGATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAAcc  <  4:533446/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                 
CAGCTGCTTTTGACAGAAATCACGCCGTGGACCTTCAGCGATGATATTCGCACCAAGACGTATGCAGATGGTCGCAAGGAAGAGATTTATATGATTTACCTGATTTTTGACTGCGTTTCTGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACGATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACC  >  NC_000913/2364734‑2365271

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