Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,111,009 A→G 9.5% S167S (AGT→AGC)  PP_5361 ← protein P47K

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,111,0090AG 9.5% 162.8 / 6.4 64S167S (AGT→AGC) PP_5361protein P47K
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base G (3/3);  total (32/32)
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.93e-01

CTTGGCCCATGACCATGGGTACGATCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCT  >  NC_002947/6110861‑6111154
                                                                                                                                                    |                                                                                                                                                 
cTTGGCCCATGACCATGGGTACGATCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCa                                                                                                                                                   <  6:83746/149‑1 (MQ=255)
   ggCCCATGACCATGGGTACGATCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCAc                                                                                                                                                  <  6:158260/147‑1 (MQ=255)
   ggCCCATGACCATGGGTACGATCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACt                                                                                                                                                 >  3:195608/1‑148 (MQ=255)
              aTGGGTACGATCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTc                                                                                                                                     <  6:44264/149‑1 (MQ=255)
               tGGGTACGATCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCg                                                                                                                                    >  8:174116/1‑149 (MQ=255)
                gggTACGATCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGt                                                                                                                                   >  2:159100/1‑149 (MQ=255)
                 ggTACGATCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGtc                                                                                                                                  <  3:245085/149‑1 (MQ=255)
                    aCGATCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGTGATTTGCTGATCGagc                                                                                                                                         >  1:12568/1‑137 (MQ=255)
                    aCGATCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGTGATTTGCTGATCGagc                                                                                                                                         <  2:12568/139‑3 (MQ=255)
                     cGATCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGtcttct                                                                                                                              <  5:276618/149‑1 (MQ=255)
                      gATCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGTGATTTGCTGATCGagcag                                                                                                                                       >  7:177494/1‑135 (MQ=255)
                      gATCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGTGATTTGCTGATCGagcag                                                                                                                                       <  8:177494/139‑5 (MQ=255)
                        tCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGTGATTTGCTGATCGag                                                                                                                                          <  7:154344/134‑2 (MQ=255)
                        tCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGTGATTTGCTGATCGag                                                                                                                                          >  8:154344/1‑133 (MQ=255)
                             gcTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGAACAGCAAATAACTGATCGACCGCTCGTCTTCTTCGCCCaa                                                                                                                      >  4:271598/1‑149 (MQ=255)
                              cTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAAc                                                                                                                     <  6:304081/149‑1 (MQ=255)
                                  ggCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGttt                                                                                                                 >  6:114815/1‑149 (MQ=255)
                                   gCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTc                                                                                                                <  5:114815/149‑1 (MQ=255)
                                         agaCTGCGCAAAATAGCGATGAGTTCTTCACGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCg                                                                                                          >  1:209193/1‑149 (MQ=255)
                                                  aaaaTAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAg                                                                                                 >  6:369742/1‑149 (MQ=255)
                                                   aaaTAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAgg                                                                                                <  7:174116/149‑1 (MQ=255)
                                                       aGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGt                                                                                            >  3:279002/1‑149 (MQ=255)
                                                         cGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTct                                                                                          >  4:284842/1‑149 (MQ=255)
                                                           aTGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTc                                                                                                                <  1:313404/125‑1 (MQ=255)
                                                           aTGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTc                                                                                                                >  2:313404/1‑125 (MQ=255)
                                                              aGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCtgctgc                                                                                     <  4:126177/149‑1 (MQ=255)
                                                               gTTCTTCGCGCTCATGTTGGCTGATCAGGTCTATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGcc                                                                                    <  3:271598/149‑1 (MQ=255)
                                                               gTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGcc                                                                                    <  5:143863/149‑1 (MQ=255)
                                                                   ttcGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGa                                                                                 >  1:317046/1‑148 (MQ=255)
                                                                     cgcgCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGa                                                                                 <  1:159100/146‑1 (MQ=255)
                                                                          tCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTcgcg                                                                         >  3:342324/1‑149 (MQ=255)
                                                                            aTGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCgcgc                                                                        <  4:342324/148‑1 (MQ=255)
                                                                            aTGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCa                                                                       <  4:195608/149‑1 (MQ=255)
                                                                                  gCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGt                                                                 >  4:352183/1‑149 (MQ=255)
                                                                                              gATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGa                                                     >  3:24676/1‑149 (MQ=255)
                                                                                                   tGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACg                                                >  7:210843/1‑149 (MQ=255)
                                                                                                    gCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGa                                                     <  4:279002/143‑1 (MQ=255)
                                                                                                     cTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGt                                              <  8:210843/149‑1 (MQ=255)
                                                                                                       cagcagcagCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGa                                            >  6:24230/1‑149 (MQ=255)
                                                                                                           agcagcaCATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCAt                                        <  2:209193/149‑1 (MQ=255)
                                                                                                             cagcaCATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTcgcg                                                                         <  6:320126/114‑1 (MQ=255)
                                                                                                             cagcaCATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTcgcg                                                                         >  5:320126/1‑114 (MQ=255)
                                                                                                                   aTCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGa                                                                                 >  2:148293/1‑100 (MQ=255)
                                                                                                                   aTCTCCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGa                                                                                 <  1:148293/100‑1 (MQ=255)
                                                                                                                      tGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATc                                                       <  4:18344/123‑1 (MQ=255)
                                                                                                                      tGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATc                                                       >  3:18344/1‑123 (MQ=255)
                                                                                                                       gCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGgt                                     >  2:165718/1‑140 (MQ=255)
                                                                                                                       gCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGgt                                     <  1:165718/140‑1 (MQ=255)
                                                                                                                        cGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTcgc                           <  4:24676/149‑1 (MQ=255)
                                                                                                                         gAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTcgcg                          >  8:317829/1‑149 (MQ=255)
                                                                                                                         gAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTcgcg                          >  5:202480/1‑149 (MQ=255)
                                                                                                                          aaCTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCgcgc                         <  2:317046/149‑1 (MQ=255)
                                                                                                                               gACCTGCTCGATCAGCAAATCTCTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCGGGCTGGCCAGGCTGTCTGCTGCCTGATAGTcgc                                                                          <  5:313966/95‑1 (MQ=255)
                                                                                                                                    gCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGa                                                                                 <  4:273922/83‑1 (MQ=255)
                                                                                                                                    gCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGa                                                                                 >  3:273922/1‑83 (MQ=255)
                                                                                                                                    gCTCGATCAGCAAATCACTGAACGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGa                >  4:286680/1‑148 (MQ=255)
                                                                                                                                          tCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTg         <  5:24230/149‑1 (MQ=255)
                                                                                                                                          tCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTg         <  3:352183/149‑1 (MQ=255)
                                                                                                                                            aGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACCGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGc        >  1:322675/1‑148 (MQ=255)
                                                                                                                                             gCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATc                                                       >  6:284510/1‑100 (MQ=255)
                                                                                                                                             gCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATc                                                       <  5:284510/100‑1 (MQ=255)
                                                                                                                                             gCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGAcc                                          <  4:277528/113‑1 (MQ=255)
                                                                                                                                             gCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGAcc                                          >  3:277528/1‑113 (MQ=255)
                                                                                                                                               aaaTCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACCGTGACCATGGTGTCCAGTCCCGCCATATCGGACAGACTGCGAcc    >  7:153376/1‑149 (MQ=255)
                                                                                                                                                 aTCACTGATCGACCGCTCGTCTTCTTCGCCAAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGTGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCt  >  3:49346/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CTTGGCCCATGACCATGGGTACGATCTGCGCTCGGGCATTGAGACTGCGCAAAATAGCGATGAGTTCTTCGCGCTCATGTTGGCTGATCAGGTCGATCTTGCTCAGCAGCAGCACATCTGCGAACTCGACCTGCTCGATCAGCAAATCACTGATCGACCGCTCGTCTTCTTCGCCCAACGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCT  >  NC_002947/6110861‑6111154

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