Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,020,512 A→C 5.7% intergenic (+23/‑28) ssnA → / → ygfM putative chlorohydrolase/aminohydrolase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,020,5120AC 5.7% 424.4 / 15.6 124intergenic (+23/‑28)ssnA/ygfMputative chlorohydrolase/aminohydrolase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (63/52);  new base C (2/5);  total (65/59)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

Only 100 of 127 total aligned reads displayed.

AGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCGTA  >  NC_000913/3020364‑3020660
                                                                                                                                                     |                                                                                                                                                    
  cGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAgg                                                                                                                                                    <  3:113654/149‑1 (MQ=255)
  cGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAgg                                                                                                                                                    <  2:537120/149‑1 (MQ=255)
    tGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGa                                                                                                                                                  >  1:315679/1‑149 (MQ=255)
        ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGC‑GGGAagag                                                                                                                                              <  4:296224/149‑1 (MQ=255)
        ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAagag                                                                                                                                              <  2:55275/149‑1 (MQ=255)
        ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAagag                                                                                                                                              <  1:640781/149‑1 (MQ=255)
         gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCCTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAAGAgg                                                                                                                                             <  1:162956/149‑1 (MQ=255)
         gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAAGAgg                                                                                                                                             <  3:575130/149‑1 (MQ=255)
         gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAAGAgg                                                                                                                                             <  3:64382/149‑1 (MQ=255)
          tCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAggg                                                                                                                                            <  1:579437/149‑1 (MQ=255)
          tCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAggg                                                                                                                                            <  4:451800/149‑1 (MQ=255)
                gtgtGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAAt                                                                                                                                      >  3:108761/1‑149 (MQ=255)
                  gtgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATaa                                                                                                                                    >  2:55891/1‑149 (MQ=255)
                   tgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAg                                                                                                                                   <  4:425984/149‑1 (MQ=255)
                     aTGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGa                                                                                                                                 >  4:171660/1‑149 (MQ=255)
                       ggTCTATGAAGACCGTCAGTTTAACTTCGTTTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTGAAATAAGGaaa                                                                                                                               >  2:255063/1‑149 (MQ=255)
                        gTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAaag                                                                                                                              >  4:235669/1‑149 (MQ=255)
                                aGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgatgat                                                                                                                      >  2:227981/1‑149 (MQ=255)
                                 gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                     <  1:255063/149‑1 (MQ=255)
                                 gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                     <  4:35051/149‑1 (MQ=255)
                                 gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                     <  4:108761/149‑1 (MQ=255)
                                    cGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGaa                                                                                                                  <  2:403163/149‑1 (MQ=255)
                                    cGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGaa                                                                                                                  <  2:891699/149‑1 (MQ=255)
                                     gTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAagaaga                                                                                                                          <  3:625596/140‑1 (MQ=255)
                                     gTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAagaaga                                                                                                                          >  4:625596/1‑140 (MQ=255)
                                     gTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgatgat                                                                                                                      <  1:519640/144‑1 (MQ=255)
                                     gTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAAc                                                                                                                 <  1:407876/149‑1 (MQ=255)
                                          tttAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTaaa                                                                                                                                       >  3:63757/1‑122 (MQ=255)
                                             aaCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAaga                                                                                                                                               <  4:485334/111‑1 (MQ=255)
                                             aaCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAaga                                                                                                                                               >  3:485334/1‑111 (MQ=255)
                                              aCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCa                                                                                                        <  3:309860/149‑1 (MQ=255)
                                               cTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGtt                                                                                                                                          >  2:248360/1‑114 (MQ=255)
                                               cTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGtt                                                                                                                                          <  1:248360/114‑1 (MQ=255)
                                               cTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAg                                                                                                       <  3:235669/149‑1 (MQ=255)
                                                ttCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAgg                                                                                                      >  4:511589/1‑149 (MQ=255)
                                                ttCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAgg                                                                                                      >  4:648601/1‑149 (MQ=255)
                                                ttCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAgg                                                                                                      >  4:757561/1‑149 (MQ=255)
                                                ttCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAgg                                                                                                      <  1:163599/149‑1 (MQ=255)
                                                   gATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGccc                                                                                                   >  1:702794/1‑149 (MQ=255)
                                                   gATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGccc                                                                                                   >  4:192713/1‑149 (MQ=255)
                                                    aTTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCg                                                                                                  >  3:581250/1‑149 (MQ=255)
                                                    aTTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCg                                                                                                  >  4:702065/1‑149 (MQ=255)
                                                    aTTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCg                                                                                                  >  3:812696/1‑149 (MQ=255)
                                                           ttCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGt                                                                                           <  2:199658/149‑1 (MQ=255)
                                                            tCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTc                                                                                          <  4:676635/149‑1 (MQ=255)
                                                               aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGaa                                                                                       >  1:113047/1‑149 (MQ=255)
                                                               aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGaa                                                                                       >  4:19090/1‑149 (MQ=255)
                                                               aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGaa                                                                                       <  2:554944/149‑1 (MQ=255)
                                                                   aTGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAagag                                                                                                                                              >  3:601603/1‑90 (MQ=255)
                                                                   aTGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAagag                                                                                                                                              <  4:601603/90‑1 (MQ=255)
                                                                   aTGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAgg                                                                                   >  3:286131/1‑149 (MQ=255)
                                                                    tGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGgc                                                                                  <  1:716996/149‑1 (MQ=255)
                                                                        caAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTg                                                                              >  3:853303/1‑149 (MQ=255)
                                                                           gCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGaa                                                                           >  1:786388/1‑149 (MQ=255)
                                                                           gCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGaa                                                                           >  2:88820/1‑149 (MQ=255)
                                                                            ccAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAAc                                                                          >  4:33574/1‑149 (MQ=255)
                                                                             cAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACt                                                                         >  3:359635/1‑149 (MQ=255)
                                                                             cAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACt                                                                         <  1:376638/149‑1 (MQ=255)
                                                                              aGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACa                                                                                     <  2:298556/136‑1 (MQ=255)
                                                                              aGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTg                                                                        >  2:440207/1‑149 (MQ=255)
                                                                                aaaaGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCA‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCGGGAACTGaa                                                                      >  3:480184/1‑149 (MQ=255)
                                                                                  aaGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAgc                                                                    >  3:757037/1‑149 (MQ=255)
                                                                                    gCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAgcgc                                                                  >  1:829033/1‑149 (MQ=255)
                                                                                            cAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGAc                                                                                                >  2:387156/1‑111 (MQ=255)
                                                                                            cAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCg                                                                                         >  4:827456/1‑118 (MQ=255)
                                                                                            cAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCg                                                                                         <  3:827456/118‑1 (MQ=255)
                                                                                            cAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACCGGCGCTGGAACTGAAGCGCCGCTACCa                                                          <  4:313239/149‑1 (MQ=255)
                                                                                             aGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAg                                                         >  3:410387/1‑149 (MQ=255)
                                                                                                    ggCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAg                                                  <  4:451131/149‑1 (MQ=255)
                                                                                                     gCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                     >  1:119749/1‑81 (MQ=255)
                                                                                                     gCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                     <  2:119749/81‑1 (MQ=255)
                                                                                                        tCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGt                                              >  1:365592/1‑149 (MQ=255)
                                                                                                         cGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTc                                             <  4:581250/149‑1 (MQ=255)
                                                                                                         cGGATGGATGCGCTGGCATAAATGACAGATGCCATCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTc                                             <  4:812696/149‑1 (MQ=255)
                                                                                                              ggatgCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCa                                                                                                        <  1:351603/85‑1 (MQ=255)
                                                                                                                atgCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCg                                      >  1:435536/1‑149 (MQ=255)
                                                                                                                      tGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCgggg                                <  1:227981/149‑1 (MQ=255)
                                                                                                                        gCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgc                              >  4:548777/1‑149 (MQ=255)
                                                                                                                          aTAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGaa                                                   <  1:391079/126‑1 (MQ=255)
                                                                                                                          aTAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGaa                                                   >  2:391079/1‑126 (MQ=255)
                                                                                                                          aTAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgcag                            >  4:54827/1‑149 (MQ=255)
                                                                                                                          aTAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgcag                            >  4:521297/1‑149 (MQ=255)
                                                                                                                            aaaTGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgcagca                          >  1:530648/1‑149 (MQ=255)
                                                                                                                                gACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACt                      >  4:578722/1‑149 (MQ=255)
                                                                                                                                 aCAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTc                     >  3:499760/1‑149 (MQ=255)
                                                                                                                                    gATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAAt                  <  2:568883/149‑1 (MQ=255)
                                                                                                                                     aTGGCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAGCAGTTTTTCAGGCGCGACTCTGTCGAACAGGCGCTGGAACTGAGGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGGCGGGGGCAGCAAACTCAATg                 <  2:204957/149‑1 (MQ=255)
                                                                                                                                      tGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGc                >  4:545687/1‑149 (MQ=255)
                                                                                                                                        cccTCTTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAAc                                                                          <  3:336686/89‑1 (MQ=39)
                                                                                                                                        cccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTa              <  3:33574/149‑1 (MQ=255)
                                                                                                                                         ccTCTTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAGACTCAATGCTac             >  1:580612/1‑149 (MQ=255)
                                                                                                                                         ccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac             >  1:498211/1‑149 (MQ=255)
                                                                                                                                         ccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac             >  1:788006/1‑149 (MQ=255)
                                                                                                                                         ccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac             >  4:907929/1‑149 (MQ=255)
                                                                                                                                           tctTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTacac           >  4:459314/1‑149 (MQ=255)
                                                                                                                                           tctTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTacac           >  1:462937/1‑149 (MQ=255)
                                                                                                                                             ttCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCa         >  1:290137/1‑149 (MQ=255)
                                                                                                                                                   gCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCg    >  2:591423/1‑149 (MQ=255)
                                                                                                                                                     aGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCGTa  <  2:579015/149‑1 (MQ=255)
                                                                                                                                                     aGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCGTa  <  3:545687/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCGTA  >  NC_000913/3020364‑3020660

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