Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,906,640 A→T 9.8% L54Q (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,6400AT 9.8% 171.4 / 7.3 61L54Q (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 = 8.85e-01

GACTGGAACATCATGTTGATTGGCCGCTCGTAGGGCGGCATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCTTTTTTATAGGCAC  >  NC_002947/5906496‑5906784
                                                                                                                                                |                                                                                                                                                    
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                                                                                                                                            <  2:344654/96‑1 (MQ=255)
                                                           aTCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGaa                                                                                                                                            >  1:344654/1‑96 (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                                                                                >  6:73896/1‑101 (MQ=255)
                                                                                                               agcaAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGaa                                                                                <  5:73896/101‑1 (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                                           >  8:199654/1‑121 (MQ=255)
                                                                                                                               cGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTg                                           <  7:199654/121‑1 (MQ=255)
                                                                                                                               cGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTg                                           <  5:269037/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                                                <  4:295432/111‑1 (MQ=255)
                                                                                                                                    ccGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGc                                                >  3:295432/1‑111 (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                               <  1:62085/121‑1 (MQ=255)
                                                                                                                                       tcttcgccCTGCTGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTg                               >  2:62085/5‑125 (MQ=255)
                                                                                                                                       gatccgccCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCTTTTTtata       <  5:39122/146‑1 (MQ=255)
                                                                                                                                            gccCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCTTTTTTATAGGCAc  >  7:79781/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GACTGGAACATCATGTTGATTGGCCGCTCGTAGGGCGGCATGTCGGTGATGTCGACGCCATCGAGGAAGATCCGCCCTTCGGACGGGCGCTCGAAGCCGGCCAGCATGCGCAGCAAGGTGGACTTGCCGGAACCGGAGCCGCCCAGCAGGGCGAAGATCTCACCCTTGCGGAT‑‑‑TTCCAGGGACACATCGTCCACGGCTACCGTTTCGTCGAACTTTTTCGTGACCCGGTCGATCTTGACCAGCACCTGCTTGGGTTGCTGGCCACCCTCGAGGGCTTTTTTATAGGCAC  >  NC_002947/5906496‑5906784

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: