Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,020,511 C→G 8.7% intergenic (+22/‑29) ssnA → / → ygfM putative aminohydrolase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,020,5110CG 8.7% 287.9 / 13.6 91intergenic (+22/‑29)ssnA/ygfMputative aminohydrolase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (36/46);  new base G (4/4);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

CAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCG  >  NC_000913/3020363‑3020658
                                                                                                                                                    |                                                                                                                                                   
caGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCa                                                                                                                                                    >  1:663927/1‑150 (MQ=255)
 aGCGGGATGGTCAATGGTGTGACTGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACCAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAg                                                                                                                                                   <  2:1336692/150‑1 (MQ=255)
  gCGGGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTACGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTAGGTGGCGTCGGATGGAAGCGCTGGCATAAATGACAGATGCCCTCTACCCGCAgg                                                                                                                                                  >  2:1613407/1‑150 (MQ=255)
   cGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAggg                                                                                                                                                 <  2:534707/150‑1 (MQ=255)
   cGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAggg                                                                                                                                                 <  2:165792/150‑1 (MQ=255)
   cGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAggg                                                                                                                                                 <  2:1446267/150‑1 (MQ=255)
      gATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAag                                                                                                                                              <  2:1299722/150‑1 (MQ=255)
        tGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAagag                                                                                                                                            <  1:310538/150‑1 (MQ=255)
          gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAggg                                                                                                                                          <  2:417724/150‑1 (MQ=255)
           tCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAgggc                                                                                                                                         <  2:274689/150‑2 (MQ=255)
               tGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAgggcatc                                                                                                                                      >  2:1298965/1‑145 (MQ=255)
               tGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAgggcatc                                                                                                                                      <  1:1298965/149‑5 (MQ=255)
                ggggtgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAACCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAAt                                                                                                                                    <  1:184103/147‑1 (MQ=255)
                ggTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAAt                                                                                                                                    <  1:1061324/150‑1 (MQ=255)
                ggTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAACCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAAt                                                                                                                                    <  1:178248/150‑1 (MQ=255)
                 gtgtGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATa                                                                                                                                   <  1:192491/150‑1 (MQ=255)
                   gtgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAATCCATAAAATCCGCTGCCATTATGTGGCGTCTGTTGGATGCGCTGTCATAAATGACAGATGCCATCTTACCGAAGGGAATATGGTTAAATAAg                                                                                                                                 <  2:166766/150‑1 (MQ=255)
                   gtgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAg                                                                                                                                 <  1:846918/150‑1 (MQ=255)
                    tgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAgg                                                                                                                                >  1:1042208/1‑150 (MQ=255)
                     gATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGa                                                                                                                               >  1:635723/1‑150 (MQ=255)
                      aTGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGaa                                                                                                                              <  2:611594/150‑1 (MQ=255)
                      aTGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGaa                                                                                                                              <  1:31354/150‑1 (MQ=255)
                        ggTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAaag                                                                                                                            >  2:1448513/1‑150 (MQ=255)
                             attaataCCGTCAGTTTAACTTCGATTGCGATTGCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGTCGTCGGATGGATTCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGTGTTAAATAAGGAAAGAAgat                                                                                                                       <  1:794731/144‑1 (MQ=255)
                              tGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgatg                                                                                                                      >  2:1094062/1‑150 (MQ=255)
                               gAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgatga                                                                                                                     >  2:628414/1‑150 (MQ=255)
                               gAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgatga                                                                                                                     >  1:1576211/1‑150 (MQ=255)
                                 aGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                   <  2:1042208/150‑1 (MQ=255)
                                  gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAATCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTg                                                                                                                  <  2:1628262/150‑1 (MQ=255)
                                  gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTg                                                                                                                  >  2:893521/1‑150 (MQ=255)
                                  gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTg                                                                                                                  >  2:887675/1‑150 (MQ=255)
                                  gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTg                                                                                                                  <  1:717290/150‑1 (MQ=255)
                                    ccGTCAGTTTAACTTCAATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGaa                                                                                                                >  1:735900/1‑150 (MQ=255)
                                      gTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACa                                                                                                              >  1:362514/1‑150 (MQ=255)
                                      gTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACa                                                                                                              >  2:1273627/1‑150 (MQ=255)
                                         aGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAtt                                                                                                           >  2:1247496/1‑150 (MQ=255)
                                          gTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttt                                                                                                          >  1:1483742/1‑150 (MQ=255)
                                            ttAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttttt                                                                                                        <  1:1084330/150‑1 (MQ=255)
                                             tAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTc                                                                                                       <  2:635723/150‑1 (MQ=255)
                                              aaCTTCGATTGCGATTCCATTTTTGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAg                                                                                                                                 >  1:867279/1‑123 (MQ=255)
                                              aaCTTCGATTGCGATTCCATTTTTGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAg                                                                                                                                 <  2:867279/123‑1 (MQ=255)
                                                cTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAgg                                                                                                    >  1:1417957/1‑150 (MQ=255)
                                                   cGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGccc                                                                                                 >  1:721174/1‑150 (MQ=255)
                                                        gCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGActc                                                                                            <  1:990070/150‑1 (MQ=255)
                                                            ttCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTc                                                                                        <  1:39222/150‑1 (MQ=255)
                                                            ttCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTc                                                                                        <  2:712516/150‑1 (MQ=255)
                                                             tCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCg                                                                                       <  1:300468/150‑1 (MQ=255)
                                                              ccATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGa                                                                                      <  2:1151190/150‑1 (MQ=255)
                                                              ccATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGa                                                                                      <  1:1062358/150‑1 (MQ=255)
                                                                      gCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGgcgc                                                                              <  1:299995/150‑1 (MQ=255)
                                                                       cacaAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCt                                                                             >  2:652962/1‑150 (MQ=255)
                                                                       cacaAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCt                                                                             <  2:1417957/150‑1 (MQ=255)
                                                                              cAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTTAATAAGGAAAGAAGATGCTTGAACAATTATTCAGGCCCGACTCTGTCGACCAGGCGCTGGATCTg                                                                      >  1:584171/1‑150 (MQ=255)
                                                                               aGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTg                                                                                                                  >  2:640465/1‑105 (MQ=255)
                                                                               aGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTg                                                                                                                  <  1:640465/105‑1 (MQ=255)
                                                                                gAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGaa                                                                    >  1:376706/1‑150 (MQ=255)
                                                                                           gCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTAcc                                                         <  2:1576211/150‑1 (MQ=255)
                                                                                              aGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAgg                                                      <  1:958640/150‑1 (MQ=255)
                                                                                                   gtgGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGaa                                                 <  2:362514/150‑1 (MQ=255)
                                                                                                     ggCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGc                                               <  1:893521/150‑1 (MQ=255)
                                                                                                     ggCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGc                                               <  1:887675/150‑1 (MQ=255)
                                                                                                         tCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTc                                           <  1:1083271/150‑1 (MQ=255)
                                                                                                           ggatggatGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTg                                         >  1:205960/1‑150 (MQ=255)
                                                                                                            gatggatgCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTgg                                        >  2:1362020/1‑150 (MQ=255)
                                                                                                              tggatgCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGtt                                      <  2:663927/150‑1 (MQ=255)
                                                                                                                gatgCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCg                                    <  2:386252/150‑1 (MQ=255)
                                                                                                                gaagcgcTGGCATAAATGACAGATGCCCTCTTCCCGCAGCGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCAACCAGCATCAACCCGTCTGGTTCg                                    >  2:561140/4‑150 (MQ=255)
                                                                                                                  tgCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGcc                                  <  1:1301315/150‑1 (MQ=255)
                                                                                                                     gcTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCggg                               <  1:1179927/150‑1 (MQ=255)
                                                                                                                       tGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCggggg                             <  1:668880/150‑1 (MQ=255)
                                                                                                                          cATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAATTTTCAGGCCCGACACTGTCGAACAGGCGCTGGAACTCAACCGCCGCTACCAGCAACAACCCCTCAGCTTCGCCCGGCgcac                          >  2:1385866/1‑149 (MQ=255)
                                                                                                                            tAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCTGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgcagca                        >  1:430896/1‑150 (MQ=255)
                                                                                                                            tAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgcagca                        >  1:1056772/1‑150 (MQ=255)
                                                                                                                                tGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACt                    >  2:1209540/1‑150 (MQ=255)
                                                                                                                                 gACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTc                   <  2:1483742/150‑1 (MQ=255)
                                                                                                                                 gACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGGCTGGTTCGCCGGGGGCAGCAAACTc                   <  1:115852/150‑1 (MQ=255)
                                                                                                                                   cAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCaa                 <  2:1512361/150‑1 (MQ=255)
                                                                                                                                    aGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAAt                <  1:1448513/150‑1 (MQ=255)
                                                                                                                                     gATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATg               >  1:908769/1‑150 (MQ=255)
                                                                                                                                      aTGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGc              <  1:830196/150‑1 (MQ=255)
                                                                                                                                      aTGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGc              <  1:73977/150‑1 (MQ=255)
                                                                                                                                       tGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCt             >  2:1535248/1‑150 (MQ=255)
                                                                                                                                       tGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCt             >  2:878752/1‑150 (MQ=255)
                                                                                                                                       tGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCt             <  2:1249997/150‑1 (MQ=255)
                                                                                                                                         cccTCTTCCCTGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac           >  1:1101732/1‑150 (MQ=255)
                                                                                                                                         cccTCTTCCCGCAGGGAAGAGGTTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCACCAAAAGCAATGCTac           >  2:790484/1‑150 (MQ=255)
                                                                                                                                          ccTCTTCCCTGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGCATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTaca          >  1:1375120/1‑150 (MQ=255)
                                                                                                                                          ccTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTaca          >  1:1064326/1‑150 (MQ=255)
                                                                                                                                           ctctTCCCTGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTacac         >  1:1291144/1‑150 (MQ=255)
                                                                                                                                                cccGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAAcc    >  2:352882/1‑150 (MQ=255)
                                                                                                                                                  cGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCAACACCAACCCg  >  2:1297929/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCG  >  NC_000913/3020363‑3020658

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