Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 543,733 C→T 8.4% A657V (GCT→GTT)  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,7330CT 8.4% 196.6 / 7.8 72A657V (GCT→GTT) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base T (4/2);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01

TGCTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGAAGGT  >  NC_002947/543585‑543880
                                                                                                                                                      |                                                                                                                                                   
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                                                                                                                                         <  2:490840/132‑4 (MQ=255)
                             cGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAGCGTTACCGAtgatgc                                                                                                                                         >  1:490840/1‑129 (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:1093734/1‑150 (MQ=255)
                                                        cacCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAggg                                                                                            <  1:482611/150‑1 (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:401451/1‑150 (MQ=255)
                                                                        tCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAg                                                                            >  2:1253043/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                                                                      <  2:894820/141‑1 (MQ=255)
                                                                                       gTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATc                                                                      >  1:894820/1‑141 (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                                           <  1:1253043/150‑1 (MQ=255)
                                                                                                           gACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGaa                                           >  2:1270886/1‑150 (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)
                                                                                                                                                    cGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGAAGGt  >  2:440264/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TGCTACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTT‑‑AACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGAAGGT  >  NC_002947/543585‑543880

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: