Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,487,1630CT41.5% 64.1 / 24.1 53intergenic (+200/‑72)hrpA/ydcFputative ATP‑dependent helicase/DUF218 superfamily protein, SAM‑binding
Reads supporting (aligned to +/- strand):  ref base C (8/23);  new base T (22/0);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 1.43e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCA  >  NC_000913/1487016‑1487311
                                                                                                                                                    |                                                                                                                                                          
tttttATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTtcct                                                                                                                                                          <  1:345941/150‑1 (MQ=255)
 ttttATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTtcctc                                                                                                                                                         <  2:329440/150‑1 (MQ=255)
   ttattaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGATATAAGATTtcctcc                                                                                                                                                        >  1:249505/1‑149 (MQ=255)
            cTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAATAAGATTTTCTCCtttttcaggg                                                                                                                                              >  1:95200/1‑143 (MQ=255)
            cTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAAAGATTTTCTCCTTTTACAGTg                                                                                                                                              >  1:637805/1‑150 (MQ=255)
              aaaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAGGAAAAAACATTTTCTCCTTTTAGAGGGCt                                                                                                                                            >  1:665641/1‑150 (MQ=255)
              aaaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGc                                                                                                                                             <  2:635014/149‑1 (MQ=255)
              aaaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCt                                                                                                                                            <  1:11585/150‑1 (MQ=255)
              aaaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGACAAACGAATTTCTCCTTTTTGAGTg                                                                                                                                              >  1:429098/1‑148 (MQ=255)
               aaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTTTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTcc                                                                                                                                          <  1:27708/151‑1 (MQ=255)
               aaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAATACAATTTTtttctttcttggggcgcc                                                                                                                                          >  2:370526/1‑141 (MQ=255)
                 aCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGAAAAACAATTttctcc                                                                                                                                                        >  2:762224/1‑135 (MQ=255)
                      aaGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTtct                                                                                                                                    <  1:92672/150‑1 (MQ=255)
                         ggTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGATTAGGAAATACAATTTTCTTCTTTCTGAGGGCCCCGcgtttccc                                                                                                                                 >  2:512076/1‑143 (MQ=255)
                          gTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGACAAAACAATTCCTCCTTTTAAAGTGCTccccct                                                                                                                                      >  1:118419/1‑140 (MQ=255)
                                aTTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTAtct                                                                                                                         <  2:511636/151‑1 (MQ=255)
                                aTTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACAAAAGAATTTCTCCTTTTTGAGGGCTCCGCTTCTCCCtattat                                                                                                                           >  2:829528/1‑149 (MQ=255)
                                 ttCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAATACAATTTTCTCCTTTCTGAGGTGTCCGCCTCTCACtatt                                                                                                                             >  2:698646/1‑146 (MQ=255)
                                   cAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCAc                                                                                                                      <  1:298522/151‑1 (MQ=255)
                                      gAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAAAAATTTTCTTCCTTCTAAGGTGGCCGCTTCTTACGATTATCTCCCGc                                                                                                                    >  1:116082/1‑150 (MQ=255)
                                       aaTTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAg                                                                                                                 <  2:680764/151‑1 (MQ=255)
                                           cAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTAt                                                                                                             <  2:386029/150‑1 (MQ=255)
                                            aGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAAAAAATTTCTCCTTTTTTAGTGGTTCGCCTTTCCCTTTTTTTTCCCGCCGTAtt                                                                                                            >  2:664126/1‑150 (MQ=255)
                                                  aTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAgg                                                                                                    <  1:902703/151‑1 (MQ=255)
                                                   tACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGGAATAAGAATTTCTCCTTTTAGAGGGCTCCCCCTCTCCCTATTTTTTCACCCAGGATtttttatg                                                                                                    >  1:437865/1‑146 (MQ=255)
                                                     cAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGaa                                                                                                 <  2:541005/151‑1 (MQ=255)
                                                       aTTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGACATAAAAATTCCTTCTTTCTGAGTGCGCCGCCTCTCACTATTTTTTTACCCACTATTTTTTAAGGGa                                                                                                  >  2:282244/1‑149 (MQ=255)
                                                                        gAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATg                                                                              <  2:75486/150‑1 (MQ=255)
                                                                           taatGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAAAAAATTTCTTCTTTTTGAGTGGGCCGCTTCTTCCTATTATTTTACCCCGTATTTTTAAAGGGACGAAAAAGAGGAACAATGaaa                                                                           >  1:439281/1‑149 (MQ=255)
                                                                           taatGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGGTAATGACATAAGAATTCTTTCTTTCATAGTGGTCCGCCTCTCCCTTTTTTCTTACCCCATATTTTTTATGGGACGACAAAGAGGGAACATAAAc                                                                           >  2:943307/1‑150 (MQ=255)
                                                                            aatGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGACAAAAGATTTTCTCCCTTTTTAGTGCTCCCCCTCTCCCTATTATTTTACCCAGTATTTTTTAAGGAAAGAAAAAGAGGGAACCagaaa                                                                            >  1:824416/1‑145 (MQ=255)
                                                                             atGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACAt                                                                        <  1:935870/150‑1 (MQ=255)
                                                                                       aCAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTAACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGttt                                                             <  2:519850/150‑1 (MQ=255)
                                                                                        cAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAATAAAAATTTCTTCTTTTATAGGGCTCCGCTTCTCTCTATTTTTTCTCGCAGTATTTTTTAAGGGAAGATAAAGAGGAAACATGAACATTAACCCCttt                                                             >  2:127087/1‑149 (MQ=255)
                                                                                         aGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCg                                                          <  1:631392/151‑1 (MQ=255)
                                                                                         aGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGGTAAAGAAATAAGAATTTCTTCTTTTTGGGGGGTCCCCGTCTCACTATTATCTCACCCCGGATTCTTTAGGGAAAGAGAAAGAGGGAACATAAACCATAACCCCttt                                                             >  2:785426/1‑148 (MQ=255)
                                                                                         aGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGCTTAAGACAAACTATTTTCTCCTTTTCGAACGGTTCCCTTTTCTCTTTTATTTCAAGCCGAATTTCTTAAGGGACGAAAAGGAGGGACCATGAACAATAACCCCGtt                                                              >  1:452495/1‑148 (MQ=255)
                                                                                           cATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAAAACCATTTTCTCCCTTTTTAGGGCTCCGCCTCTCCCTATTATTTTACCCCGTATTTTTTAAGGGAAGAAAAAGAGGGACCATAAAAAATAACCCCCTTTCGAc                                                        >  1:494763/1‑151 (MQ=255)
                                                                                                ccTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCttt                                                   <  2:60252/151‑1 (MQ=255)
                                                                                                 cTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAAAAAGATTTTCTCCCTTTCGAGTGGTCCCCCTCTCCCTATTTTTTCCCCCAGTATTTTTTAAGGGAAGATAAAGAGGGAACAATAACAATAACCCCCTTcccgcgcctt                                                   >  1:353353/1‑142 (MQ=255)
                                                                                                 cTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAACGATTTTCTCCTTTTATAGTGCTCCCCCTCTCCCTATTTTTTCCCCCCGTATTTTTTAAGGGAAGAAAAAGAGGGAACATGAAAATTAcccccctttccacgcctt                                                    >  1:946156/1‑136 (MQ=255)
                                                                                                   tAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTc                                                  <  2:258522/149‑1 (MQ=255)
                                                                                                    aaaaGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTCGc                                                <  1:69167/150‑1 (MQ=255)
                                                                                                     aaaGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAACGATTTTCTCCCTTTAAAGGGCTCCCCTTCTCCCTATTATTTCACGCCGGAATCTTTAAGGGAAGAAAAGGAGGAAACATGGACAATACCCCGGTTTccccgcctttccc                                                >  1:44042/1‑139 (MQ=255)
                                                                                                        gATATGTCAGGCTTGCGGAGTGGCGGTTAAGGAAAAACAATTTTCTCCCTTTTGAGTGCTCCCCCTCTCACTCTTATCTCCCCCCGCGTTTTTTAGGGAGAGCTAAAGGGGGACAAAGGACATTAACCCCCTTTCGCCGCTTTCCCCCGc                                            >  2:746789/1‑150 (MQ=255)
                                                                                                             gTCAGGCTTGCGGAGTGGCGGTTATGGAAATACAAATTTCTCCCTTTTGAGTGGTCCCCCTTTCCCTATTTTTTCCCCCAGGAATTTTTAAGGGAAGACAAAGGGGGAACAATGAAATTATCCCCTTTcccccgccttttgc                                               >  2:354509/1‑130 (MQ=255)
                                                                                                                  gCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATg                                 <  1:614913/151‑1 (MQ=255)
                                                                                                                      gCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCa                              <  1:729364/150‑1 (MQ=255)
                                                                                                                               cGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTAt                    >  1:525242/1‑151 (MQ=255)
                                                                                                                                      aGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGAc               >  1:530349/1‑150 (MQ=255)
                                                                                                                                      aGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACa              <  2:335097/151‑1 (MQ=255)
                                                                                                                                         aCATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAg             <  2:480941/149‑1 (MQ=255)
                                                                                                                                           aTACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtggc         >  2:550074/1‑151 (MQ=255)
                                                                                                                                             aCGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtggctg       >  1:705321/1‑151 (MQ=255)
                                                                                                                                             aCGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtggct        <  1:827535/150‑1 (MQ=255)
                                                                                                                                               gATTTCCTCCTTTCAGAGTGCTCCGCTTCTCCCTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtggctggc     <  2:249469/151‑1 (MQ=255)
                                                                                                                                                aTTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGgcg    >  2:374970/1‑151 (MQ=255)
                                                                                                                                                   tcctccTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCa  <  1:370587/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                          
TTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTATTCTTAAGGGAACGATAAGGAGGAACCATGAACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCA  >  NC_000913/1487016‑1487311

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: