Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 875,531 C→T 10.2% T66M (ACG→ATG)  yliF → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913875,5310CT10.2% 204.7 / 8.9 67T66M (ACG→ATG) yliFputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (35/25);  new base T (3/4);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTG  >  NC_000913/875383‑875678
                                                                                                                                                    |                                                                                                                                                   
aTCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAacac                                                                                                                                                     <  4:111872/149‑1 (MQ=255)
  cATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGc                                                                                                                                                   >  3:133799/1‑149 (MQ=255)
   aTGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCt                                                                                                                                                  >  4:9117/1‑149 (MQ=255)
     gATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTaa                                                                                                                                                <  1:215666/149‑1 (MQ=255)
      aTTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACATGCTaaa                                                                                                                                               <  1:318002/149‑1 (MQ=255)
          cagcagTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATccc                                                                                                                                           <  4:403411/149‑1 (MQ=255)
              agTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACgg                                                                                                                                       >  1:50642/1‑149 (MQ=255)
               gTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACgga                                                                                                                                      <  4:428625/149‑1 (MQ=255)
                ttGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACggag                                                                                                                                     <  2:50642/149‑1 (MQ=255)
                    cTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAAc                                                                                                                                 >  1:356906/1‑149 (MQ=255)
                    cTGCGGGATCTAAACTCAAATGGTAAAGGAACGTGTGTAGAGCCTGAAACAGTCCGTTATTGATCCTGCTTTTGCGGTGGGAAATATTGCTGAATATCGGCGTAGCGTTGCAATTGATCTTATCAACATGCTAAATCCCACGGAGGAAc                                                                                                                                 <  2:331384/149‑1 (MQ=255)
                     tGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACa                                                                                                                                <  1:244399/149‑1 (MQ=255)
                        gggATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCt                                                                                                                             <  3:292943/149‑1 (MQ=255)
                          gATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTgt                                                                                                                           >  1:82566/1‑149 (MQ=255)
                           aTCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTgtt                                                                                                                          <  2:557769/149‑1 (MQ=255)
                           aTCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTgtt                                                                                                                          <  4:183692/149‑1 (MQ=255)
                            tCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTgttg                                                                                                                         >  3:83652/1‑149 (MQ=255)
                             cTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTgttgg                                                                                                                        <  3:341640/149‑1 (MQ=255)
                               acacTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTgttggtt                                                                                                                      <  4:268496/149‑1 (MQ=255)
                                 acTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACATGCTAAATCCCACGGAGGAACAGCTGttggttgg                                                                                                                    >  4:42839/1‑149 (MQ=255)
                                 acTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGttggttgg                                                                                                                    <  1:462462/149‑1 (MQ=255)
                                  cTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTtggttggt                                                                                                                   <  4:441369/149‑1 (MQ=255)
                                   tCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTggttggtt                                                                                                                  >  3:493741/1‑149 (MQ=255)
                                   tCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTggttggtt                                                                                                                  >  2:317093/1‑149 (MQ=255)
                                   tCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTggttggtt                                                                                                                  >  3:170868/1‑149 (MQ=255)
                                   tCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTggttggtt                                                                                                                  >  3:276528/1‑149 (MQ=255)
                                     aaaTGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTg                                                                                                                >  3:325330/1‑149 (MQ=255)
                                          gTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGcac                                                                                                           <  3:9117/149‑1 (MQ=255)
                                           tAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGcaca                                                                                                          >  2:273671/1‑149 (MQ=255)
                                               ggAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCtt                                                                                                      >  2:199396/1‑149 (MQ=255)
                                               ggAACGGGTGTATAGACTGAAACAGTGCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCtt                                                                                                      >  2:245818/1‑149 (MQ=255)
                                                    gggTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGcc                                                                                                 <  2:516862/149‑1 (MQ=255)
                                                           taGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACgga                                                                                                                                      <  1:477509/105‑1 (MQ=255)
                                                           taGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACgga                                                                                                                                      >  2:477509/1‑105 (MQ=255)
                                                                tGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTcc                                                                                     <  4:145331/149‑1 (MQ=255)
                                                                   aaCAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTccccc                                                                                  <  4:83652/149‑1 (MQ=255)
                                                                     cAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCctc                                                                                <  2:341872/149‑1 (MQ=255)
                                                                        tCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTa                                                                             >  3:235660/1‑149 (MQ=255)
                                                                        tCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTa                                                                             >  3:61464/1‑149 (MQ=255)
                                                                        tCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTa                                                                             >  1:287499/1‑149 (MQ=255)
                                                                               ttGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTAc                                                                      >  4:63439/1‑149 (MQ=255)
                                                                                 gATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGa                                                                    >  1:442640/1‑149 (MQ=255)
                                                                                   tACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACATGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATg                                                                  <  4:319590/149‑1 (MQ=255)
                                                                                      tGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCg                                                               <  2:26629/149‑1 (MQ=255)
                                                                                                    cAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGAtt                                                 >  2:438651/1‑149 (MQ=255)
                                                                                                     aaaTATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGAttt                                                >  1:33288/1‑149 (MQ=255)
                                                                                                            gCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGa                                         <  2:356906/149‑1 (MQ=255)
                                                                                                            gCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGa                                         >  1:124061/1‑149 (MQ=255)
                                                                                                                  tatCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATg                                   >  1:233837/1‑149 (MQ=255)
                                                                                                                  tatCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATg                                   >  2:132188/1‑149 (MQ=255)
                                                                                                                  tatCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGAAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATg                                   >  2:386933/1‑149 (MQ=255)
                                                                                                                   atCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGt                                  >  1:310215/1‑149 (MQ=255)
                                                                                                                   atCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGt                                  >  1:84293/1‑149 (MQ=255)
                                                                                                                         gTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAg                            <  2:279894/149‑1 (MQ=255)
                                                                                                                              gTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTaaa                       <  2:33288/149‑1 (MQ=255)
                                                                                                                              gTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTAGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTaaa                       >  2:174104/1‑149 (MQ=255)
                                                                                                                               tGGCAATTGATCTTATCAACATGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAg                      <  3:42839/149‑1 (MQ=255)
                                                                                                                                  cAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAg                   >  2:498959/1‑149 (MQ=255)
                                                                                                                                      tGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCg               <  2:495534/149‑1 (MQ=255)
                                                                                                                                        aTCTTATCAACATGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCgag             >  4:323671/1‑149 (MQ=255)
                                                                                                                                        aTCTTATCAACATGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCgag             >  4:88627/1‑149 (MQ=255)
                                                                                                                                         tCTTATCAACACGCTAAATCGCACGGAAGAACAACTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCAAGTGACGAATGTATACAGGCAAAGGAGTGCgaga            >  2:41669/1‑149 (MQ=255)
                                                                                                                                           ttATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGaaa          >  1:97744/1‑149 (MQ=255)
                                                                                                                                           ttATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGaaa          >  2:56930/1‑149 (MQ=255)
                                                                                                                                           ttATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGaaa          >  2:92273/1‑149 (MQ=255)
                                                                                                                                             aTCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGaaaaa        <  2:372668/149‑1 (MQ=255)
                                                                                                                                                aaCACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAgg                     <  2:409380/133‑1 (MQ=255)
                                                                                                                                                aaCACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAgg                     >  1:409380/1‑133 (MQ=255)
                                                                                                                                                 acgcGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAattat    >  3:557841/4‑149 (MQ=255)
                                                                                                                                                  cgcGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAattatt   <  4:557841/147‑1 (MQ=255)
                                                                                                                                                   acGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATtg  >  2:20218/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAACAGCTGTTGGTTGGTTTGCGCACAGCTTACGCCGACTCGGTTTCCCCCTCTTATTTGTACGATGTCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTG  >  NC_000913/875383‑875678

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