Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑1/+1) gatC ← / ← gatC pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC/pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3610C.89.2% 265.2 / 15.3 65pseudogene (2/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/7);  new base . (36/22);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 2.24e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3620C.90.6% 260.2 / 11.0 65pseudogene (1/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/6);  new base . (36/22);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 5.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCA  >  NC_000913/2173212‑2173509
                                                                                                                                                      ||                                                                                                                                                   
gCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATaccacc                                                                                                                                                     >  2:1021285/1‑151 (MQ=255)
gCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATaccac                                                                                                                                                      <  2:779779/150‑1 (MQ=255)
       aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATaccaccgccgtat                                                                                                                                              <  1:520466/151‑8 (MQ=255)
          tttCAACGATGCCCGGGGACCGCCATCGCCACGAAGAAGCAGGTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACATCAATTAAATTGGTGTGTGGTTTAAAAGCCT‑GGCTTGCCGATTCCA‑‑CCGCCGTATGt                                                                                                                                          <  2:602498/150‑1 (MQ=255)
          ttcccccGATGCCCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATg                                                                                                                                           <  2:1143610/143‑1 (MQ=255)
             ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTccc                                                                                                                                       >  2:1143274/1‑150 (MQ=255)
             ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                                      <  1:1160153/151‑1 (MQ=255)
                  atccACGGCGACGGCCAGCGCCACGAAGAAGCCGATGGTGGCAAGAGCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagcag                                                                                                                                  <  2:759099/147‑1 (MQ=255)
                  aTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                                   <  2:233364/149‑1 (MQ=255)
                           gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCgg                                                                                                                         >  1:36317/1‑150 (MQ=255)
                             ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCgg                                                                                                                         <  2:326730/148‑1 (MQ=255)
                              cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                      >  1:275091/1‑150 (MQ=255)
                              cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                      >  1:429704/1‑150 (MQ=255)
                               gCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCa                                                                                                                     >  2:521144/1‑150 (MQ=255)
                               gCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCa                                                                                                                     >  1:1172607/1‑150 (MQ=255)
                                   tCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                >  1:770487/1‑151 (MQ=255)
                                       cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                            >  2:1117938/1‑151 (MQ=255)
                                         cGAAGAAGCCGATGGGGGCAAGATCGCCAAACGGGAGCACCGGATTCCCCGGCACACAAACAGAAATTAAAATGGTGAGTGGGATAAAAATCA‑GACTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                                          <  2:382783/151‑1 (MQ=255)
                                          gaagaaGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                                         >  1:1066030/1‑151 (MQ=255)
                                          gaagaaGCCGATGGTGGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                                          <  2:24978/150‑1 (MQ=255)
                                           aagaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                                          <  2:89186/149‑1 (MQ=255)
                                                 ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                                  <  2:103596/151‑1 (MQ=255)
                                                         ggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                          <  1:376251/151‑1 (MQ=255)
                                                             aGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                                      >  2:624137/1‑151 (MQ=255)
                                                               aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                                     <  1:704483/150‑1 (MQ=255)
                                                                 cGCCAAACGGCAGCCCCTGATTCCCCGGCACACAAACAGAAATAAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                                    <  2:722982/149‑1 (MQ=255)
                                                                 cGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                                  >  1:84402/1‑151 (MQ=255)
                                                                   ccAACCGGCAGCACCTGATTACCCGGCCCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAACCA‑GGCTTCCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                                <  1:599937/151‑1 (MQ=255)
                                                                    cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                                >  1:129405/1‑150 (MQ=255)
                                                                       aCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTa                                                                             >  2:1085026/1‑150 (MQ=255)
                                                                           cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACg                                                                         >  2:49531/1‑150 (MQ=255)
                                                                            agcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGAc                                                                       >  2:848767/1‑151 (MQ=255)
                                                                                ccTGATTACCCGCCACAAAAACAGCAATAAAAAGGGTGAGTGGGAAAAAAACCA‑GGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                                     <  2:669851/149‑1 (MQ=255)
                                                                                  tGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                   >  2:1060049/1‑149 (MQ=255)
                                                                                   gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                >  2:1095114/1‑151 (MQ=255)
                                                                                    aTTACCCGGCACACAAAAACCAATTAAAAGGGTGAGTGGGATAAAAATCA‑GGTTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                                 <  2:506404/149‑1 (MQ=255)
                                                                                     ttccccGGCACACAACCACCAATTAAAATGGAGAGGGGGAAAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                                 <  2:39269/145‑1 (MQ=255)
                                                                                     ttACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCt                                                              >  2:463754/1‑151 (MQ=255)
                                                                                        cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                                           >  2:1144094/1‑151 (MQ=255)
                                                                                        cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                                           >  1:1127236/1‑151 (MQ=255)
                                                                                        cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                                           >  2:747557/1‑151 (MQ=255)
                                                                                           ggCACACAAACAGAAATTAAAATGGGGAGGGGAATAAAAACCA‑GGTTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                                         <  1:775619/150‑1 (MQ=255)
                                                                                             cacacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                      >  1:782263/1‑151 (MQ=255)
                                                                                              acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                      >  2:176839/1‑150 (MQ=255)
                                                                                                acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                    >  1:30300/1‑150 (MQ=255)
                                                                                                  aaaCAGCAATAAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATggg                                                 <  1:973094/151‑1 (MQ=255)
                                                                                                      aGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                             >  1:328224/1‑151 (MQ=255)
                                                                                                      aGCAATTAAAATGGTGAGTGGAATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                             <  1:717826/151‑1 (MQ=255)
                                                                                                        cAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                                            >  1:1188577/1‑150 (MQ=255)
                                                                                                         aaTTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                                           >  1:635988/1‑150 (MQ=255)
                                                                                                         aaTTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                                           >  2:231898/1‑150 (MQ=255)
                                                                                                         aaTTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                                           >  2:1042796/1‑150 (MQ=255)
                                                                                                         aaTAAAAATGGTGAGGGGAATAAAAATCG‑GGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                             <  2:1097452/148‑1 (MQ=255)
                                                                                                             aaaaagggtggggggAAAAAAAACAGGGCTGGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcca                                       <  2:680440/139‑1 (MQ=255)
                                                                                                             aaaaTGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑ACGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                                      >  2:140868/1‑151 (MQ=255)
                                                                                                             aaaaTGGGAAGGGGGATAAAAACCA‑GGCTTCCGGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcca                                       <  1:67879/150‑1 (MQ=255)
                                                                                                                  ggtggggggaaaaaaaaTCG‑GGTTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                                 <  2:822131/139‑1 (MQ=255)
                                                                                                                   gTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                                 >  2:99999/1‑150 (MQ=255)
                                                                                                                    tGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                               >  1:1142388/1‑151 (MQ=255)
                                                                                                                    tGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                               >  1:1087601/1‑151 (MQ=255)
                                                                                                                     gAGTGGGATAAAAATCA‑GGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                               <  1:1075811/150‑1 (MQ=255)
                                                                                                                         gggATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt                          >  2:150373/1‑151 (MQ=255)
                                                                                                                            aaaaaaaCCG‑GGCTTGCCAATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa                       <  2:1081472/149‑1 (MQ=255)
                                                                                                                              aaaaaTCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                      >  2:1188297/1‑150 (MQ=255)
                                                                                                                               aaaaaccgcgtttGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATAAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                      <  1:671985/139‑1 (MQ=255)
                                                                                                                                aaaTCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAc                   <  1:325409/151‑1 (MQ=255)
                                                                                                                                 aaTCA‑GGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc                  <  1:140217/151‑1 (MQ=255)
                                                                                                                                  aTCA‑GGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                 >  1:1113054/1‑151 (MQ=255)
                                                                                                                                          gcttgccaaaccccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg              <  2:1185725/137‑1 (MQ=255)
                                                                                                                                                    ccccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAgcagca  <  2:1192865/149‑1 (MQ=255)
                                                                                                                                                      ccccACCGTATGTCCCAGCAGCAACGCGGGATCAAGGCCAATCAGGAACTCCTAACCGCCGAACTTCTCCTGTAGACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCa           <  2:485135/140‑1 (MQ=255)
                                                                                                                                                      ||                                                                                                                                                   
GCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCA  >  NC_000913/2173212‑2173509

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