Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 160,648 (G)5→6 intergenic (‑44/+134) dksA ← / ← sfsA transcriptional regulator of rRNA transcription, DnaK suppressor protein/sugar fermentation stimulation protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913160,6431.G94.4% 302.7 / 8.3 89intergenic (‑39/+139)dksA/sfsAtranscriptional regulator of rRNA transcription, DnaK suppressor protein/sugar fermentation stimulation protein A
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base G (39/45);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

ATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTC‑GGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACT  >  NC_000913/160506‑160774
                                                                                                                                           |                                                                                                                                   
attcatATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTC‑g                                                                                                                                    >  1:1080962/1‑139 (MQ=255)
attcatATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTC‑g                                                                                                                                    >  2:253296/1‑139 (MQ=255)
          tctTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTC‑GGGGGgaaaaa                                                                                                                          >  1:957623/1‑138 (MQ=255)
            ttCGCCCGGCTTCTCCTGATATGGTTCCCCCCCAGCGATGGCGAGAATGCTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATg                                                                                                                         >  2:772986/1‑139 (MQ=255)
            ttCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATg                                                                                                                         >  2:474606/1‑139 (MQ=255)
            ttCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATg                                                                                                                         >  2:530832/1‑139 (MQ=255)
            ttCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTC‑gg                                                                                                                                   >  2:383922/1‑128 (MQ=255)
            ttCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTC‑gg                                                                                                                                   <  1:383922/128‑1 (MQ=255)
            ttCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTC‑GGGGGAAAAATGa                                                                                                                        <  1:1354892/139‑1 (MQ=255)
               gCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATGCTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      <  1:772986/139‑1 (MQ=255)
               gCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      <  2:1020331/139‑1 (MQ=255)
                  cGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCg                                                                                                                   <  1:486731/139‑1 (MQ=255)
                  cGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCg                                                                                                                   <  2:1080962/139‑1 (MQ=255)
                    gCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCt                                                                                                                 <  2:137970/139‑1 (MQ=255)
                    gCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCt                                                                                                                 >  1:1384333/1‑139 (MQ=255)
                              atatGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      <  1:700626/124‑1 (MQ=255)
                              atatGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      >  2:700626/1‑124 (MQ=255)
                              atatGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCa                                                                                                       >  1:1250592/1‑139 (MQ=255)
                                atGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGTCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGa                                                                                                     <  2:957623/139‑1 (MQ=255)
                                atGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGa                                                                                                     <  2:1116367/139‑1 (MQ=255)
                                   gTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGc                                                                                                  >  1:1595216/1‑139 (MQ=255)
                                   gTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGc                                                                                                  >  1:470889/1‑139 (MQ=255)
                                    ttCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCt                                                                                                 >  1:802329/1‑139 (MQ=255)
                                     tCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      <  2:315166/117‑1 (MQ=255)
                                     tCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      >  1:315166/1‑117 (MQ=255)
                                         ccccAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      >  1:621379/1‑113 (MQ=255)
                                         ccccAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      <  2:621379/113‑1 (MQ=255)
                                         ccccAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCg                                                                                                                   >  1:1664371/1‑116 (MQ=255)
                                         ccccAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCg                                                                                                                   <  2:1664371/116‑1 (MQ=255)
                                         ccccAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTcc                                                                                            >  1:506266/1‑139 (MQ=255)
                                         ccccAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACACACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTcc                                                                                            >  1:491471/1‑139 (MQ=255)
                                            cAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGa                                                                                         >  1:885911/1‑139 (MQ=255)
                                            cAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACACACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGa                                                                                         <  2:491471/139‑1 (MQ=255)
                                             aGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      >  2:313134/1‑109 (MQ=255)
                                             aGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      <  1:313134/109‑1 (MQ=255)
                                             aGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATg                                                                                                   >  1:1307274/1‑128 (MQ=255)
                                             aGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATg                                                                                                   <  2:1307274/128‑1 (MQ=255)
                                              gCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATa                                                                                       >  2:476802/1‑139 (MQ=255)
                                              gCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATa                                                                                       >  1:1141268/1‑139 (MQ=255)
                                                   ggCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCt                                                                                                                 >  2:783312/1‑108 (MQ=255)
                                                   ggCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCt                                                                                                                 <  1:783312/108‑1 (MQ=255)
                                                             cTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACg                                                                        <  2:802329/139‑1 (MQ=255)
                                                              tCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGt                                                                       <  1:1569311/139‑1 (MQ=255)
                                                              tCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGt                                                                       <  2:470889/139‑1 (MQ=255)
                                                                aGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTaa                                                                     <  1:419352/139‑1 (MQ=255)
                                                                aGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTaa                                                                     <  2:1595216/139‑1 (MQ=255)
                                                                aGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTaa                                                                     <  2:1141268/139‑1 (MQ=255)
                                                                  ggACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      <  2:767713/88‑1 (MQ=255)
                                                                  ggACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      >  1:767713/1‑88 (MQ=255)
                                                                  ggACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTC‑gg                                                                                                                                   >  1:676728/1‑74 (MQ=255)
                                                                  ggACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTC‑gg                                                                                                                                   <  2:676728/74‑1 (MQ=255)
                                                                   gACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAacac                                                                  >  1:337597/1‑139 (MQ=255)
                                                                     cGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAAAACtt                                                                >  1:586021/1‑139 (MQ=255)
                                                                       aTGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGa                                                              >  2:318273/1‑139 (MQ=255)
                                                                        tGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGAc                                                             >  1:1484253/1‑139 (MQ=255)
                                                                          ttttACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTTG‑GGGGGAAAAATGAGGCCGCTATAAATATCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACaa                                                           <  1:1722595/139‑1 (MQ=255)
                                                                          ttttACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACaa                                                           <  2:337597/139‑1 (MQ=255)
                                                                          ttttACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACaa                                                           <  1:253296/139‑1 (MQ=255)
                                                                            ttACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAACAATTATCTAAACGTAACACTTGACAACt                                                         >  2:1217053/1‑139 (MQ=255)
                                                                                            cttTCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGa                                                                                                     <  2:1081094/79‑1 (MQ=255)
                                                                                            cttGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGa                                                                                                     >  1:1081094/1‑79 (MQ=255)
                                                                                               gCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACg                                                                        <  2:93203/105‑1 (MQ=255)
                                                                                               gCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACg                                                                        >  1:93203/1‑105 (MQ=255)
                                                                                                     tGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCa                                <  2:1384333/139‑1 (MQ=255)
                                                                                                     tGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTC‑GGGGGAAAAATGGGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCAc                               >  1:1231782/1‑139 (MQ=255)
                                                                                                      gCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      >  2:505926/1‑52 (MQ=255)
                                                                                                      gCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      <  1:271821/52‑1 (MQ=255)
                                                                                                      gCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      <  1:505926/52‑1 (MQ=255)
                                                                                                      gCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAgg                                                                                                                      >  2:271821/1‑52 (MQ=255)
                                                                                                       cTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGAt                                                                                        <  2:1287672/81‑1 (MQ=255)
                                                                                                       cTTCTCCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCAGGAt                                                                                        >  1:1287672/1‑81 (MQ=255)
                                                                                                         tctcCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCAcgcg                            >  2:863755/1‑139 (MQ=255)
                                                                                                           tcCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGAc                          <  2:1484253/139‑1 (MQ=255)
                                                                                                           tcCTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGAc                          <  1:1713760/139‑1 (MQ=255)
                                                                                                            ccTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGAcc                         <  2:1250592/139‑1 (MQ=255)
                                                                                                             cTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCa                        <  1:474606/139‑1 (MQ=255)
                                                                                                             cTTAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCa                        <  1:318273/139‑1 (MQ=255)
                                                                                                               tAACACGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGa                      >  2:672965/1‑139 (MQ=255)
                                                                                                                 acacGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGAc                          >  2:1161577/1‑133 (MQ=255)
                                                                                                                 acacGCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGAc                          <  1:1161577/133‑1 (MQ=255)
                                                                                                                     gCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCt                                                                                                                 <  2:1329026/42‑1 (MQ=255)
                                                                                                                     gCACTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCt                                                                                                                 >  1:1329026/1‑42 (MQ=38)
                                                                                                                       aCTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGa                                                                                                     <  2:1113160/52‑1 (MQ=255)
                                                                                                                       aCTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGa                                                                                                     >  1:1113160/1‑52 (MQ=255)
                                                                                                                       aCTATCGATCCCCATG‑TTC‑GGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAAtt             >  2:1447484/1‑139 (MQ=255)
                                                                                                                        cTATCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAAAACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAAtt             <  2:586021/139‑1 (MQ=255)
                                                                                                                           tCGATCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAAc          <  1:476802/139‑1 (MQ=255)
                                                                                                                              aTCCCCATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAg       <  1:1217053/139‑1 (MQ=255)
                                                                                                                                 cccATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAATACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCAcg                              >  1:596695/1‑113 (MQ=255)
                                                                                                                                 cccATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCAcg                              <  2:596695/113‑1 (MQ=255)
                                                                                                                                   cATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACtgtgtg                                                    <  2:533268/89‑1 (MQ=255)
                                                                                                                                   cATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACtgtgtg                                                    >  1:533268/1‑89 (MQ=255)
                                                                                                                                   cATG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGttact  <  1:672965/139‑1 (MQ=255)
                                                                                                                                     tG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCtttt                                                                                              >  2:740912/1‑45 (MQ=255)
                                                                                                                                     tG‑TTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCtttt                                                                                              <  1:740912/45‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                   
ATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATG‑TTC‑GGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACT  >  NC_000913/160506‑160774

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑