Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,767,817 C→T 6.5% T619I (ACC→ATC)  PP_1577 → lambda family tail tape meausure protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,767,8170CT 6.5% 183.0 / 2.0 62T619I (ACC→ATC) PP_1577lambda family tail tape meausure protein
Reads supporting (aligned to +/- strand):  ref base C (33/24);  new base T (2/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

AGGAGAAGTATTACCAGGATGTCGACGCGATGCAGCAGAACGGCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTGGCCATGGCCGGCCTGAACGCATTCACCGCGA  >  NC_002947/1767673‑1767960
                                                                                                                                                |                                                                                                                                                
aGGAGAAGTATTACCAGGATGTCGACGCGATGCAGCAGAACGGCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGAc                                                                                                                                              <  2:302217/149‑1 (MQ=255)
  gagaAGTATTACCAGGATGTCGACGCGGTGCAGCAGAACGGCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGAGGCCTTCAGCTCGTTGACCGACAt                                                                                                                                            >  1:219132/1‑149 (MQ=255)
     aaGTATTACCAGGATGTCGACGCGATGCAGCAGAACGGCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCt                                                                                                                                         >  4:42460/1‑149 (MQ=255)
       gTATTACCAGGATGTCGACGCGATGCAGCAGAACGGCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACGTGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGa                                                                                                                                       >  6:137462/1‑149 (MQ=255)
        tATTACCAGGATGTCGACGCGATGCAGCAGAACGGCACCACTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGAc                                                                                                                                      >  4:108002/1‑149 (MQ=255)
          ttACCAGGATGTCGACGCGATGCAGCAGAACGGCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGAcc                                                                                                                                     <  5:284258/148‑1 (MQ=255)
                ggATGTCGACGCGATGCAGCAGAACGGCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATggg                                                                                                                              >  1:33635/1‑149 (MQ=255)
                 gATGTCGACGCGATGCAGCAGAACGGCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGc                                                                                                                             >  2:75759/1‑149 (MQ=255)
                  aTGTCGACGCGATGCAGCAGAACGGCACCACTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCt                                                                                                                            <  3:108002/149‑1 (MQ=255)
                     tCGACGCGATGCAGCAGAACGGCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTga                                                                                                                          <  6:14542/148‑1 (MQ=255)
                       gACGCGATGCAGCAGAACGGCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGAc                                                                                                                       <  2:33635/149‑1 (MQ=255)
                          gcgATGCAGCAGAACGGCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACGTGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCgg                                                                                                                    <  5:137462/149‑1 (MQ=255)
                              tGCAGCAGAACGGCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAg                                                                                                                <  2:117370/149‑1 (MQ=255)
                                    agAACGGCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCTGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAAc                                                                                                          <  7:352080/149‑1 (MQ=255)
                                     gAACGTCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGATCGACATGCTGACCGAATGGGCTGAGACCGGGAATCTCAACg                                                                                                         >  8:75222/1‑149 (MQ=255)
                                     gAACGTCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGATCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACg                                                                                                         >  8:80586/1‑149 (MQ=255)
                                        cGGCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCa                                                                                                      <  3:106979/149‑1 (MQ=255)
                                         ggCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGcc                                                                                                       >  2:318902/1‑147 (MQ=255)
                                          gCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑g                                                                                                   <  3:42460/149‑1 (MQ=255)
                                            aCCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGa                                                                                                 >  6:280176/1‑149 (MQ=255)
                                               gCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACtt                                                                                              <  3:230784/149‑1 (MQ=255)
                                                cTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCCTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTc                                                                                             >  7:154417/1‑149 (MQ=255)
                                                cTGGCTTAATCAGCGGCTTCGCAACCCAGTCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTc                                                                                             <  8:373129/149‑1 (MQ=255)
                                                   gCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCt                                                                                          >  3:85910/1‑149 (MQ=255)
                                                    cTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTg                                                                                         >  7:30354/1‑149 (MQ=255)
                                                      caatccGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTg                                                                                         <  5:192662/145‑1 (MQ=255)
                                                              ggCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCat                                                                               <  2:151222/149‑1 (MQ=255)
                                                              ggCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGAAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCat                                                                               >  2:162440/1‑149 (MQ=255)
                                                                cTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCatc                                                                              >  3:188762/1‑148 (MQ=255)
                                                                   cGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCatccaa                                                                           <  8:30354/148‑1 (MQ=255)
                                                                    gCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCatccaat                                                                          <  1:318902/148‑1 (MQ=255)
                                                                        cccAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATc                                                                     >  7:49765/1‑149 (MQ=255)
                                                                         ccAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCatc                                                                              >  6:276330/1‑139 (MQ=255)
                                                                         ccAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCatc                                                                              <  5:276330/139‑1 (MQ=255)
                                                                          cAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGTGTTGGTTCGGCGGCCTTCAGCTCGTTGAGCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCgg                                                                   <  7:374169/149‑1 (MQ=255)
                                                                             gCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTAGCTGCAACCTTCAGCCAATCCATCGGTCa                                                                >  8:249467/1‑149 (MQ=255)
                                                                                    cTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCgctg                                                         >  3:273231/1‑149 (MQ=255)
                                                                                       ccATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAACGGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCgctg                                                         <  1:75759/147‑1 (MQ=255)
                                                                                       ccATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCgctgc                                                        <  4:85910/147‑1 (MQ=255)
                                                                                          tGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCg                                                   >  2:320657/1‑149 (MQ=255)
                                                                                                   aCAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATgcggcg                                          >  2:76755/1‑149 (MQ=255)
                                                                                                     aGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGcggcgg                                         >  3:353067/1‑148 (MQ=255)
                                                                                                       cAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGca                                       >  8:118050/1‑148 (MQ=255)
                                                                                                           aTGCAGAGCGTTGGTTCGGCGGCCTCCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGt                                  >  1:190416/1‑149 (MQ=255)
                                                                                                            tGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTg                                 >  6:11945/1‑149 (MQ=255)
                                                                                                            tGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTg                                 >  5:284916/1‑149 (MQ=255)
                                                                                                              cAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGATCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTGGc                               <  7:75222/149‑1 (MQ=255)
                                                                                                              cAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGATCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTGGc                               <  7:80586/149‑1 (MQ=255)
                                                                                                               agagCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCgg                                                                   >  1:162052/1‑112 (MQ=255)
                                                                                                               agagCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCgg                                                                   <  2:162052/112‑1 (MQ=255)
                                                                                                                gagCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTGGCCa                             >  8:289009/1‑149 (MQ=255)
                                                                                                                     ttGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTGGCCATggcc                        >  8:79759/1‑149 (MQ=255)
                                                                                                                             gcggcCTTCAGCTCGTTGACCGACATGCTTACCGAATGTGCTGAGACCGGTAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCTGCACAGGTGGCCATGGCCGGCCTGaa                <  8:154417/149‑1 (MQ=255)
                                                                                                                             gcggcCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTGGCCATGGCCGGCCTGaa                >  1:159524/1‑149 (MQ=255)
                                                                                                                                 ccTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTACATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTGGCCATGGCCGGCCTGAACGCa            <  1:162440/149‑1 (MQ=255)
                                                                                                                                   ttCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCGTCATCCAATCCATTGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTGGCCATGGCCGGCCTGAACGCAtt          >  6:169257/1‑149 (MQ=255)
                                                                                                                                     cAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTGGCCATGGCCGGCCTGAACGCAtt          <  6:394925/147‑1 (MQ=255)
                                                                                                                                       gCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTGGCCATGGCCGGCCTGAACGCATTCAcc      >  5:114357/1‑149 (MQ=255)
                                                                                                                                       gCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTACATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTGGCCATGGCCGGCCTGAACGCATTCAcc      >  3:316348/1‑149 (MQ=255)
                                                                                                                                           gTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTGGCCATGGCCGGCCTGAACGCATTCACCGCga  >  2:44237/1‑149 (MQ=255)
                                                                                                                                              gACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTGGCCATggccggc                     >  4:316418/1‑127 (MQ=255)
                                                                                                                                              gACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTGGCCATggccggc                     <  3:316418/127‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AGGAGAAGTATTACCAGGATGTCGACGCGATGCAGCAGAACGGCACCGCTGGCTTCATCAGCGGCTTCGCAACCCAGGCAGAAGCTGCCATGGATCTGTACAGCAACATGCAGAGCGTTGGTTCGGCGGCCTTCAGCTCGTTGACCGACATGCTGACCGAATGGGCTGAGACCGGGAAGCTCAACGCCAA‑GGACTTCGCTGCAACCTTCATCCAATCCATCGGGCACTCGCTGCTGGCGTATGCGGCGGCACAGGTGGCCATGGCCGGCCTGAACGCATTCACCGCGA  >  NC_002947/1767673‑1767960

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑