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.94.5% 224.0 / 4.0 55pseudogene (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/3);  new base . (24/28);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
*NC_0009132,173,3620C.96.3% 222.1 / 1.3 55pseudogene (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/2);  new base . (24/28);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACC  >  NC_000913/2173219‑2173493
                                                                                                                                              ||                                                                                                                                   
aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCACCCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTCAAATGGTGAGTGGGATACAAATCAGGCTCGCCGATaccac                                                                                                                                      >  2:1154376/1‑143 (MQ=255)
 gATTTCCACGATGCACGGCGAGCGCCAGCGCCACGAAGAGGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccacc                                                                                                                                     <  1:177010/143‑1 (MQ=255)
  aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                                                                                                    >  2:739726/1‑142 (MQ=255)
   tttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCCGCACCTGCTTACCCGGCACACAAACAGCAATTCAAAAGCCGAGTGGGATACAAACCAGGCTTGCCAATaccaccgc                                                                                                                                   >  2:53845/1‑141 (MQ=255)
       cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACACCAATCAAAATGGCGACTCAGATAAAAATCAGGCTTGCCGATaccaccgccgta                                                                                                                               >  2:445594/1‑137 (MQ=255)
              cACGGCGACCGCCATCGCCACGAAGAGGGCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAGCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                      <  2:721773/143‑1 (MQ=255)
                 ggcggcCGGCATGGCGACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCAGCTGATTACCCGGCTCACAAGCAGCAATTAAAGTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                   <  2:626662/143‑1 (MQ=255)
                    gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAgcagca                                                                                                                >  2:314406/1‑143 (MQ=255)
                     aCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAgcagcag                                                                                                               >  1:614203/1‑142 (MQ=255)
                      ccGCCATCGCCACGAACAAGCCGATGGTGGCAATATCGCCAAACGGCAGCACCTGACTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                                                                                                              <  1:204988/143‑1 (MQ=255)
                       cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCGGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCGGCAGCAACg                                                                                                             >  2:997388/1‑143 (MQ=255)
                       cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                                                                                             >  1:730765/1‑143 (MQ=255)
                         ccATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                                                                                                           >  2:404811/1‑143 (MQ=255)
                          cATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCg                                                                                                          <  2:311987/143‑1 (MQ=255)
                                    aagaagCCGAGGGGGGCAAGATCACCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                <  2:927994/143‑1 (MQ=255)
                                       tggCCGTTGGGGGCAAGGCGGCCAGACGGCAACACCTGACTGCCCGGCACACAAACAGCAGTTGAAATGGTGGGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAGCGCCGGATCAAGGCCAAt                                                                                             <  2:493200/141‑1 (MQ=255)
                                       aggCCGATGGTGGCCAGGTCGCCAAACGGCAGCACCTGATTACCCGGGACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCTAAt                                                                                             <  2:208555/141‑1 (MQ=255)
                                           cgttggtggCAAGTTCGCCAAACGGCAGCACCTGATTGCCCGTCACACAAACAGCAATTAAAATGGTGAGTTGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                         <  2:454917/140‑1 (MQ=255)
                                             aTGGTGGCAAGATCGCCAAACGGCACCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTCGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                       <  2:950912/143‑1 (MQ=255)
                                                 tggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc                                                                                   <  2:1102587/143‑1 (MQ=255)
                                                   gCAAGATCGCCAAACGGCAGCACCTGCTTACCCAGAACACAAACAGCAATCAAAATGCTGAGTGGGATAAAAATCAAGCTTGCCGATACCA‑‑CCGCCGTATCTCCCAACCGCAACTCCCAATCAAACCCACTCACCAACTCCtc                                                                                 >  2:826008/1‑142 (MQ=255)
                                                   gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGGCTGTCCCAGCAGCAACGCCGCATCAAGGCCAACCAGGCACTCCTg                                                                                 >  2:77412/1‑143 (MQ=255)
                                                   gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                 >  1:54784/1‑143 (MQ=255)
                                                    cAAGATCGCCACACGCCAGCACCTGATTACCCGGCACACAACCAGCAATTAAAATGGTGAGTGTGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGa                                                                                <  2:416740/143‑1 (MQ=255)
                                                        aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                            >  1:557825/1‑143 (MQ=255)
                                                        aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATACAGGCCCACTCGGAACTCCTGAccgc                                                                            >  2:1055343/1‑143 (MQ=255)
                                                         tCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                           <  2:1034183/143‑1 (MQ=255)
                                                           gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGa                                                                         <  1:1129649/143‑1 (MQ=255)
                                                            caaaaCGTCAGCACGTCATTACCGGGCACACAAACGGCAATGACAATGGCGAGTGGGATAAAAATCCGGCTTGCCGATACCA‑‑CCGCCGGATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                        <  2:941999/141‑1 (MQ=255)
                                                             cagaCGGCAGCTCCTGATTCCCCGGCACTCGAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGTCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                       <  2:368155/140‑1 (MQ=255)
                                                             cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCGAGCAGGAACTCCTCAccgcccgac                                                                       >  2:836235/1‑139 (MQ=255)
                                                             cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGGCGGATCAAGGCCAATCAGGAACTCCTGACCGGCGAAc                                                                       >  1:403400/1‑143 (MQ=255)
                                                                        aCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACCCCTGACCGCCGAACTTAGCCTGTag                                                            >  2:546010/1‑142 (MQ=255)
                                                                         gcgtactACCCGTAACACAAACTGCTATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAGCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaa                                                           <  2:653931/137‑1 (MQ=255)
                                                                            gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                        <  1:404811/143‑1 (MQ=255)
                                                                            gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                        >  1:5121/1‑143 (MQ=255)
                                                                            gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                        <  1:1061456/143‑1 (MQ=255)
                                                                               tACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGAAAACGACTa                                                     >  2:976228/1‑143 (MQ=255)
                                                                                aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAc                                                    >  1:694349/1‑143 (MQ=255)
                                                                                           ccgatctaaTTAAAATGGTGAGTGGGATAACAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGGCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACTGCCGAACTTCGCCTGTAAACGACTACGAGCCGGCTTa                                         <  1:517628/136‑1 (MQ=255)
                                                                                           aaaCAGCAATTAAAATAGTGAGGGGGGTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                         <  1:1085757/143‑1 (MQ=255)
                                                                                              cAGCAATTACCATGGTGAGTGGGATAAAAATCAGGCTTGCCGACACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                      <  1:88734/143‑1 (MQ=255)
                                                                                              cAGCAATTAAACTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                      <  1:53845/143‑1 (MQ=255)
                                                                                                  aaTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                  >  2:1142581/1‑143 (MQ=255)
                                                                                                   aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCAGAAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAAGACTACGAGCCTGCTTAGAGATcgg                                 >  2:517628/1‑143 (MQ=255)
                                                                                                    ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CAGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACACCTACGAGCCAGCCCAGCCATgcac                                >  2:758322/1‑140 (MQ=255)
                                                                                                         aTGGTGAGTGGGATAACAATCAGGCTGGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCCGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                           <  2:1106670/143‑1 (MQ=255)
                                                                                                          tGGTGAGTGGGAGAAAACTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACCGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                          <  2:159006/143‑1 (MQ=255)
                                                                                                            gTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                        <  1:739726/143‑1 (MQ=255)
                                                                                                             tGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTCAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcca                       >  1:245789/1‑143 (MQ=255)
                                                                                                               aGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGAAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                     >  2:164375/1‑143 (MQ=255)
                                                                                                                  gcgATAAACACCAGGCTTGCCGTTAC‑‑CCGCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatccat                  <  2:849862/141‑1 (MQ=255)
                                                                                                                     aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCCGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat               >  2:818046/1‑143 (MQ=255)
                                                                                                                       aaaacTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCTGCTTCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg             <  1:314406/143‑1 (MQ=255)
                                                                                                                        aaaaTCAGCCTTGCCGATATCACCC‑‑CCGTATGCCCCAGCAGCAACGCCGGCTGAACGCAAATCAGGTACTCCTGCCCGCCGAACTGCGCCTGTAAACGACTACGATCCTGCTTATCGATGGGCGTTAAACCATCCATGATggg            <  2:295274/143‑1 (MQ=255)
                                                                                                                          aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCGGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt          <  1:144411/143‑1 (MQ=255)
                                                                                                                          aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGACATAGCACCAACGCCGTATCAAGGCCAATCAGGAACACCTGACCGCGGAACTTAGCCCGTAAACGACTACGAGCCTGCCTAGCGATGGGCGTTAAACCATCCATGATGGGtt          >  2:139430/1‑143 (MQ=255)
                                                                                                                          aaTCAGGCTGGCCGATACCA‑‑CCGCCGTATGCCCCAGCAGCAACGCCGGAGCAAGGCCAATCAGGAACTCCTGACCGCCCAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt          <  2:578807/143‑1 (MQ=255)
                                                                                                                               ggCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATc     >  1:516084/1‑143 (MQ=255)
                                                                                                                                  gtGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCCACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc  <  2:342355/142‑1 (MQ=255)
                                                                                                                                              ||                                                                                                                                   
AGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACC  >  NC_000913/2173219‑2173493

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: