Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,508,800 G→A 11.1% intergenic (+58/+34) 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,8000GA11.1% 274.5 / 17.1 99intergenic (+58/+34)ydcP/yncJputative peptidase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (48/39);  new base A (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 = 9.33e-01

CATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATG  >  NC_000913/1508655‑1508944
                                                                                                                                                 |                                                                                                                                                
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)
                                                                                                                                                 |                                                                                                                                                
CATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATG  >  NC_000913/1508655‑1508944

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: