Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,886,180 G→C 7.7% intergenic (+14/+109) PP_4295 → / ← PP_4296 TetR family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,886,1800GC 7.7% 224.3 / 7.9 74intergenic (+14/+109)PP_4295/PP_4296TetR family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/40);  new base C (3/3);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

ACCTATGCCGACTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACC  >  NC_002947/4886035‑4886301
                                                                                                                                                 |                                                                                                                         
aCCTATGCCGACTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGTCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCAT                                                                                                                       >  2:816763‑M1/1‑150 (MQ=255)
aCCTATGCCGACTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCAT                                                                                                                       >  2:112771‑M1/1‑150 (MQ=255)
 ccTATGCCGACTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATT                                                                                                                      <  1:1502892‑M1/150‑1 (MQ=255)
  cTATGCCGACTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTC                                                                                                                     <  2:68181‑M1/150‑1 (MQ=255)
      gCCGACTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGG                                                                                                                 <  1:226379‑M1/150‑1 (MQ=255)
      gCCGACTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGG                                                                                                                 <  1:1444311‑M1/150‑1 (MQ=255)
      gCCGACTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGG                                                                                                                 >  1:2418389‑M1/1‑150 (MQ=255)
        cGACTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAggg                                                                                                                           <  1:2036933/138‑1 (MQ=255)
        cGACTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAggg                                                                                                                           >  2:2036933/1‑138 (MQ=255)
                aCTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGC                                                                                                       <  1:974540‑M1/150‑1 (MQ=255)
                    gCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCC                                                                                                   >  2:1733852‑M1/1‑150 (MQ=255)
                    gCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGC                                                                                                           <  1:374650‑M1/142‑1 (MQ=255)
                    gCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGC                                                                                                           >  2:374650‑M1/1‑142 (MQ=255)
                     cAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCTTGCCCGCgaatg                                                                                                  <  1:2013145‑M1/150‑6 (MQ=255)
                      aGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCTTGCCCGCgaatgg                                                                                                 <  1:262301‑M1/150‑7 (MQ=255)
                       gATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACA                                                                                                <  1:2175414‑M1/150‑1 (MQ=255)
                       gATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACA                                                                                                <  2:1994761‑M1/150‑1 (MQ=255)
                          ttCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCTCGGGCAAGCCCGCTCCCACAGGa                                                                                             <  1:793163/150‑1 (MQ=255)
                              gCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATcac                                                                                         <  1:165889/150‑1 (MQ=255)
                                ggTGACCGGTAAGGCCAAGCTGGCAGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATcacca                                                                                       <  2:506688/150‑1 (MQ=255)
                                  tGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACcaca                                                                                     <  1:777935/150‑1 (MQ=255)
                                      cgatctggCCAAGCTGGCCGACAGTGATTATGACGCGGGATGCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGcc                                                                                 <  1:1204696/144‑1 (MQ=255)
                                          aaGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGCACCGTCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGgactccacaaccttcag                                                                             <  1:426336‑M1/150‑18 (MQ=255)
                                          aaGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGaa                                                                             <  1:898903/150‑1 (MQ=255)
                                           atgCCAGGCTGGCCGACAGTGACTCTGCCGTGGTGGCCGACTGTATCATTGGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCAGCTCCCACAGGATCACCACATGCCTGAAg                                                                            <  2:624483/148‑1 (MQ=255)
                                             gCCAAGCTGGCCGACAGTGATTATGATGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAAACATCTCGAAGGCCCAATGCCGCGATAGCCCACTTCCACAAGAGCACCACATGCCtagatct                                                                          >  2:1204696/1‑144 (MQ=255)
                                             gCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGt                                                                          >  1:1930940/1‑150 (MQ=255)
                                             gCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGt                                                                          >  1:81337/1‑150 (MQ=255)
                                                 aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                      <  1:112771/150‑1 (MQ=255)
                                                 aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                      <  1:1733852/150‑1 (MQ=255)
                                                   cgggCCGACTGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                    <  2:561387/148‑1 (MQ=255)
                                                     gtccgaCAGTGATTATGCCGCGGCGGCAGCGACCATCATTCGCATGGCCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGatcaccacaaccttcaggcatgtggtga                                                                  <  2:501959‑M1/148‑29 (MQ=255)
                                                     ggCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGatcaccacaaccttcaggcatgtggtga                                                                  <  2:499542‑M1/150‑29 (MQ=255)
                                                      gCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGCGGCAGCCTACAACCTCCAGGAGGGCAATTCAAACGCGAGTCAACTGCCGCTACATCACCACTTGCCTGAAGGGTGTGGAGTc                                                                 <  2:1843357/150‑1 (MQ=255)
                                                         gACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcct                                                              >  1:1279970‑M2/1‑147 (MQ=255)
                                                           cAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGCACTATCTCGAAGGCCCATTCGCGGGCAAGCCCGCTCCCACAGGgactccccaaccttcaggcatgtggtgatcctgt                                                            <  1:2236659‑M1/150‑35 (MQ=255)
                                                           cAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgt                                                            >  1:919665‑M2/1‑145 (MQ=255)
                                                           cAGTGATTATGACGCGGCGGCCGAGACAATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGTGCCATTCGCGGGCAGGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgt                                                            >  2:439380‑M2/1‑145 (MQ=255)
                                                            aGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgtg                                                           <  1:1178112‑M2/150‑7 (MQ=255)
                                                             gTGATTATGACGCGGCGGCCGAGATCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGTAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGgactccacagccttcaggcatgtggtgatcctgtgg                                                          <  2:1209091‑M1/150‑37 (MQ=255)
                                                             gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAATGTTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCAATTCGCGGGCACGCCCGAACCCACAGGATCACCAGGTTCTTGAAGGTTGTTGAGTCcctgtgg                                                          >  2:1844130‑M2/1‑143 (MQ=255)
                                                             gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTTTGGAGTCcctgtgg                                                          >  2:81337‑M2/1‑143 (MQ=255)
                                                              tGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgtgg                                                          <  2:1835884‑M2/149‑8 (MQ=255)
                                                              tGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgtgg                                                          >  1:1835884‑M2/1‑142 (MQ=255)
                                                               gATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgtggga                                                        >  2:1680016‑M2/1‑141 (MQ=255)
                                                                   aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGgactccacaaccttcaggcatgtggtgatcctgtgggagcgg                                                    <  2:341897‑M1/150‑43 (MQ=255)
                                                                   aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgtgggagcgg                                                    >  2:2418389‑M2/1‑137 (MQ=255)
                                                                     gACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAGATGTGGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgtgggagcgggc                                                  >  2:564975‑M2/1‑135 (MQ=255)
                                                                      aCGCGGCGGCCGAGACCATCATTCGTATGGTATTGCAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgtgggagcgggct                                                 <  1:564975‑M2/150‑17 (MQ=255)
                                                                         cggcggCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgtgggagcgggcttgc                                              >  2:1178112‑M2/1‑131 (MQ=255)
                                                                          ggcggcCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgtgggagcgggcttgcc                                             >  2:1930940‑M2/1‑130 (MQ=255)
                                                                             ggcCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgtgggagcgggctt                                                <  1:2401468‑M2/144‑18 (MQ=255)
                                                                             ggcCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgtgggagcgggctt                                                >  2:2401468‑M2/1‑127 (MQ=255)
                                                                                cGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGCCCATTCGCGGGCAAGCCCGCTCCCACAGGgactccacaaccttcaggcatgtggtgatcctgt                                                            <  1:1397678‑M1/129‑35 (MQ=255)
                                                                                cGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGCCCATTCGCGGGCAAGCCCGCTCCCACAGGgactccacaaccttcaggcatgtggtgatcctgt                                                            >  2:1397678‑M1/1‑95 (MQ=255)
                                                                                 gagaCCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgtgggagcgggcttgcccgcgaat                                      <  2:74583‑M2/150‑28 (MQ=255)
                                                                                       atcatTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACA                                                                                                >  2:2131464‑M1/1‑86 (MQ=255)
                                                                                       atcatTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACA                                                                                                <  1:2131464‑M1/86‑1 (MQ=255)
                                                                                               tATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgtgg                                                          >  2:197908‑M2/1‑109 (MQ=255)
                                                                                               tATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCcctgtgg                                                          <  1:197908‑M2/116‑8 (MQ=255)
                                                                                               tATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGCCCATTCGCGGGCAAGCCCGCTCCCACAGGgactccacaaccttcaggcatgtggtgatcctgtgggagc                                                      >  2:852344‑M1/1‑80 (MQ=255)
                                                                                               tATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGCCCATTCGCGGGCAAGCCCGCTCCCACAGGgactccacaaccttcaggcatgtggtgatcctgtgggagc                                                      <  1:852344‑M1/120‑41 (MQ=255)
                                                                                                        gAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGaa                                                                             >  1:2228126/1‑88 (MQ=255)
                                                                                                        gAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGaa                                                                             <  2:2228126/88‑1 (MQ=255)
                                                                                                          aaGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGgactccacaaccttcaggcatgtggtgatcctgtg                                                           >  1:1877974‑M1/1‑69 (MQ=255)
                                                                                                          aaGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGgactccacaaccttcaggcatgtggtgatcctgtg                                                           <  2:1877974‑M1/104‑36 (MQ=255)
                                                                                                              ttGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGc         <  1:74583/150‑1 (MQ=255)
                                                                                                                  gAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                    >  2:745339/1‑87 (MQ=255)
                                                                                                                  gAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                    <  1:745339/87‑1 (MQ=255)
                                                                                                                     ccTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTAcc  >  1:1759546/1‑150 (MQ=255)
                                                                                                                             aGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                       >  1:1968372/1‑73 (MQ=255)
                                                                                                                             aGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                       <  2:1968372/73‑1 (MQ=255)
                                                                                                                                  gAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtg                                                                        <  2:478534/67‑1 (MQ=255)
                                                                                                                                  gAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtg                                                                        >  1:478534/1‑67 (MQ=255)
                                                                                                                                     ccATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAG                                                                                               >  1:22802‑M1/1‑41 (MQ=255)
                                                                                                                                     ccATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAG                                                                                               <  2:22802‑M1/41‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                         
ACCTATGCCGACTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACC  >  NC_002947/4886035‑4886301

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