Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,323,313 G→C 5.7% intergenic (‑32/+23) phnE ← / ← phnD phosphonate ABC transporter membrane subunit/phosphonate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,323,3130GC 5.7% 174.4 / 2.6 60intergenic (‑32/+23)phnE/phnDphosphonate ABC transporter membrane subunit/phosphonate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base C (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

GGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTC  >  NC_000913/4323171‑4323454
                                                                                                                                              |                                                                                                                                             
ggCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACGTTTTCCGATGCGACGCCgg                                                                                                                                               >  1:793349/1‑143 (MQ=255)
   ctaggtgaGAAGCGGGAACCAGGTGCGCTCGGGTGCGGCGACGGTGGTGGTTTGGATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCGTGCGCCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGc                                                                                                                                            <  2:99603/136‑1 (MQ=255)
       cTCAGAAGCGAGAACCAGCTGCGCTGGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCt                                                                                                                                        <  1:527046/143‑1 (MQ=255)
         cggAAGCTAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGGTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATAAGACGCCCCGGCGTcgca                                                                                                                                      <  2:509626/141‑4 (MQ=255)
          agagGCGGGACCGAGCTGGGCTTGGGTGGTGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCCCGGCGTcgcat                                                                                                                                     <  2:235749/143‑5 (MQ=255)
             ggCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCg                                                                                                                                  <  2:463787/142‑1 (MQ=255)
                 atacgcAGCTGCGCTTGGGTGGCCGATGGTCATGGTTTGCATGTATAGCAGATACCGGGACGCCGGTGCCGGGGCTTGCTCGGCATCCGGCATCAACGCCTGCACCTTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTa                                                                                                                             <  2:870817/138‑1 (MQ=255)
                 agaACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCAAcaa                                                                                                                              >  1:766425/1‑141 (MQ=255)
                       aGCTGCGCTTGGGTGGGGCGGTGGTGATGGTTTGGATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGAtt                                                                                                                        <  1:106850/143‑1 (MQ=255)
                        gCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGAttt                                                                                                                       >  2:515118/1‑143 (MQ=255)
                         cTGCGCTGGGGTGGGGCGATGGTGATGGTTTGCATGTGTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTa                                                                                                                      <  2:743303/143‑1 (MQ=255)
                             gcTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCTGCATCAACGCCTGCACATTGTCGGATGCGACGCCGCGGCGACGTATCCGACCTACGATTAACttc                                                                                                                  >  2:696641/1‑141 (MQ=255)
                             gcTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGc                                                                                                                  >  2:26900/1‑143 (MQ=255)
                                     ggggCCATGGTGATGGCTTGCATGTTTGGCTCCTGTCTGTAGTCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCttt                                                                                                          <  1:287475/143‑1 (MQ=255)
                                       ggCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTAc                                                                                                        >  1:610000/1‑143 (MQ=255)
                                        gCGATGCTGCTGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACt                                                                                                       <  1:817627/143‑1 (MQ=255)
                                          gatggtgatggtTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCa                                                                                                     <  1:667559/143‑1 (MQ=255)
                                          gatggtgatggtTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCAGGGGCGCCTTATCCGACCTACAATTTACTGCACCGCTTTACTCa                                                                                                     >  2:352010/1‑143 (MQ=255)
                                           gtggtgatggtTTGCATGTTGGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCAc                                                                                                    <  2:28553/142‑1 (MQ=255)
                                           ctggtgatggtTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCAc                                                                                                    <  1:270166/142‑1 (MQ=255)
                                               tgatggtTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCgca                                                                                                >  2:252676/1‑141 (MQ=255)
                                                     tttGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATc                                                                                          >  1:29145/1‑143 (MQ=255)
                                                     tttGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACCCATc                                                                                          >  2:506608/1‑143 (MQ=255)
                                                      gtGCATGTTTGGCTCCGGCCTGTAGGCAGGATAAGGTCCTTGCACCGCATCCGGCATCAACGCCTGCACACTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATcg                                                                                         <  1:239666/142‑1 (MQ=255)
                                                                  cTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCgt                                                                             >  1:872458/1‑143 (MQ=255)
                                                                   tCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCTGGTGTGGGCGCCTCCACATTGTCGGATGCGCCGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCACCGCGCTTAGCGCgtt                                                                            >  2:191177/1‑143 (MQ=255)
                                                                   tCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCgtt                                                                            >  1:429846/1‑143 (MQ=255)
                                                                   tCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCgtt                                                                            >  1:468673/1‑143 (MQ=255)
                                                                   tCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCgtt                                                                            >  2:131836/1‑143 (MQ=255)
                                                                     cGGGCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTAAGCGCgatga                                                                          <  2:408748/143‑2 (MQ=255)
                                                                       gTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTc                                                                        <  1:786964/143‑1 (MQ=255)
                                                                       gTCTGTAGGCCGGATAAGGTGCTGGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTc                                                                        >  2:110844/1‑143 (MQ=255)
                                                                            tAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGCTGCGACGCCGGGGCGTCTTATCCGACCTACGAGTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAACGCGTTGTACAGGCg                                                                   >  2:93432/1‑143 (MQ=255)
                                                                            tAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGAGCGTTGTTCAgcct                                                                   >  1:789739/1‑140 (MQ=255)
                                                                                ccGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCGCATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGAGTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTcc                                                               <  2:513276/143‑1 (MQ=255)
                                                                                   gATAGGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAgg                                                            <  2:796993/143‑1 (MQ=255)
                                                                                   gATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAgg                                                            <  2:23213/143‑1 (MQ=255)
                                                                                   gATAAGGCGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGAACTACGAGTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAgg                                                            <  2:12866/143‑1 (MQ=255)
                                                                                    aTAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGt                                                           <  1:848486/143‑1 (MQ=255)
                                                                                      aaGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCa                                                         <  1:617301/143‑1 (MQ=255)
                                                                                         gTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTCCAGGCGGTCCCGGTccacc                                                      >  2:289781/1‑139 (MQ=255)
                                                                                          tGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTCTCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGGCATCCa                                                     >  1:839844/1‑143 (MQ=255)
                                                                                           gCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAg                                                    >  2:214189/1‑143 (MQ=255)
                                                                                                  ccGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGcc                                             >  2:364799/1‑143 (MQ=255)
                                                                                                   cGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTCTCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCCGGTCATCCCGCTGCGccc                                            >  2:622802/1‑142 (MQ=255)
                                                                                                   cGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCt                                            >  2:316647/1‑143 (MQ=255)
                                                                                                         cGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCCGGCGGTCCAGGTCATCCAGCTGCGCCTGAAtct                                      >  1:855738/1‑143 (MQ=255)
                                                                                                         cGGCATCAACGCCTGCACATTGACGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAAtct                                      <  2:327256/143‑1 (MQ=255)
                                                                                                           gCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCg                                    <  1:874015/143‑1 (MQ=255)
                                                                                                            cATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGCCATCCAGCTGCGCCTGAATCTCgg                                   <  1:811437/143‑1 (MQ=255)
                                                                                                             aTCCACGCCTGCACATGGCCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCACGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCggt                                  <  1:594817/143‑1 (MQ=255)
                                                                                                                 aCGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGtt                              >  1:229751/1‑143 (MQ=255)
                                                                                                                   gCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGtttt                            <  1:678677/143‑1 (MQ=255)
                                                                                                                         acaTTGTCGGATGCGACGCCGGGGCGTCTTATCAGACCTAAGATCTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAg                      <  1:743279/143‑1 (MQ=255)
                                                                                                                              gTCGGATAAGACGCCCCGGCGTCGCATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTGTTGCCAGCTTGt                 >  2:520055/1‑143 (MQ=255)
                                                                                                                                   aTGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGGACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGc            <  1:79120/143‑1 (MQ=255)
                                                                                                                                   aTGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACCCACCGAACTCATCGCGCTTATCTCGTTGTTCAGTCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGc            <  1:79270/143‑1 (MQ=255)
                                                                                                                                       gACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGACTCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCAt        >  1:696267/1‑143 (MQ=255)
                                                                                                                                       gACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACCCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTCGTTGCCAGCTTATCCTGCTCAt        >  2:494224/1‑143 (MQ=255)
                                                                                                                                        aCGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCAtt       <  1:391887/143‑1 (MQ=255)
                                                                                                                                          gCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTCCGCCTCGATCTCGGTGGGGTTTGCCAGCTTGTCCTGCTCATTCa     <  1:248286/143‑1 (MQ=255)
                                                                                                                                             ggggCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTc  <  1:214189/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTC  >  NC_000913/4323171‑4323454

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