Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,508,809 T→A 12.5% intergenic (+67/+25) 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,8090TA12.5% 235.4 / 21.1 96intergenic (+67/+25)ydcP/yncJputative peptidase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (46/38);  new base A (7/5);  total (53/43)
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.43e-01

TGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTGA  >  NC_000913/1508660‑1508954
                                                                                                                                                     |                                                                                                                                                 
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                                                                     >  2:87421/1‑150 (MQ=255)
                                                                              ggTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTcgcg                                                                     >  1:1017064/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)
                                                                                                                                         cggatgcggCGTAACGCCTTACCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGAc          >  2:673319/8‑150 (MQ=255)
                                                                                                                                              aaGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGc     <  2:1539019/150‑1 (MQ=255)
                                                                                                                                               aGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCt    <  1:1108026/150‑1 (MQ=255)
                                                                                                                                                 gCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCtga  >  2:1615082/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTGA  >  NC_000913/1508660‑1508954

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: