Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,906,643 A→T 9.8% L53Q (CTG→CAG)  spuF ← spermidine/putrescine ABC transporter ATP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,906,6430AT 9.8% 168.8 / 7.7 61L53Q (CTG→CAG) spuFspermidine/putrescine ABC transporter ATP binding protein
Reads supporting (aligned to +/- strand):  ref base A (31/24);  new base T (3/3);  total (34/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.68e-01

GACTGGAACATCATGTTGATTGGCCGCTCGTAGGGCGGCATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCTTTTTTATAGGCACCGGAGGC  >  NC_002947/5906496‑5906791
                                                                                                                                                   |                                                                                                                                                        
gACTGGAACATCATGTTGATTGGCCGCTCGTAGGGCGGCATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCcagcag                                                                                                                                                         >  6:321794/1‑149 (MQ=255)
gACTGGAACATCATGTTGATTGGCCGCTCGTAGGGCGGCATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCcagcag                                                                                                                                                         <  1:152298/149‑1 (MQ=255)
         atcatGTTGATTGGCCGCTCGTAGGGCGGCATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCTGCTGGGCggctcc                                                                                                                                                <  7:290804/149‑6 (MQ=255)
             tGTTGATTGGCCGCTCGTAGGGCGGCATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCa                                                                                                                                            >  8:229506/1‑149 (MQ=255)
              gTTGATTGGCCGCTCGTAGGGCGGCATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCAc                                                                                                                                           >  4:5564/1‑149 (MQ=255)
               ttGATTGGCCGCTCGTAGGGCGGCATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAATATCTCAcc                                                                                                                                          >  2:106839/1‑149 (MQ=255)
                   ttGGCCGCTCGTAGGGCGGCATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGTATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTg                                                                                                                                      <  6:199489/149‑1 (MQ=255)
                   ttGGCCGCTCGTAGGGCGGCATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTg                                                                                                                                      >  1:186520/1‑149 (MQ=255)
                     ggCCGCTCGTAGGGCGGCATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGc                                                                                                                                     <  7:59917/148‑1 (MQ=255)
                       ccGCTCGTAGGGCGGCATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGa                                                                                                                                  >  4:199204/1‑149 (MQ=255)
                             gTAGGGCGGCATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTc                                                                                                                           <  7:233816/147‑1 (MQ=255)
                              tAGGGCGGTATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCGAg                                                                                                                        <  8:69784/149‑1 (MQ=255)
                                  gcggcATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGac                                                                                                                    <  5:321794/149‑1 (MQ=255)
                                       aTGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATc                                                                                                               <  3:5564/149‑1 (MQ=255)
                                       aTGTCGGTGATGTCGACGCCATCGAGGAAGATACGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACAtg                                                                                                               >  2:73956/1‑148 (MQ=255)
                                                   tCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTAc                                                                                                   >  6:339653/1‑149 (MQ=255)
                                                          cATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTtcg                                                                                            >  8:319366/1‑149 (MQ=255)
                                                           aTCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGaa                                                                                                                                                   >  1:344654/1‑96 (MQ=255)
                                                           aTCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGaa                                                                                                                                                   <  2:344654/96‑1 (MQ=255)
                                                           aTCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTtcgt                                                                                           >  5:29160/1‑149 (MQ=255)
                                                               aGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGaa                                                                                       <  5:339653/149‑1 (MQ=255)
                                                                 gAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACt                                                                                     >  3:81700/1‑149 (MQ=255)
                                                                   aGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAt‑‑‑tt                                                                                                                            <  7:266118/108‑1 (MQ=255)
                                                                   aGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAt‑‑‑tt                                                                                                                            >  8:266118/1‑108 (MQ=255)
                                                                       ccGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCAc                                                                                                         <  7:190534/123‑1 (MQ=255)
                                                                       ccGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCAc                                                                                                         >  8:190534/1‑123 (MQ=255)
                                                                        cGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAATATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGt                                                                              <  1:106839/149‑1 (MQ=255)
                                                                         gCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGt                                                                              >  6:39122/1‑148 (MQ=255)
                                                                         gCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGt                                                                              >  4:78334/1‑148 (MQ=255)
                                                                           ccTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGAc                                                                          >  4:258184/1‑149 (MQ=255)
                                                                                       cgcTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCCGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCtt                                                              <  3:199204/149‑1 (MQ=255)
                                                                                            gAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTc                                                                                                            <  1:137104/99‑1 (MQ=255)
                                                                                            gAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTc                                                                                                            >  2:137104/1‑99 (MQ=255)
                                                                                            gAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGAc                                                           <  7:217086/147‑1 (MQ=255)
                                                                                                  ggCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGGCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCt                                                   >  3:29292/1‑149 (MQ=255)
                                                                                                   gCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTg                                                  >  1:202303/1‑149 (MQ=255)
                                                                                                        cATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGATAGGGGACAGGGACCCATCCTACACGGCAACCGTTTCGTCGAACTTTTTCGTTACCCCGTCGAAAATGACCAGGACCTGCt                                                >  1:294161/1‑149 (MQ=255)
                                                                                                         aTGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTggg                                            <  4:81700/149‑1 (MQ=255)
                                                                                                            cgcAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGa                                                            <  8:308764/130‑1 (MQ=255)
                                                                                                            cgcAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGa                                                            >  7:308764/1‑130 (MQ=255)
                                                                                                              cagcaAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCCCCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCATGACAACCACCAGCATgggatcct                                       >  6:108964/1‑144 (MQ=255)
                                                                                                               agcaAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGaa                                                                                       <  5:73896/101‑1 (MQ=255)
                                                                                                               agcaAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGaa                                                                                       >  6:73896/1‑101 (MQ=255)
                                                                                                               agcaAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTg                                      >  2:364350/1‑149 (MQ=255)
                                                                                                                 caAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGc                                    <  6:29160/149‑1 (MQ=255)
                                                                                                                        gACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATGGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCt                              <  1:73956/149‑1 (MQ=255)
                                                                                                                             gCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAggg                        >  2:2700/1‑149 (MQ=255)
                                                                                                                               cGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTg                                                  <  5:269037/121‑1 (MQ=255)
                                                                                                                               cGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTg                                                  >  8:199654/1‑121 (MQ=255)
                                                                                                                               cGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTg                                                  <  7:199654/121‑1 (MQ=255)
                                                                                                                               cGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTg                                                  >  6:269037/1‑121 (MQ=255)
                                                                                                                               cGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCt                      >  7:166835/1‑149 (MQ=255)
                                                                                                                               cGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCt                      >  5:75653/1‑149 (MQ=255)
                                                                                                                                    ccGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGc                                                       >  3:295432/1‑111 (MQ=255)
                                                                                                                                    ccGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGc                                                       <  4:295432/111‑1 (MQ=255)
                                                                                                                                    ccGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCtttttt                 <  2:202303/149‑1 (MQ=255)
                                                                                                                                     gatcttcgccCTGCTGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCTTTTTta                >  3:206971/7‑149 (MQ=255)
                                                                                                                                     gatcttcgccCTGCTGGGCGAAGATATCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCTTTTTta                <  4:206971/143‑1 (MQ=255)
                                                                                                                                      atcttcgccCTGCTGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCTTTTTtat               >  3:134700/6‑149 (MQ=255)
                                                                                                                                       tcttcgccCTGCTGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTg                                      >  2:62085/5‑125 (MQ=255)
                                                                                                                                       tcttcgccCTGCTGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTg                                      <  1:62085/121‑1 (MQ=255)
                                                                                                                                       gatccgccCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCTTTTTtata              <  5:39122/146‑1 (MQ=255)
                                                                                                                                            gccCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCTTTTTTATAGGCAc         >  7:79781/1‑149 (MQ=255)
                                                                                                                                                  cagGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCTAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCTTTTTTATAGGCACCGGAgg   >  2:39828/1‑149 (MQ=255)
                                                                                                                                                   agGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCTTTTTTATAGGCACCGGAGGc  <  6:75653/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                        
GACTGGAACATCATGTTGATTGGCCGCTCGTAGGGCGGCATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCTTTTTTATAGGCACCGGAGGC  >  NC_002947/5906496‑5906791

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