Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 543,730 A→G 8.5% 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,7300AG 8.5% 193.9 / 6.1 70N656S (AAC→AGC) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base G (4/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CCAGTGCTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGAAGG  >  NC_002947/543581‑543879
                                                                                                                                                       |                                                                                                                                                     
ccAGTGCTACCGCGTGGTTGGTCTGAAAGAGATCGTTATTTTCGCAGACCAACTGATATACGCCGGTTTCGCATACTCGACCATTTCCGGCGATTCCATAAGTGTT‑‑CACGACTTCGTTATCCCGCACGAGAAAGCCCGCATCATCGGTaa                                                                                                                                                       <  2:1020306/150‑1 (MQ=255)
ccAGTGCTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTCCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTaa                                                                                                                                                       >  1:121058/1‑150 (MQ=255)
  aGTGCTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACg                                                                                                                                                     >  1:315470/1‑150 (MQ=255)
    tGCTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCt                                                                                                                                                   >  2:896239/1‑150 (MQ=255)
          cgcgTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAt                                                                                                                                             <  2:1233396/150‑1 (MQ=255)
                 ttggtCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaagtga                                                                                                                                      <  1:896239/150‑1 (MQ=255)
                  tggtCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCTGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaagtgaa                                                                                                                                     <  1:740660/150‑1 (MQ=255)
                            gagaCCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGGCGACAAAGCCCGGAGCATCGGTAACGCTACCGATGAAGTGAAGGAAATCgag                                                                                                                           <  1:739890/150‑1 (MQ=255)
                             agaCCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑TTTAACTTCGTGATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCgaga                                                                                                                          >  1:1241338/1‑150 (MQ=255)
                                 cGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAGCGTTACCGAtgatgc                                                                                                                                        >  1:490840/1‑129 (MQ=255)
                                 cGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAGCGTTACCGAtgatgc                                                                                                                                        <  2:490840/132‑4 (MQ=255)
                                 cGTTATCTTCGCCGACCAAATGATGTCAAGCGGTTTCGGATAGTCGACCACTTTCGTCTTTTCCCTCGGGGGTCCAAAGGCCAAGGACTCATTGAGGAG‑‑AGCCCGCATCATCGGTAACGCTACCGATGAAGTAGTGGAAGGCGAGAGCCa                                                                                                                      >  1:678944/1‑150 (MQ=255)
                                  gTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAg                                                                                                                     <  1:29167/150‑1 (MQ=255)
                                     aTCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaagt                                                                                                                                        <  2:144686/128‑1 (MQ=255)
                                     aTCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaagt                                                                                                                                        >  1:144686/1‑128 (MQ=255)
                                     aTCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTAc                                                                                                                  <  1:896634/150‑1 (MQ=255)
                                        ttCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGcc                                                                                                               >  1:559369/1‑150 (MQ=255)
                                                 cAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaa                                                                                                                                          >  1:880949/1‑114 (MQ=255)
                                                 cAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaa                                                                                                                                          <  2:880949/114‑1 (MQ=255)
                                                   aCTGATGTACACCGGTTTCGCATACTCGACCATTTCCTGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtga                                                                                                                                           >  1:822612/1‑111 (MQ=255)
                                                   aCTGATGTACACCGGTTTCGCATACTCGACCATTTCCTGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtga                                                                                                                                           <  2:822612/111‑1 (MQ=255)
                                                   aCTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGt                                                                                                    >  1:271256/1‑150 (MQ=255)
                                                     tGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTaa                                                                                                  >  1:938153/1‑150 (MQ=255)
                                                        tGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTg                                                                                                      >  1:204129/1‑143 (MQ=255)
                                                        tGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTg                                                                                                      <  2:204129/143‑1 (MQ=255)
                                                           acacCGGTTTCGCATACTCGACAATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAgg                                                                                            >  1:1376015/1‑150 (MQ=255)
                                                            cacCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAggg                                                                                           <  1:482611/150‑1 (MQ=255)
                                                            cacCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAggg                                                                                           >  1:1093734/1‑150 (MQ=255)
                                                             acCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCa                                                                                                                      >  2:950370/1‑122 (MQ=255)
                                                             acCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCa                                                                                                                      <  1:950370/122‑1 (MQ=255)
                                                             acCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGc                                                                                          >  2:1291891/1‑150 (MQ=255)
                                                             acCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGGT‑‑AGCGACTTCGTTACCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTTCGCCTCCGGCCTGGTAACCCAGGGc                                                                                          >  2:1342416/1‑150 (MQ=255)
                                                               cGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCga                                                                                        <  1:1395165/150‑1 (MQ=255)
                                                                  tttCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAGCGTTACCGAtgatgcgg                                                                                                                                      >  2:967856/1‑96 (MQ=39)
                                                                  tttCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAGCGTTACCGAtgatgcgg                                                                                                                                      <  1:967856/101‑6 (MQ=39)
                                                                   ttCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAg                                                                                    <  1:929562/150‑1 (MQ=255)
                                                                      gCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTac                                                                                 <  2:75357/150‑1 (MQ=255)
                                                                      gCATACTCGACCATTTCCGGCGTTTCAATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGACAAGTac                                                                                 >  2:386872/1‑150 (MQ=255)
                                                                            tCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCgag                                                                                       >  2:717829/1‑138 (MQ=255)
                                                                            tCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCgag                                                                                       <  1:717829/138‑1 (MQ=255)
                                                                            tCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAg                                                                           >  2:1253043/1‑150 (MQ=255)
                                                                            tCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAg                                                                           >  2:401451/1‑150 (MQ=255)
                                                                              gACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGt                                                                         >  2:90655/1‑150 (MQ=255)
                                                                              gACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGt                                                                         <  2:559369/150‑1 (MQ=255)
                                                                                  aTTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATc                                                                     <  1:1291891/150‑1 (MQ=255)
                                                                                   tttCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACATCAAGAGCATCg                                                                    >  2:443101/1‑150 (MQ=255)
                                                                                     tCCGGCGTTTCCATCGGGGTT‑‑AACGGCTTCGGTATCCCGGAGGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGAc                                                                  <  2:678944/150‑1 (MQ=255)
                                                                                       cGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACtt                                                                <  2:315470/150‑1 (MQ=255)
                                                                                           gTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATc                                                                     >  1:894820/1‑141 (MQ=255)
                                                                                           gTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATc                                                                     <  2:894820/141‑1 (MQ=255)
                                                                                                  tCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCt                                                                                                              <  2:781976/93‑1 (MQ=255)
                                                                                                  tCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCt                                                                                                              >  1:781976/1‑93 (MQ=255)
                                                                                                      tgtT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCa                                                                                              <  1:1190185/105‑1 (MQ=255)
                                                                                                      tgtT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCa                                                                                              >  2:1190185/1‑105 (MQ=255)
                                                                                                        tt‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAacg                                               >  1:111212/1‑150 (MQ=255)
                                                                                                             aCGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAac                                                >  2:12633/1‑146 (MQ=255)
                                                                                                             aCGACTTCGTTATCCCGGACGAGAAAGCCCGCATAATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAac                                                <  1:12633/146‑1 (MQ=255)
                                                                                                               gACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaagtgaa                                                                                                                                     >  2:93340/1‑59 (MQ=255)
                                                                                                               gACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtgaagtgaa                                                                                                                                     <  1:93340/59‑1 (MQ=255)
                                                                                                               gACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGaa                                          >  2:1270886/1‑150 (MQ=255)
                                                                                                               gACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGaa                                          <  1:1253043/150‑1 (MQ=255)
                                                                                                                aCTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATc                                                                     >  2:564732/1‑122 (MQ=255)
                                                                                                                aCTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATc                                                                     <  1:564732/122‑1 (MQ=255)
                                                                                                                 cTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAgt                                        >  2:1349133/1‑150 (MQ=255)
                                                                                                                  ttCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAgtg                                       >  1:1147151/1‑150 (MQ=255)
                                                                                                                        aTCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGGCCAAg                                 <  1:1270886/150‑1 (MQ=255)
                                                                                                                                   gaAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGc                      >  1:1292527/1‑150 (MQ=255)
                                                                                                                                    acaGCCCGCATCTTCAGCAACGATATCGCTGAAGGTTAAGAAATCGATAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGcc                     <  2:1016782/148‑1 (MQ=255)
                                                                                                                                         ccGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACct                >  1:1301345/1‑150 (MQ=255)
                                                                                                                                            catcatCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCcctc             <  2:1241338/150‑1 (MQ=255)
                                                                                                                                             ttcatcGGTAGCGTTACCGATGAAGTGAAGGAAATCGAGAGCCAGTAAGGCTAAGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCg            >  2:1416495/2‑150 (MQ=255)
                                                                                                                                                atcGGTAGCGTTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGaaa         >  2:12309/1‑150 (MQ=255)
                                                                                                                                                      aaCGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGAAg   <  2:1292527/150‑1 (MQ=255)
                                                                                                                                                       aCGCTACCGATGAAGTGAACGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGAAgg  <  2:938153/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CCAGTGCTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGAAGG  >  NC_002947/543581‑543879

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