Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 543,727 A→G 6.7% N656S (AAC→AGC)  rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947543,7270AG 6.7% 214.9 / 4.1 74N656S (AAC→AGC) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base A (32/37);  new base G (2/3);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGA  >  NC_002947/543581‑543873
                                                                                                                                                  |                                                                                                                                                  
gTGCTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACg                                                                                                                                                  <  6:176387/149‑1 (MQ=255)
   cTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTa                                                                                                                                               >  8:36352/1‑149 (MQ=255)
    tACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTAc                                                                                                                                              <  3:137225/149‑1 (MQ=255)
      ccGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCCACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTAc                                                                                                                                              <  6:17379/147‑1 (MQ=255)
      ccGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTAc                                                                                                                                              <  1:337010/147‑1 (MQ=255)
           tggttggtCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAGCGTTACCGAtgat                                                                                                                                       <  4:307332/149‑2 (MQ=255)
           tggttggtCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaa                                                                                                                                       <  3:302325/149‑1 (MQ=255)
             gttggtCTGTAACCCTCCCTTATCTTCGCCTTCCAACTTTTGTTCTCCCGTTTCGCATTATCCTCCATTTCCGGCGTTTACATCGGTGTTAACGACTTCGTTATCCCGGACTAGAAAGCCCGCATCATCGGTAACGCTACCGAtgatgt                                                                                                                                     <  8:317289/149‑1 (MQ=255)
             gttggtCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaagt                                                                                                                                     <  5:304363/149‑1 (MQ=255)
                 gtCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaagtgaa                                                                                                                                  <  4:74064/148‑1 (MQ=255)
                 gtCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaagtgaa                                                                                                                                  <  4:30845/148‑1 (MQ=255)
                  tCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAGCGTTACCGAtg                                                                                                                                         >  7:119885/1‑140 (MQ=255)
                  tCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAGCGTTACCGAtg                                                                                                                                         <  8:119885/140‑1 (MQ=255)
                           gaCCGTTATCTTCGCCGACCTACTGATTTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACTACTTCGTTATCCCGGACGATAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCgaga                                                                                                                       <  1:267779/149‑1 (MQ=255)
                                               taCTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTg                                                                                                   <  3:143593/148‑1 (MQ=255)
                                               aaCTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTg                                                                                                   <  7:191522/149‑1 (MQ=255)
                                                 cTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTg                                                                                                   >  4:143593/1‑147 (MQ=255)
                                                      gTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAAcc                                                                                             >  8:29278/1‑148 (MQ=255)
                                                        acacCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAg                                                                                          >  6:229607/1‑149 (MQ=255)
                                                         cacCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAgg                                                                                         <  2:320896/149‑1 (MQ=255)
                                                          acCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAgg                                                                                         >  5:166487/1‑148 (MQ=255)
                                                          acCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAgg                                                                                         <  6:166487/148‑1 (MQ=255)
                                                            cGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTg                                                                                                   >  2:196960/1‑136 (MQ=255)
                                                            cGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTg                                                                                                   <  1:196960/136‑1 (MQ=255)
                                                            cGGTTTCGCATACTCGACCATTTCCGGCCTTTCCATCGGTGTTAACGACTTCGTTATCCCCGACGAGAAAGCCAGCATCATCCGTAAAGCAACCGAAGAAGTGAAGGAAATCGAGAGCCAGAACGCCTCCGGCCTGCTAACCCAgcgc                                                                                       >  1:61838/1‑148 (MQ=255)
                                                              gTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCga                                                                                                                         >  4:186245/1‑112 (MQ=255)
                                                                  cGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAg                                                                                 <  4:55235/148‑1 (MQ=255)
                                                                   gCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTa                                                                               <  7:29278/149‑1 (MQ=255)
                                                                   gCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGACCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCgaga                                                                                                                       <  1:147963/109‑1 (MQ=255)
                                                                   gCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGACCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCgaga                                                                                                                       >  2:147963/1‑109 (MQ=255)
                                                                    cATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGc                                                                                       <  5:121072/140‑1 (MQ=255)
                                                                    cATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGc                                                                                       >  6:121072/1‑140 (MQ=255)
                                                                    cATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTac                                                                              <  3:8649/149‑1 (MQ=255)
                                                                       aCTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTacaac                                                                           <  7:36352/149‑1 (MQ=255)
                                                                         tCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTacaac                                                                           >  4:264421/1‑147 (MQ=255)
                                                                           gACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAg                                                                        <  1:177751/148‑1 (MQ=255)
                                                                           gACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAg                                                                        >  2:177751/1‑148 (MQ=255)
                                                                             ccATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCACGGACGAGAAAGCCCGCATCATCGGTAACGCTACAGATGAAGTGAAGGAAATCGAGAGCAAGTACGCCTCCGGCCTCGTAACCCAGGGCGAGAAGTACAACAAgga                                                                      >  2:96952/1‑147 (MQ=255)
                                                                               ctttCCGGCGCTTCCCTCGGTGTTACCGACTTCGTTATCCCGGACGAGCAAGCCCGCCTCATCGGTAACGCTACCGATGAAGTGACCGAACTCGAGAGCCACTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAACGTCAt                                                                   <  4:241446/148‑1 (MQ=255)
                                                                                      gCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTtg                                                            <  5:229607/149‑1 (MQ=255)
                                                                                        gTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCga                                                                                     <  8:155812/122‑1 (MQ=255)
                                                                                        gTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCga                                                                                     >  7:155812/1‑122 (MQ=255)
                                                                                            ccATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTc                                                       >  8:90665/1‑148 (MQ=255)
                                                                                              aacGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCg                                                      >  4:42277/3‑147 (MQ=255)
                                                                                              aTCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCg                                                      >  1:100448/1‑147 (MQ=255)
                                                                                              aTCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGaa                                                    <  3:264421/149‑1 (MQ=255)
                                                                                              aTCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGaa                                                    <  3:42277/149‑1 (MQ=255)
                                                                                              aTCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCAAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCg                                                      >  5:95186/1‑147 (MQ=255)
                                                                                                    gtTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAg                                                                                 >  5:333562/1‑114 (MQ=255)
                                                                                                    gtTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAg                                                                                 <  6:333562/114‑1 (MQ=255)
                                                                                                     ttAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGaa                                              <  7:90665/148‑1 (MQ=255)
                                                                                                             ttCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTg                                                                                                   >  8:124785/1‑87 (MQ=255)
                                                                                                             ttCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTg                                                                                                   <  7:124785/87‑1 (MQ=255)
                                                                                                             ttCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAgt                                     <  1:96952/149‑1 (MQ=255)
                                                                                                                 ttATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTtg                                                            >  6:19224/1‑122 (MQ=255)
                                                                                                                 ttATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTcc                                 <  2:100448/149‑1 (MQ=255)
                                                                                                                      cctgACGAGAAAGCCCGCATCATAGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGACTACGGCCTGGTAAACCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCTAAGGAGAACGAAGAAGTGTACAAGGc                            >  6:204681/4‑149 (MQ=255)
                                                                                                                      ccGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGCAGGAGAACGACGAAGTCTCAAAGGc                            >  8:323400/1‑149 (MQ=255)
                                                                                                                      ccGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGc                            >  8:327429/1‑149 (MQ=255)
                                                                                                                            gagaAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAAGAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCgatgat                      >  1:296353/1‑149 (MQ=255)
                                                                                                                            gagaAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTtg                                                            <  4:3665/111‑1 (MQ=255)
                                                                                                                            gagaAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTtg                                                            >  3:3665/1‑111 (MQ=255)
                                                                                                                            gagaAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCgatgat                      >  5:248498/1‑149 (MQ=255)
                                                                                                                            gagaAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCga                          >  8:143848/1‑145 (MQ=255)
                                                                                                                            gagaAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCga                          <  7:143848/145‑1 (MQ=255)
                                                                                                                               aaaGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATgg                    >  1:99544/1‑148 (MQ=255)
                                                                                                                               aaaGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGc                   >  1:291835/1‑149 (MQ=255)
                                                                                                                                       catcatCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACctct           <  2:296353/149‑1 (MQ=255)
                                                                                                                                        ttcatcGTTAGCGTTACCGATGAAGTGAATGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAACGTCATCGACTTGTGGTCGAAGGCGAACGACGaa                                       <  4:359711/119‑1 (MQ=255)
                                                                                                                                        ttcatcGGTAGCGTTACCGATGATGTGAAGGAAATCGAGAGCCAGTACGACTCCGGCCTGGTAACCAAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGaa                                       >  3:359711/2‑120 (MQ=255)
                                                                                                                                         tcatcGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCtctc          >  2:303969/1‑148 (MQ=255)
                                                                                                                                            tcGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGaaa      <  1:251673/149‑1 (MQ=255)
                                                                                                                                               gTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGaaa      >  1:100664/1‑146 (MQ=255)
                                                                                                                                                tAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAgaga  <  6:248498/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTGCTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGA  >  NC_002947/543581‑543873

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