Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,508,802 A→G 11.1% intergenic (+60/+32) ydcP → / ← yncJ putative peptidase/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,508,8020AG11.1% 261.7 / 20.8 99intergenic (+60/+32)ydcP/yncJputative peptidase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (48/40);  new base G (6/5);  total (54/45)
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.77e-01

CATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGC  >  NC_000913/1508655‑1508951
                                                                                                                                                   |                                                                                                                                                     
caTACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAg                                                                                                                                                     <  1:401326/150‑1 (MQ=255)
  tACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACCACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGc                                                                                                                                                   >  2:302920/1‑150 (MQ=255)
   aCTGAAAATGCCTTTACCGGGAAGCGTAGTGGCATCCTTAAGTCCGGATTAGCTGCTGAAAACATTGCCGAAGCGAAAATGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCg                                                                                                                                                  <  1:128991/150‑1 (MQ=255)
   aCTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCg                                                                                                                                                  <  2:935574/150‑1 (MQ=255)
    cTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGt                                                                                                                                                 >  1:162347/1‑150 (MQ=255)
    cTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGt                                                                                                                                                 <  2:1586523/150‑1 (MQ=255)
     tGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGtt                                                                                                                                                >  1:1009948/1‑150 (MQ=255)
      gAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTGATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTa                                                                                                                                               >  1:279474/1‑150 (MQ=255)
        aaaTGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTAcg                                                                                                                                             <  1:1190116/150‑1 (MQ=255)
          ataccGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTAcgcc                                                                                                                                           <  1:1139879/147‑1 (MQ=255)
           tGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTAcgccg                                                                                                                                          <  1:1317826/150‑1 (MQ=255)
           tGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTAcgccg                                                                                                                                          >  2:307355/1‑150 (MQ=255)
             ccGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCa                                                                                                                                        >  1:7299/1‑150 (MQ=255)
                ttACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGTGTAAAACGCCAGTTTTCTGTTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGATGCGGCGTAAcgccttaccc                                                                                                                                     <  2:881573/150‑7 (MQ=255)
                ttACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATcc                                                                                                                                     <  2:219686/150‑1 (MQ=255)
                 tACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGTGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCg                                                                                                                                    <  2:1571132/150‑1 (MQ=255)
                 tACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCg                                                                                                                                    <  2:1590847/150‑1 (MQ=255)
                  aCCGTGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGTGTAAGGCGTTACGCCGCATCCgg                                                                                                                                   <  2:1053011/150‑1 (MQ=255)
                   cccggAAGCGTAGTCGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGAGAAAAGTGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGATGCGGCGTAACGCCTTACCCGGc                                                                                                                                  <  1:230763/147‑1 (MQ=255)
                    cGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCa                                                                                                                                 <  2:688083/150‑1 (MQ=255)
                     gggAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCAt                                                                                                                                <  2:703974/150‑1 (MQ=255)
                     gggAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGATGCGGCGTAACGCCTTACCCGGCag                                                                                                                                <  1:340439/150‑2 (MQ=255)
                      ggAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGATGCGGCGTAACGCCTTACCCGGCagg                                                                                                                               <  2:623386/150‑3 (MQ=255)
                       gAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGATGCGGCGTAACGCCTTACCCGGCaggc                                                                                                                              <  1:89579/150‑4 (MQ=255)
                           cGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATg                                                                                                                          <  2:778028/150‑1 (MQ=255)
                             tAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGAc                                                                                                                        <  1:398814/150‑1 (MQ=255)
                                 ggCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCaa                                                                                                                    >  1:347961/1‑150 (MQ=255)
                                    aTCCGTAAGTCCTGATTAGCTGCTGAAAACATTGCCGAAGCGAAAATGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGt                                                                                                                 <  1:1331619/150‑1 (MQ=255)
                                       cGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGt                                                                                                              >  2:1237339/1‑150 (MQ=255)
                                         tAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTcc                                                                                                            >  1:363543/1‑150 (MQ=255)
                                           aGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCTCATCCGGCATCAAATGACTCAACGTTGTCCCg                                                                                                          >  1:1257647/1‑150 (MQ=255)
                                                ggATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCTAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCt                                                                                                     >  2:231343/1‑150 (MQ=255)
                                                  aTGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTgg                                                                                                   <  2:1601870/150‑1 (MQ=255)
                                                    gAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACCACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGtt                                                                                                 >  1:1142904/1‑150 (MQ=255)
                                                    gAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATTTTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGtt                                                                                                 <  2:596930/150‑1 (MQ=255)
                                                      gctgctGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTg                                                                                               >  1:138203/1‑150 (MQ=255)
                                                         gctgAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTgcgg                                                                                            >  1:346084/1‑149 (MQ=255)
                                                         gctgAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTgcgc                                                                                            >  2:1108026/1‑150 (MQ=255)
                                                             aaaaCATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTTAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCt                                                                                                     <  1:1343227/137‑1 (MQ=255)
                                                             aaaaCATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCt                                                                                                     >  2:1343227/1‑137 (MQ=255)
                                                                aCATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGAttttt                                                                                     >  2:470698/1‑150 (MQ=255)
                                                                 cATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTg                                                                                    >  2:329572/1‑150 (MQ=255)
                                                                 aaTTGCAGAAGCGAAAAGGGTAATACGCGAGTTTTATTGTTACTCACAACTAATTGAATCTGCATGATATTGCCTGCCGGTTAAGGCGTTACGACGCATCCGGCATCAAATGACTCAACGTGTTCGCTCTTCTGGTTTGCGCGATTTTTg                                                                                    <  1:275198/149‑1 (MQ=255)
                                                                     gCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAg                                                                                >  1:1544570/1‑150 (MQ=255)
                                                                       cGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTa                                                                              >  1:606177/1‑150 (MQ=255)
                                                                        ggaGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTaa                                                                             >  1:541089/3‑150 (MQ=255)
                                                                          aGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTaaaa                                                                           >  2:791529/1‑150 (MQ=255)
                                                                           gCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTaaaaa                                                                          >  1:1017729/1‑150 (MQ=255)
                                                                             gAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCCCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATg                                                                        >  2:188334/1‑150 (MQ=255)
                                                                             gAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTAGTGAATCTGCATCATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCAGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATg                                                                        >  1:447202/1‑150 (MQ=255)
                                                                                  gggTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTcgc                                                                   >  1:1475824/1‑150 (MQ=255)
                                                                                  gggTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTcgc                                                                   >  1:154721/1‑150 (MQ=255)
                                                                                  gggTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAACGCTcgc                                                                   <  1:1534552/150‑1 (MQ=255)
                                                                                  gggTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGACCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTcgc                                                                   <  2:373788/150‑1 (MQ=255)
                                                                                  gggTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCCCCTCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTcgc                                                                   >  2:92988/1‑150 (MQ=255)
                                                                                   ggTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGTCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTcgcg                                                                  >  2:143761/1‑150 (MQ=255)
                                                                                   ggTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTcgcg                                                                  >  1:1017064/1‑150 (MQ=255)
                                                                                   ggTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTcgcg                                                                  >  2:87421/1‑150 (MQ=255)
                                                                                    gTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGACTCTGCATGATATTGCCTCCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTACTGGTTTGCGCGATTTTTGACAGAAAAAATGCTCgcgc                                                                 >  2:265539/1‑150 (MQ=255)
                                                                                      aaaaCGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCCACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCaa                                                               >  2:1087608/1‑150 (MQ=255)
                                                                                        aaCGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAAc                                                             <  2:363543/150‑1 (MQ=255)
                                                                                        aaCGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAAc                                                             <  1:711014/150‑1 (MQ=255)
                                                                                         aCGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGTCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAAcc                                                            <  2:357205/150‑1 (MQ=255)
                                                                                               gTTTTCTGGATACTCACAACTAATAGAATCTGCATGATATTGCCTGCCGTGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGGCCCGCTGCCGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTccgc                                                      <  1:255054/150‑1 (MQ=255)
                                                                                                   tCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGAc                                                  <  2:1250641/150‑1 (MQ=255)
                                                                                                       gTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATCACTAAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCAAGTAAAAATGCTCGCGCAAACCTACCGCCGACTCtt                                              >  1:475970/1‑150 (MQ=255)
                                                                                                          aCTCACAACTAATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGTTTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCg                                           <  1:856430/150‑1 (MQ=255)
                                                                                                            tCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGcc                                         <  2:1017064/150‑1 (MQ=255)
                                                                                                            tCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGcc                                         >  1:1588635/1‑150 (MQ=255)
                                                                                                             cacaACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCATTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCca                                        >  1:344529/1‑150 (MQ=255)
                                                                                                             cacaACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCca                                        >  2:763588/1‑150 (MQ=255)
                                                                                                              acaACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCcac                                       >  1:788974/1‑150 (MQ=255)
                                                                                                              acaACTTATTGAATCTGCATGATATTGCCAGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCcac                                       <  2:1593448/150‑1 (MQ=255)
                                                                                                                aaCTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAg                                     <  1:231343/150‑1 (MQ=255)
                                                                                                                 aCTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGc                                    >  1:386721/1‑150 (MQ=255)
                                                                                                                   ttATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCAc                                  >  1:699298/1‑150 (MQ=255)
                                                                                                                   ttATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTAAGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCAc                                  <  1:302920/150‑1 (MQ=255)
                                                                                                                    tATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACg                                 <  1:1587679/150‑1 (MQ=255)
                                                                                                                    tATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACg                                 >  2:885490/1‑150 (MQ=255)
                                                                                                                    tATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACg                                 >  2:490470/1‑150 (MQ=255)
                                                                                                                     tttGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGc                                <  1:341157/149‑1 (MQ=255)
                                                                                                                        gAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCAGCTTCTGGTTTGCGCGATTTTTTCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCTCCACAGCACGCAAt                             >  1:785610/1‑150 (MQ=255)
                                                                                                                              gCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGCCTCTTCCGCCACAGCACGCAATTCATCg                       >  2:512421/1‑150 (MQ=255)
                                                                                                                              gCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCg                       >  1:1705427/1‑150 (MQ=255)
                                                                                                                              gCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCg                       >  1:2947/1‑150 (MQ=255)
                                                                                                                               cATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGc                      <  2:1681239/150‑1 (MQ=255)
                                                                                                                               cATGATATTACCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGc                      >  2:593347/1‑150 (MQ=255)
                                                                                                                                   atatTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTc                  >  1:1118267/1‑150 (MQ=255)
                                                                                                                                    catttgatgccGGATGCGGCGTAACGCCTTACCCGGCATCAACTGACTCGACGTTGTCCCGCTTCTGGTTTGCGCGAGTTTTGCCAGTAAAAATGCTCGCGCAAGCCGTCCGCCGACTCGTCCGCCACAGCACGCAATTCATCGCTGTcc                 >  1:69670/8‑150 (MQ=25)
                                                                                                                                     atttgatgccGGATGCGGCGTAACGCCTTACCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAACCCTTCCGCCGACGCTTCCGCCACAGCAGGCTATTCATCGCTGTCCg                >  1:833178/7‑150 (MQ=37)
                                                                                                                                       tgcctgccGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCt              <  1:661391/150‑1 (MQ=255)
                                                                                                                                        tgatgccGGATGCGGCGTAACGCCTTACCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCtt             >  1:361678/4‑150 (MQ=255)
                                                                                                                                        tgatgccGGATGCGGCGTAACGCCTTACCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATAGCTGTCCGCtt             >  2:130808/4‑150 (MQ=255)
                                                                                                                                         gatgccGGATGCGGCGTAACGCCTTACCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTc            >  1:417420/3‑150 (MQ=255)
                                                                                                                                         gatgccGGATGCGGCGTAACGCCTTACCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGACGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTc            >  2:693045/3‑150 (MQ=255)
                                                                                                                                         cctgccGGGTAATGCGTTACGCCGCTTCCGGCATCAAATGACTCACCGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAATTGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTc            <  2:665594/150‑1 (MQ=255)
                                                                                                                                         cctgccGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTc            <  2:747565/150‑1 (MQ=255)
                                                                                                                                          ctgccGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCa           <  2:768419/150‑1 (MQ=255)
                                                                                                                                            gccGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATg         <  2:1208907/150‑1 (MQ=255)
                                                                                                                                                   aaGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGc  <  2:1539019/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGC  >  NC_000913/1508655‑1508951

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