Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,173,361 Δ2 bp 100% pseudogene (917‑918/1358 nt) gatC ← galactitol‑specific PTS enzyme IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3610C.89.7% 212.4 / 14.5 58pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/6);  new base . (24/28);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 3.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.
*NC_0009132,173,3620C.89.7% 210.2 / 9.4 58pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/3);  major base . (24/28);  minor base G (0/3);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173219‑2173496
                                                                                                                                              ||                                                                                                                                      
aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccac                                                                                                                                         >  1:899444/1‑143 (MQ=255)
aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccac                                                                                                                                         >  1:387876/1‑143 (MQ=255)
  aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                                                                                                       >  2:488177/1‑142 (MQ=255)
   tttCGACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgc                                                                                                                                      <  2:997364/143‑3 (MQ=255)
          gATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                             <  1:743401/143‑1 (MQ=255)
               aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGCCCcag                                                                                                                        <  1:807940/143‑1 (MQ=255)
                 ggCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                      <  1:294093/143‑1 (MQ=255)
                  gCGTCCGCCATCGCCCCGAAGAAGCCGACGGTGGCACGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGCGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagcag                                                                                                                     <  2:952935/143‑1 (MQ=255)
                   cGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTCTGTCCCagcagc                                                                                                                    >  1:80886/1‑143 (MQ=255)
                    gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAgcagca                                                                                                                   <  1:1024231/143‑1 (MQ=255)
                      ccGTCATCGCCGCGTCGAATCCGGTGGTGGCAAGATCGCCAAACTGCAGCACCGGATTACCCTGCACACAAACAGCAATTAAAATGGTGAGTGTGATAAATAGCTCGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                                                                                                                 <  2:874621/143‑1 (MQ=255)
                           aTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCgg                                                                                                            <  2:147681/143‑1 (MQ=255)
                            tCGCCACGGAGAAGGCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                           <  1:748027/143‑1 (MQ=255)
                                cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                                       <  2:202504/143‑1 (MQ=255)
                                   gcacgagCCGGTCGTGGCAAGAGCGCCAAGCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAGCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                    <  2:146686/138‑1 (MQ=255)
                                       atgCCGATGGAGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAGACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                                <  1:450628/141‑1 (MQ=255)
                                         gCCGTTGGTGGCCAGATTCCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGTATCAAGGCCAATCa                                                                                              <  2:866879/143‑1 (MQ=255)
                                          ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                             <  2:53717/143‑1 (MQ=255)
                                                    cAAGATCGCCAAACGGCAGCACCTAAATACCCACCCACCACAAAGTGGTTAAGATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGa                                                                                   <  1:488177/143‑1 (MQ=255)
                                                        aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                               <  1:115182/143‑1 (MQ=255)
                                                        aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                               >  2:956551/1‑143 (MQ=255)
                                                           gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGGATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGa                                                                            >  1:437712/1‑143 (MQ=255)
                                                            ccAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCCGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                           <  1:547604/143‑1 (MQ=255)
                                                              aaaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGTCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCACGGCCAACCAGGAACTCCTGGCCGCCGCAcc                                                                         >  2:793762/1‑142 (MQ=255)
                                                                     agcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTg                                                                  <  1:1074733/143‑1 (MQ=255)
                                                                            gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTACACGa                                                           >  1:510336/1‑143 (MQ=255)
                                                                            gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAGCAAGGCCAATCAGGCACCCCTGACCGCCGAACTCCGCCCGCAAACGa                                                           >  2:780225/1‑143 (MQ=255)
                                                                               tccccGGCACCCAAGCAGCAATGAACATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                        <  2:357579/141‑1 (MQ=255)
                                                                                 cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACCACg                                                      >  1:799455/1‑143 (MQ=255)
                                                                                  ccGGCACGCAAACAACAATGAAAATGGTGAGTGGTATAAACATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGa                                                     <  2:630993/143‑1 (MQ=255)
                                                                                  ccGGCACACAAACAGCAATTAAAATGGTGAGTTGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCGCAGCAGCAACGCCGGATCAAGGCAAATCAGGAACTCCTGACCGCCGAACTGCGCCCGTAAACGACTACGa                                                     >  2:131281/1‑143 (MQ=255)
                                                                                  ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCACAAACGGTCAGGCCAGGGCAGTCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGa                                                     >  1:1065520/1‑143 (MQ=255)
                                                                                  ccGGCACAAAAAAAACAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGa                                                     <  2:61507/143‑1 (MQ=255)
                                                                                   cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                    >  1:504566/1‑143 (MQ=255)
                                                                                    ggCACACCAACAGCAATTAAAATGGTGAGCGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                   <  1:137432/143‑1 (MQ=255)
                                                                                        cacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACATCGCCTGTAAACGACAACAAGCCTGc                                               >  2:447020/1‑143 (MQ=255)
                                                                                           aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCAGCTTa                                            >  1:462014/1‑143 (MQ=255)
                                                                                            agcAGCACTTAAAAGGTGGCGTGGGATAACCATGCGCCTTTCCTATAC‑‑CCCCGCCGTATGTCCCAGCAGCGACCGCGGATCAAGGGCAATCAGGAACTCCTGGCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAg                                           <  2:64729/141‑1 (MQ=255)
                                                                                               aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCACGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACAACTACGAGCCTGCTTACCGa                                        >  1:423441/1‑143 (MQ=255)
                                                                                               aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGATCGCCGAACTTCGCCTGTAAACGACTACGAGACTGCTTcgcgt                                        >  1:659817/1‑142 (MQ=255)
                                                                                               aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAACCGACTACGAGCCTGCTTaccca                                        >  1:176727/1‑139 (MQ=255)
                                                                                                  aaggaaaaTGGGGAGTGGGAGAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                     <  1:982478/139‑1 (MQ=255)
                                                                                                  aaTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                     <  2:387876/143‑1 (MQ=255)
                                                                                                    ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACCGTCTACGAGCCTGCTTAGCGATggac                                   >  1:808199/1‑141 (MQ=255)
                                                                                                         aTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGGTaa                              >  1:618377/1‑143 (MQ=255)
                                                                                                          cggggAGTGGGATAAAAGTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                             <  1:920874/139‑1 (MQ=255)
                                                                                                              gAGTGGGATAAAAACCAGGCTGGCCGATACCA‑‑CCGCCGTATGTCCCAGTGGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                         <  2:500012/143‑1 (MQ=255)
                                                                                                               aGTGGGATAAAAATCAGGCTTGCGGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                        <  1:312622/143‑1 (MQ=255)
                                                                                                                 tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCAATGGGAGTTAAAccatcca                      >  2:18294/1‑143 (MQ=255)
                                                                                                                    taTAAAAATCAGGCTTGCCGATAC‑‑CCCCTCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                   <  2:168754/142‑1 (MQ=255)
                                                                                                                       gaaaaTCCGGCTTGCCGATAC‑‑CGCCGGCGTATGTCCCAGCGGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTCCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                <  2:756667/142‑1 (MQ=255)
                                                                                                                         aaaCCAGGCTTGCCGATACCA‑‑GCGCCGTATGTCCCCGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGGCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt              <  1:622481/143‑1 (MQ=255)
                                                                                                                          gaTCAGGCGTGCCGATAC‑‑CGCCGCCCTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt             <  2:560209/142‑1 (MQ=255)
                                                                                                                          aacccggCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGATCTCCTGACCGCTGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt             <  1:216733/138‑1 (MQ=255)
                                                                                                                          aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATAAGGAACTCCTAACCGCCGGACTTCGCCTGTAAACGACTACGAGCCAGCTTAGCGATGGGCGTTAAACCATCCatgatgcgat             >  1:149037/1‑139 (MQ=255)
                                                                                                                           aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt            <  1:958977/143‑1 (MQ=255)
                                                                                                                             gaGGCATGGCGATACCA‑‑CCGCTGTATGTCCCCGCAGCAACGCCGGATCAAGGCTAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTATCGATGGGCGTTAAACCATCCATGATGGGTTTaa          <  2:58127/142‑1 (MQ=255)
                                                                                                                                gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCa       >  1:631402/1‑143 (MQ=255)
                                                                                                                                gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGCTTAATCa       >  1:316008/1‑143 (MQ=255)
                                                                                                                                gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGCGCGTTAAACCATCCATGATGGGGTTAATCa       >  1:168049/1‑143 (MQ=255)
                                                                                                                                 cTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTCAACGATCTATGATGGGGTTAATCAc      >  2:166855/1‑143 (MQ=255)
                                                                                                                                     ccGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGCGCGTTCAACCACCCCTGATGGCTTTACTCCCGCGt  >  2:109267/1‑143 (MQ=255)
                                                                                                                                              ||                                                                                                                                      
AGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173219‑2173496

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: