Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,731,846 A→G 7.9% R465R (CGT→CGC)  PP_3299 ← lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,731,8460AG 7.9% 176.3 / 5.4 63R465R (CGT→CGC) PP_3299lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (30/28);  new base G (3/2);  total (33/30)
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.57e-01

ACTTGGATTACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACC  >  NC_002947/3731699‑3731992
                                                                                                                                                   |                                                                                                                                                  
aCTTGGATTACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCGTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCAc                                                                                                                                                   <  2:79236/149‑1 (MQ=255)
       ttACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCGCGtgctgc                                                                                                                                            >  8:210011/1‑149 (MQ=255)
        tACGCTGTATGAGGCAAATCCTCAACGTAGGTGAGTAGTTCTTCTGCACCAAGGCATTTACGGAGTCCCATCATCGCTGTGTGAATTATGTGTCAGTTGAGCACGCGCCTCCCAGGGGGAGGTAGGCAGCTTTAGCAGCACGCgctgc                                                                                                                                            <  4:232062/148‑1 (MQ=255)
          cgcgGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCGCGTgctgctat                                                                                                                                         <  7:210011/149‑3 (MQ=255)
                  gAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCTCCAAGGGATTTACGGAGTCCGATCATCGCTTTGCCTCATCTGAGTCAGTTGAGCACCCGCCAGGCTGGGCGAGGTAGGCTGCAATAGCAGCACGCGCTGCTTCAGTTTgc                                                                                                                                  <  4:161490/148‑1 (MQ=255)
                   aGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTgcg                                                                                                                                 >  8:259460/1‑148 (MQ=255)
                       cAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAAc                                                                                                                             >  8:9275/1‑148 (MQ=255)
                       cAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAAc                                                                                                                             <  7:9275/148‑1 (MQ=255)
                           ccTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTc                                                                                                                        >  3:269538/1‑149 (MQ=255)
                            cTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCg                                                                                                                       >  7:3352/1‑149 (MQ=255)
                                             cGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTga                                                                                                      <  4:85781/149‑1 (MQ=255)
                                                 ttctGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGAGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAg                                                                                                  >  3:127433/1‑149 (MQ=255)
                                                         cAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCGCGTgctgctattg                                                                                                                                       <  7:327094/104‑5 (MQ=255)
                                                         cAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCGCGTgctgctattg                                                                                                                                       >  8:327094/1‑100 (MQ=255)
                                                          aaGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATc                                                                                         >  5:71185/1‑149 (MQ=255)
                                                            ggCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCt                                                                                                        >  8:203533/1‑132 (MQ=255)
                                                            ggCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCt                                                                                                        <  7:203533/132‑1 (MQ=255)
                                                              cATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCg                                                                                     <  7:259460/149‑1 (MQ=255)
                                                               aTTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGt                                                                                    >  5:93055/1‑149 (MQ=255)
                                                               aTTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACAACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGt                                                                                    >  1:252956/1‑149 (MQ=255)
                                                                     gCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAgctg                                                                                               <  3:285027/132‑1 (MQ=255)
                                                                     gCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAgctg                                                                                               >  4:285027/1‑132 (MQ=255)
                                                                           ccATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCg                                                                                                                       <  7:109925/102‑1 (MQ=255)
                                                                           ccATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCg                                                                                                                       >  8:109925/1‑102 (MQ=255)
                                                                           ccATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCg                                                                                                                       <  1:249604/102‑1 (MQ=255)
                                                                           ccATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCg                                                                                                                       <  1:165638/102‑1 (MQ=255)
                                                                           ccATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCg                                                                                                                       >  2:165638/1‑102 (MQ=255)
                                                                           ccATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCg                                                                                                                       >  2:249604/1‑102 (MQ=255)
                                                                           ccATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTa                                                                        <  4:127433/149‑1 (MQ=255)
                                                                           ccATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTa                                                                        >  5:165852/1‑149 (MQ=255)
                                                                           ccATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGACGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTa                                                                        >  4:159590/1‑149 (MQ=255)
                                                                             atcatcGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAgctg                                                                                               >  1:349596/1‑124 (MQ=255)
                                                                             atcatcGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAgctg                                                                                               <  2:349596/124‑1 (MQ=255)
                                                                             atcatcGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATc                                                                      <  4:269538/149‑1 (MQ=255)
                                                                                    cTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTgg                                                               <  6:93055/149‑1 (MQ=255)
                                                                                         ccACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAg                                                                                                  >  8:249624/1‑109 (MQ=255)
                                                                                         ccACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAg                                                                                                  <  7:249624/109‑1 (MQ=255)
                                                                                           acaTCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACAcc                                                        >  5:149824/1‑149 (MQ=255)
                                                                                             aTCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAgct                                                                                                <  2:120055/107‑1 (MQ=255)
                                                                                                    tCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGaa                                               <  6:165852/149‑1 (MQ=255)
                                                                                                    tCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGACGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGaa                                               <  3:159590/149‑1 (MQ=255)
                                                                                                      aGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATg                                             >  6:300854/1‑149 (MQ=255)
                                                                                                        ttGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGc                                            >  7:317961/1‑148 (MQ=255)
                                                                                                             caaccGCCACCAAGGGGGTTGTATGGATGAATACCATCACGGGGTGCTTTAGTTTTCGCAAATTCTCGGAAATCCGTATTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGaa                                      <  2:119228/146‑1 (MQ=255)
                                                                                                              accCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAg                                     >  5:228366/1‑149 (MQ=255)
                                                                                                              accCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAg                                     >  4:153636/1‑149 (MQ=255)
                                                                                                              accCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAg                                     >  8:133008/1‑149 (MQ=255)
                                                                                                                ccGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAAt                                              <  1:323906/138‑1 (MQ=255)
                                                                                                                ccGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAAt                                              >  2:323906/1‑138 (MQ=255)
                                                                                                                ccGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGa                                   >  6:65118/1‑149 (MQ=255)
                                                                                                                       ccAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAAcgt                            >  1:95302/1‑149 (MQ=255)
                                                                                                                       ccAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAAcgt                            >  6:16555/1‑149 (MQ=255)
                                                                                                                          gggCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAAcgtcgt                         <  5:16555/149‑1 (MQ=255)
                                                                                                                            gCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTAc                                         <  4:81488/131‑1 (MQ=255)
                                                                                                                            gCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTAc                                         >  3:81488/1‑131 (MQ=255)
                                                                                                                               aGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTgctcg                    <  3:153636/149‑1 (MQ=255)
                                                                                                                                ggTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTgctcg                    <  7:133008/148‑1 (MQ=255)
                                                                                                                                  tAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTgctcgctc                 >  4:151206/1‑149 (MQ=255)
                                                                                                                                     gcagcaATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAg              <  6:149824/149‑1 (MQ=255)
                                                                                                                                      cagcaATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAgg             <  6:228366/149‑1 (MQ=255)
                                                                                                                                       agcaATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAgg             >  3:24817/1‑148 (MQ=255)
                                                                                                                                          tgaagcagcGCGTGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATa         >  4:290991/4‑149 (MQ=255)
                                                                                                                                          aaTAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATa         <  8:3352/149‑1 (MQ=255)
                                                                                                                                               cagcaCGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCa    >  4:201154/1‑149 (MQ=255)
                                                                                                                                                agcaCGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCAc   >  8:166554/1‑149 (MQ=255)
                                                                                                                                                 gcaCGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCAcc  <  3:151206/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACTTGGATTACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACC  >  NC_002947/3731699‑3731992

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