Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 53,376 G→A 8.0% P14L (CCC→CTC)  pdxA ← 4‑hydroxythreonine‑4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091353,3760GA 8.0% 135.0 / 2.8 60P14L (CCC→CTC) pdxA4‑hydroxythreonine‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (32/23);  new base A (3/2);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

GGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATG  >  NC_000913/53236‑53513
                                                                                                                                            |                                                                                                                                         
gggAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTGAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCGGGACAACTAAGACCGGGCGAATCCCGGCGGGc                                                                                                                                         >  2:32746/1‑143 (MQ=255)
  gAATATGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCGGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTc                                                                                                                                       <  1:823887/143‑1 (MQ=255)
  gAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTc                                                                                                                                       >  2:858536/1‑143 (MQ=255)
    aTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGGGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGc                                                                                                                                     <  1:190413/143‑1 (MQ=255)
    aTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCACAGGCACAAACAAGCGGTTCGACCGGCCCCTCAACCCGGCCAAGCCGGACAACTAAGTCCGGGCCAGTCCCGGCGGGCTCGc                                                                                                                                     >  2:178647/1‑143 (MQ=255)
       aGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCgg                                                                                                                                  <  2:728280/143‑1 (MQ=255)
         ggcggcGGGTGAGCGGGCACCCGCGCCTCGCTCGCCGGTTGCGAAGGACAGTGGCATCGGCACAAACACCCAGTTCGACCGTCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCGGGa                                                                                                                                <  2:545492/143‑1 (MQ=255)
          gCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACGCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCCATCCCGGCGGGCTCGCCGGGAg                                                                                                                               >  2:400274/1‑143 (MQ=255)
                  tGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGCGTGATccca                                                                                                                       >  2:958346/1‑143 (MQ=255)
                          caaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGCGTGATCACACCACGCTg                                                                                                               >  2:434135/2‑143 (MQ=255)
                           aaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGGGATCACAACACGTTgg                                                                                                              >  1:428431/1‑143 (MQ=255)
                           aaCCGGGCAGCGATGCTCGGCTGGTAAAGAGAGTGGCATCGGCACAAACAACCAGTTCGACCAGCCACTCACGCAGCGCAAGGTGGACAACCAAGTCCGGGCCAATCCCGGCCACCCCGCCGCCAGCGAGCACAACACGtgtg                                                                                                              >  2:757637/1‑140 (MQ=255)
                            aCCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTggg                                                                                                             <  2:347892/143‑1 (MQ=255)
                                           ccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGGGTGATCACAACCCGTTGGGTTTTAACCCTTAGtt                                                                                              >  1:320529/1‑143 (MQ=255)
                                                  gTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTGACCATTAGTTGCTCAgc                                                                                       >  2:8284/1‑142 (MQ=255)
                                                   tAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTCGTTGCTCAggg                                                                                      >  1:45908/1‑142 (MQ=255)
                                                       gagaGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGAtttt                                                                                  >  1:739058/1‑143 (MQ=255)
                                                       gagaGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGAAAACTAAGTACGGGCCAATCACGGCGGGCTCGCCGGGAGTGATCACAAGACGTTGGGTTTTAACCATTAGTTGCTCAGGAtttt                                                                                  >  1:800328/1‑143 (MQ=255)
                                                          agTGGCATCGGCACAAACCACCCGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAAc                                                                               <  2:167384/143‑1 (MQ=255)
                                                            tGGCATCGGCACAAACAACCAGTTCGCCCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGt                                                                             >  1:468811/1‑143 (MQ=255)
                                                               cATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTaaactcga                                                                          >  1:758394/1‑137 (MQ=255)
                                                                   ggCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGAAAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGAGATCACCACACGTTGGGTTTTAACCATGAGTTGATCAGGATTTTAACGCAGGCGCt                                                                      >  2:163495/1‑143 (MQ=255)
                                                                    gCACAAGCCACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTg                                                                     <  2:30855/143‑1 (MQ=255)
                                                                    gCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAAAGTAGGCGCTg                                                                     >  1:297939/1‑143 (MQ=255)
                                                                     cacaAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACGACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTgg                                                                    >  2:297434/1‑143 (MQ=255)
                                                                        aaacaaCCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGCTTTTAACCCTTAGTTGCTCAGGATTTTAACGTAGCCGCTGGCAc                                                                 >  1:721657/1‑143 (MQ=255)
                                                                        aaacaaCCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGACCGGGCCAATCCCGGCGGGCTCGCCGGGAATGATCACAACACGTTGGGTCTTAACCCACAGTTGCTCAGGATTTTAACCTAGGCGCTGGCAc                                                                 >  2:446482/1‑143 (MQ=255)
                                                                         aacaacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACg                                                                <  2:719043/143‑1 (MQ=255)
                                                                               cAGTTGGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTc                                                          <  1:432013/143‑1 (MQ=255)
                                                                                  ttCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTg                                                       <  1:349597/143‑1 (MQ=255)
                                                                                    cGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCAGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCa                                                     <  1:598434/143‑1 (MQ=255)
                                                                                      aCCGGCCAGGCACGCTGTGCGTGCGGGATGACTAGGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATc                                                   <  2:166652/143‑1 (MQ=255)
                                                                                        cGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCGATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGGTCTTCACGTAAGCAAAGGTACGTTGTTCCCGCGTCCa                                                 >  2:681915/1‑143 (MQ=255)
                                                                                        cGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCa                                                 <  1:858536/143‑1 (MQ=255)
                                                                                          gCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACCCGTTGGGTTTTAACCATTAGTTGCGCAGGCTGTTAACGTCGGCCCTGGCACGTTGTTCCTGCATCCAGc                                               >  2:338232/1‑143 (MQ=255)
                                                                                               tcccGCTGCGCAGGCTGGACAACTAATTCCGGACCAATCCCTGCGGGCTCGCCGGGCGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTtgcg                                          <  2:488220/140‑2 (MQ=255)
                                                                                                gaggCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCACGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAGGGTCAGCGCAAGCAGGTTGTTTTTGTTTTGAGCTtgctg                                         <  2:228814/140‑1 (MQ=255)
                                                                                                cACGCTGTGCAAGCTGGACAACTCGCTCCGGGCCAGTCCCGGCGGGCTCGCCGGGAGTGATCACCACGCGTTGGGTTTTAGCCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTtgctg                                         <  2:763255/143‑1 (MQ=255)
                                                                                                    cTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTGCCTGCATCCAGCTTGCGGcttc                                     >  2:883870/1‑143 (MQ=255)
                                                                                                        gCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTcc                                 <  1:223010/143‑1 (MQ=255)
                                                                                                             cTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGaa                            <  2:439997/143‑1 (MQ=255)
                                                                                                              tGGACAACTAAGTCCGGGCCAATCCCGGCGGGATCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGACTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAAc                           >  2:874356/1‑143 (MQ=255)
                                                                                                                gACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACtt                         >  1:180326/1‑143 (MQ=255)
                                                                                                                gACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGGTTTAACCAGCAGTTGCGCAGGATGTTAACGGAGGCGCTGGCACGTGGTGCCTGCCTGCCTCTTGCTGCCTCTTCCCAGCACtt                         >  2:414550/1‑143 (MQ=255)
                                                                                                                 aCAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTa                        <  1:325329/143‑1 (MQ=255)
                                                                                                                    aCTAAGTCGGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACgg                     <  2:794113/143‑1 (MQ=255)
                                                                                                                    aCTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACgg                     <  1:1007652/143‑1 (MQ=255)
                                                                                                                     cTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGGACTCACGGt                    >  1:798172/1‑143 (MQ=255)
                                                                                                                     cTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGGTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCCGCATCCCGCTTGCTGCTCCTTCCGGGAACTTCCGGt                    >  2:336418/1‑143 (MQ=255)
                                                                                                                         gTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTtcat                <  1:400274/143‑1 (MQ=255)
                                                                                                                           ccGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTtcatca              >  1:1241/1‑143 (MQ=255)
                                                                                                                                tcactCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATg         >  1:997739/5‑143 (MQ=255)
                                                                                                                                  tttCCTGGCGGGCTCGCGGGGAGTTATTACCACACGTTGGGTTTTATCCATTAGTTGCTCAGGATTTTAACGTGGGCGCTGGCACGTGGTTTCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCg       <  2:189175/141‑1 (MQ=255)
                                                                                                                                  attCCCGGCTGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTATCCATTAGTTGCTCAGTATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCCGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCg       <  2:355570/141‑1 (MQ=255)
                                                                                                                                  actCCCGGCGGGCTCGCCGGGAGTGCTCACAACCCGTTGGGTTTTAACCATTAGTTGCCCAGGCTTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTCCCGCGACCTTACGCTCCATCAGCATGCg       >  2:275266/3‑143 (MQ=255)
                                                                                                                                  aaTCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCGGAGAACTTACGGTTCATCAGCATGCg       >  1:301692/1‑143 (MQ=255)
                                                                                                                                   aTCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCGCAGGATTTTAACCTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCgg      >  1:344581/1‑143 (MQ=255)
                                                                                                                                    tCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTGACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGGTTGCTGTGTCGTCCGAGAGCTCACGGTTCGTCAGCATGCgct     >  1:193153/1‑141 (MQ=255)
                                                                                                                                    tCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCACCATCATGCGGt     >  1:367226/1‑143 (MQ=255)
                                                                                                                                     cccGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCTAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTa    <  1:202315/143‑1 (MQ=255)
                                                                                                                                      ctgtcggGCTCGGCGGGAGTGATCACAACACGTTGGGTGTTAACCATTAGTTGCTCGGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGTTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTAt   <  2:908120/139‑1 (MQ=255)
                                                                                                                                       cggcgaGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATg  >  1:427136/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                         
GGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATG  >  NC_000913/53236‑53513

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: ‑