Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,508,816 C→T 12.2% intergenic (+74/+18) 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,8160CT12.2% 245.3 / 19.1 90intergenic (+74/+18)ydcP/yncJputative peptidase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (44/35);  new base T (7/4);  total (51/39)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

CCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTGATGATCCAC  >  NC_000913/1508668‑1508962
                                                                                                                                                    |                                                                                                                                                  
ccGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCa                                                                                                                                                   >  1:7299/1‑150 (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:154721/1‑150 (MQ=255)
                                                                     gggTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTcgc                                                                              >  1:1475824/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                                                                        <  1:711014/150‑1 (MQ=255)
                                                                           aaCGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAAc                                                                        <  2:363543/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:2947/1‑150 (MQ=255)
                                                                                                                 gCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCg                                  >  1:1705427/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)
                                                                                                                                               cgccgcATCCGGCATCAAATGCCTCAACGTTGTCCCGCTTCTGGTTTGCGCCAATTTTGGCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACTGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTGATGCTcc    >  1:1028549/1‑150 (MQ=255)
                                                                                                                                               cgccgcATCCGGCATCAAATGACTCAACTTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTGATGATcc    <  1:1268082/150‑1 (MQ=255)
                                                                                                                                               cgccgcATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTGATGATcc    <  2:1548712/150‑1 (MQ=255)
                                                                                                                                                gccgcATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTGATGATCca   >  2:980781/1‑150 (MQ=255)
                                                                                                                                                gccgcATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTGATGATCca   >  1:17330/1‑150 (MQ=255)
                                                                                                                                                gccgcATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTGATGATCca   <  1:470698/150‑1 (MQ=255)
                                                                                                                                                 ccgcATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTGATGATCcac  >  1:560037/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTGATGATCCAC  >  NC_000913/1508668‑1508962

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: