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

ATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACT  >  NC_000913/160506‑160774
                                                                                                                                           |                                                                                                                                   
attcatATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCg                                                                                                                                    >  1:1080962/1‑139 (MQ=255)
attcatATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCg                                                                                                                                    >  2:253296/1‑139 (MQ=255)
          tctTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGgaaaaa                                                                                                                          >  1:957623/1‑138 (MQ=255)
            ttCGCCCGGCTTCTCCTGATATGGTTCCCCCCCAGCGATGGCGAGAATGCTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATg                                                                                                                         >  2:772986/1‑139 (MQ=255)
            ttCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATg                                                                                                                         >  2:474606/1‑139 (MQ=255)
            ttCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATg                                                                                                                         >  2:530832/1‑139 (MQ=255)
            ttCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCgg                                                                                                                                   >  2:383922/1‑128 (MQ=255)
            ttCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCgg                                                                                                                                   <  1:383922/128‑1 (MQ=255)
            ttCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGa                                                                                                                        <  1:1354892/139‑1 (MQ=255)
               gCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATGCTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      <  1:772986/139‑1 (MQ=255)
               gCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      <  2:1020331/139‑1 (MQ=255)
                  cGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCg                                                                                                                   <  1:486731/139‑1 (MQ=255)
                  cGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCg                                                                                                                   <  2:1080962/139‑1 (MQ=255)
                    gCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCt                                                                                                                 <  2:137970/139‑1 (MQ=255)
                    gCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCt                                                                                                                 >  1:1384333/1‑139 (MQ=255)
                              atatGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      <  1:700626/124‑1 (MQ=255)
                              atatGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      >  2:700626/1‑124 (MQ=255)
                              atatGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCa                                                                                                       >  1:1250592/1‑139 (MQ=255)
                                atGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGTCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGa                                                                                                     <  2:957623/139‑1 (MQ=255)
                                atGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGa                                                                                                     <  2:1116367/139‑1 (MQ=255)
                                   gTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGc                                                                                                  >  1:1595216/1‑139 (MQ=255)
                                   gTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGc                                                                                                  >  1:470889/1‑139 (MQ=255)
                                    ttCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCt                                                                                                 >  1:802329/1‑139 (MQ=255)
                                     tCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      <  2:315166/117‑1 (MQ=255)
                                     tCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      >  1:315166/1‑117 (MQ=255)
                                         ccccAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      >  1:621379/1‑113 (MQ=255)
                                         ccccAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      <  2:621379/113‑1 (MQ=255)
                                         ccccAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCg                                                                                                                   >  1:1664371/1‑116 (MQ=255)
                                         ccccAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCg                                                                                                                   <  2:1664371/116‑1 (MQ=255)
                                         ccccAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTcc                                                                                            >  1:506266/1‑139 (MQ=255)
                                         ccccAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACACACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTcc                                                                                            >  1:491471/1‑139 (MQ=255)
                                            cAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGa                                                                                         >  1:885911/1‑139 (MQ=255)
                                            cAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACACACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGa                                                                                         <  2:491471/139‑1 (MQ=255)
                                             aGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      >  2:313134/1‑109 (MQ=255)
                                             aGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      <  1:313134/109‑1 (MQ=255)
                                             aGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATg                                                                                                   >  1:1307274/1‑128 (MQ=255)
                                             aGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATg                                                                                                   <  2:1307274/128‑1 (MQ=255)
                                              gCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATa                                                                                       >  2:476802/1‑139 (MQ=255)
                                              gCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATa                                                                                       >  1:1141268/1‑139 (MQ=255)
                                                   ggCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCt                                                                                                                 >  2:783312/1‑108 (MQ=255)
                                                   ggCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCt                                                                                                                 <  1:783312/108‑1 (MQ=255)
                                                             cTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACg                                                                        <  2:802329/139‑1 (MQ=255)
                                                              tCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGt                                                                       <  1:1569311/139‑1 (MQ=255)
                                                              tCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGt                                                                       <  2:470889/139‑1 (MQ=255)
                                                                aGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTaa                                                                     <  1:419352/139‑1 (MQ=255)
                                                                aGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTaa                                                                     <  2:1595216/139‑1 (MQ=255)
                                                                aGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTaa                                                                     <  2:1141268/139‑1 (MQ=255)
                                                                  ggACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      <  2:767713/88‑1 (MQ=255)
                                                                  ggACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      >  1:767713/1‑88 (MQ=255)
                                                                  ggACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCgg                                                                                                                                   >  1:676728/1‑74 (MQ=255)
                                                                  ggACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCgg                                                                                                                                   <  2:676728/74‑1 (MQ=255)
                                                                   gACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAacac                                                                  >  1:337597/1‑139 (MQ=255)
                                                                     cGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAAAACtt                                                                >  1:586021/1‑139 (MQ=255)
                                                                       aTGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGa                                                              >  2:318273/1‑139 (MQ=255)
                                                                        tGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGAc                                                             >  1:1484253/1‑139 (MQ=255)
                                                                          ttttACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTTGGGGGGAAAAATGAGGCCGCTATAAATATCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACaa                                                           <  1:1722595/139‑1 (MQ=255)
                                                                          ttttACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACaa                                                           <  2:337597/139‑1 (MQ=255)
                                                                          ttttACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACaa                                                           <  1:253296/139‑1 (MQ=255)
                                                                            ttACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAACAATTATCTAAACGTAACACTTGACAACt                                                         >  2:1217053/1‑139 (MQ=255)
                                                                                            cttTCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGa                                                                                                     <  2:1081094/79‑1 (MQ=255)
                                                                                            cttGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGa                                                                                                     >  1:1081094/1‑79 (MQ=255)
                                                                                               gCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACg                                                                        <  2:93203/105‑1 (MQ=255)
                                                                                               gCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACg                                                                        >  1:93203/1‑105 (MQ=255)
                                                                                                     tGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCa                                <  2:1384333/139‑1 (MQ=255)
                                                                                                     tGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGGGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCAc                               >  1:1231782/1‑139 (MQ=255)
                                                                                                      gCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      >  2:505926/1‑52 (MQ=255)
                                                                                                      gCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      <  1:271821/52‑1 (MQ=255)
                                                                                                      gCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      <  1:505926/52‑1 (MQ=255)
                                                                                                      gCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAgg                                                                                                                      >  2:271821/1‑52 (MQ=255)
                                                                                                       cTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGAt                                                                                        <  2:1287672/81‑1 (MQ=255)
                                                                                                       cTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCAGGAt                                                                                        >  1:1287672/1‑81 (MQ=255)
                                                                                                         tctcCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCAcgcg                            >  2:863755/1‑139 (MQ=255)
                                                                                                           tcCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGAc                          <  2:1484253/139‑1 (MQ=255)
                                                                                                           tcCTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGAc                          <  1:1713760/139‑1 (MQ=255)
                                                                                                            ccTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGAcc                         <  2:1250592/139‑1 (MQ=255)
                                                                                                             cTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCa                        <  1:474606/139‑1 (MQ=255)
                                                                                                             cTTAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCa                        <  1:318273/139‑1 (MQ=255)
                                                                                                               tAACACGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGa                      >  2:672965/1‑139 (MQ=255)
                                                                                                                 acacGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGAc                          >  2:1161577/1‑133 (MQ=255)
                                                                                                                 acacGCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGAc                          <  1:1161577/133‑1 (MQ=255)
                                                                                                                     gCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCt                                                                                                                 <  2:1329026/42‑1 (MQ=255)
                                                                                                                     gCACTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCt                                                                                                                 >  1:1329026/1‑42 (MQ=38)
                                                                                                                       aCTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGa                                                                                                     <  2:1113160/52‑1 (MQ=255)
                                                                                                                       aCTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGa                                                                                                     >  1:1113160/1‑52 (MQ=255)
                                                                                                                       aCTATCGATCCCCATGTTCGGGGGGAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAAtt             >  2:1447484/1‑139 (MQ=255)
                                                                                                                        cTATCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAAAACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAAtt             <  2:586021/139‑1 (MQ=255)
                                                                                                                           tCGATCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAAc          <  1:476802/139‑1 (MQ=255)
                                                                                                                              aTCCCCATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAg       <  1:1217053/139‑1 (MQ=255)
                                                                                                                                 cccATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAATACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCAcg                              >  1:596695/1‑113 (MQ=255)
                                                                                                                                 cccATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCAcg                              <  2:596695/113‑1 (MQ=255)
                                                                                                                                   cATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACtgtgtg                                                    <  2:533268/89‑1 (MQ=255)
                                                                                                                                   cATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACtgtgtg                                                    >  1:533268/1‑89 (MQ=255)
                                                                                                                                   cATGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGttact  <  1:672965/139‑1 (MQ=255)
                                                                                                                                     tGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCtttt                                                                                              >  2:740912/1‑45 (MQ=255)
                                                                                                                                     tGTTCGGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCtttt                                                                                              <  1:740912/45‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                   
ATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACT  >  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: