Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,010,864 G→T 6.1% L67M (CTG→ATG)  PP_0871 ← glycine betaine/carnitine/choline ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,010,8640GT 6.1% 209.7 / 3.6 66L67M (CTG→ATG) PP_0871glycine betaine/carnitine/choline ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base T (2/2);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

GCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCC  >  NC_002947/1010716‑1011011
                                                                                                                                                    |                                                                                                                                                   
gcgcCGATGCCCGGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAg                                                                                                                                                     >  2:213508/1‑149 (MQ=255)
     gATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGcgccg                                                                                                                                                >  4:353781/1‑149 (MQ=255)
             ggATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATc                                                                                                                                        >  1:5801/1‑149 (MQ=255)
                   tCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGGCGGCGGCGCTGGCCGTGGGCATCCCCGCCGGCAttgc                                                                                                                                      <  1:126742/145‑4 (MQ=37)
                   tCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGGCGGCGGCGCTGGCCGTGGGCATCCCCGCCGGCAttgc                                                                                                                                      >  2:126742/1‑142 (MQ=37)
                   tCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGc                                                                                                                                   >  1:262107/1‑148 (MQ=255)
                      gggCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAg                                                                                                                                >  3:74283/1‑148 (MQ=255)
                       ggCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCa                                                                                                                              <  5:238596/149‑1 (MQ=255)
                        gCGTTGGCCAGAACGGCCAGAGGGGGGAGGGGGGTGCCAACGGTGCAGAACTGCCTGACGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCGCTCGGCCAGCGCCGCCGCCATCGAAGCCAGCAc                                                                                                                             <  6:258504/149‑1 (MQ=255)
                        gCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCAc                                                                                                                             <  2:262107/149‑1 (MQ=255)
                             ggCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATa                                                                                                                        >  8:295774/1‑149 (MQ=255)
                               ccAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATAcc                                                                                                                      <  3:353781/149‑1 (MQ=255)
                               ccAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATAcc                                                                                                                      <  7:295774/149‑1 (MQ=255)
                                   aaCGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGt                                                                                                                  >  4:203738/1‑149 (MQ=255)
                                     cGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTgc                                                                                                                <  3:203738/149‑1 (MQ=255)
                                        ccAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCt                                                                                                             >  3:105010/1‑149 (MQ=255)
                                        ccAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGCGCTGGCCGTGGGCATccccgccgg                                                                                                                                            <  1:114381/118‑2 (MQ=37)
                                        ccAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGCGCTGGCCGTGGGCATccccgccgg                                                                                                                                            >  2:114381/1‑117 (MQ=37)
                                            agGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCa                                                                                                         <  2:5801/149‑1 (MQ=255)
                                               ggggATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGt                                                                                                      >  8:326897/1‑149 (MQ=255)
                                                            ccAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAg                                                                                                        >  2:353251/1‑134 (MQ=255)
                                                            ccAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAg                                                                                                        <  1:353251/134‑1 (MQ=255)
                                                             cAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTa                                                                                                     >  1:108876/1‑136 (MQ=255)
                                                             cAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGAtcg                                                                                        >  4:64839/1‑149 (MQ=255)
                                                             cAACGTTGAAGAACTGCATGAACCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTa                                                                                                     <  2:108876/136‑1 (MQ=255)
                                                              aaCGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGAtcgt                                                                                       >  5:366292/1‑149 (MQ=255)
                                                                      agaaCTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCATTACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATaa                                                                               <  1:249023/149‑1 (MQ=255)
                                                                             cATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGa                                                                        >  2:171597/1‑149 (MQ=255)
                                                                                   gcgcTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCgtgt                                                                  <  4:105010/149‑1 (MQ=255)
                                                                                    cgcTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTc                                                                 <  7:92977/149‑1 (MQ=255)
                                                                                     gcTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCCCCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTc                                                                 >  8:69081/1‑148 (MQ=255)
                                                                                        cGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCgtgt                                                                  >  1:272498/1‑144 (MQ=255)
                                                                                        cGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCgtgt                                                                  <  2:272498/144‑1 (MQ=255)
                                                                                            tttGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGa                                                         >  3:309533/1‑149 (MQ=255)
                                                                                              tGTCGACCCGGTGCGGTCGCCTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCGTGCAGGTAGAACCCAAGATCGTCGCGAGACCGCGCGATCGTGTCGATGCCGCTc                                                       >  6:362244/1‑149 (MQ=255)
                                                                                              tGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGa                                                         <  1:171597/147‑1 (MQ=255)
                                                                                               gTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATc                                                       <  4:309533/148‑1 (MQ=255)
                                                                                               gTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATc                                                       <  6:366292/148‑1 (MQ=255)
                                                                                               gTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATc                                                       <  3:64839/148‑1 (MQ=255)
                                                                                                tCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGc                                                            >  5:340454/1‑142 (MQ=255)
                                                                                                tCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGc                                                            <  6:340454/142‑1 (MQ=255)
                                                                                                  gACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGc                                                            <  1:87709/140‑1 (MQ=255)
                                                                                                  gACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGc                                                            >  2:87709/1‑140 (MQ=255)
                                                                                                  gACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAg                                                    <  7:326897/148‑1 (MQ=255)
                                                                                                          cgggTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCCGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACcagc                                           <  7:189030/147‑1 (MQ=255)
                                                                                                           cGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACcagc                                           <  1:213508/148‑1 (MQ=255)
                                                                                                             gTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAgg                                        <  4:74283/149‑1 (MQ=255)
                                                                                                              tCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCAAGTGGACCAGCAggg                                       >  1:242480/1‑149 (MQ=255)
                                                                                                                   ttAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAg                                   >  1:153960/1‑148 (MQ=255)
                                                                                                                   ttAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAg                                   >  6:7075/1‑148 (MQ=255)
                                                                                                                   ttAACGCAATGCCGGCGGGGATGCCCACGGCCAGCACCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAg                                   >  8:219155/1‑148 (MQ=255)
                                                                                                                     aaCGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGAt                                >  4:232201/1‑149 (MQ=255)
                                                                                                                        gCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGAc                             >  8:287105/1‑149 (MQ=255)
                                                                                                                        gCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGAc                             >  2:201951/1‑149 (MQ=255)
                                                                                                                        gCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGAc                             >  4:148967/1‑149 (MQ=255)
                                                                                                                          aaTGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGCTGCCGATCCAGTGGACCAGCAGGGCCAGGCTGAc                             >  2:154763/1‑147 (MQ=255)
                                                                                                                          aaTGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGa                           <  5:7075/149‑1 (MQ=255)
                                                                                                                                ggcggGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCAcgg                     >  3:106752/1‑149 (MQ=255)
                                                                                                                                     ggATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGcc                >  2:94466/1‑149 (MQ=255)
                                                                                                                                         gCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATc            >  3:252502/1‑149 (MQ=255)
                                                                                                                                         gCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATc            >  3:247968/1‑149 (MQ=255)
                                                                                                                                           ccACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATccc          <  8:14343/149‑1 (MQ=255)
                                                                                                                                           ccACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATcc           >  7:14343/1‑148 (MQ=255)
                                                                                                                                             aCGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCaa        <  3:232201/149‑1 (MQ=255)
                                                                                                                                                 ccAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGc    <  4:252502/149‑1 (MQ=255)
                                                                                                                                                 ccAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGc    <  4:247968/149‑1 (MQ=255)
                                                                                                                                                   aGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGccc  <  4:106752/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAGCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCGCGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCC  >  NC_002947/1010716‑1011011

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