Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,915,868 A→T 5.3% intergenic (+254/‑35) fklB‑II → / → PP_1715 FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,915,8680AT 5.3% 236.7 / 3.9 76intergenic (+254/‑35)fklB‑II/PP_1715FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (43/29);  new base T (2/2);  total (45/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

GAAGTAGGATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTCG  >  NC_002947/1915720‑1916014
                                                                                                                                                    |                                                                                                                                                  
gaagTAGGATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAAt                                                                                                                                                   <  1:178246/150‑1 (MQ=255)
   gTAGGATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATaaa                                                                                                                                                >  2:920145/1‑150 (MQ=255)
      ggATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCTTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAAtgt                                                                                                                                             <  1:1031646/150‑1 (MQ=255)
      ggATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAAtgt                                                                                                                                             <  1:920145/150‑1 (MQ=255)
      ggATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAAtgt                                                                                                                                             <  1:1131218/150‑1 (MQ=255)
      ggATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAAtgt                                                                                                                                             <  1:1934603/150‑1 (MQ=255)
      ggATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAAtgt                                                                                                                                             >  1:1735137/1‑150 (MQ=255)
        aTGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTgg                                                                                                                                           <  1:1866180/150‑1 (MQ=255)
        aTGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTgg                                                                                                                                           >  1:1685226/1‑150 (MQ=255)
             tCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAAc                                                                                                                                      <  1:1200189/150‑1 (MQ=255)
               tGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCt                                                                                                                                    >  1:2022842/1‑150 (MQ=255)
                  cGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACATTTATTGTGGATAACCttgt                                                                                                                                 <  1:600633/150‑3 (MQ=255)
                   gCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTc                                                                                                                                >  1:583334/1‑150 (MQ=255)
                      aaTCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCctc                                                                                                                             >  1:514263/1‑150 (MQ=255)
                           tttGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGa                                                                                                                        <  2:514263/150‑1 (MQ=255)
                                gccgcAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCg                                                                                                                          >  1:1017274/1‑143 (MQ=255)
                                gccgcAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCg                                                                                                                          <  2:1017274/143‑1 (MQ=255)
                                 ccgcAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTcc                                                                                                                  >  1:1814388/1‑150 (MQ=255)
                                   gcAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCAt                                                                                                                <  2:2313967/150‑1 (MQ=255)
                                   gcAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCAt                                                                                                                <  1:1041807/150‑1 (MQ=255)
                                       gTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAgc                                                                                                            <  1:105930/150‑1 (MQ=255)
                                           tCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCAc                                                                                                        >  2:300079/1‑150 (MQ=255)
                                             aTGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGTGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAgc                                                                                                            <  2:315092/144‑1 (MQ=255)
                                             aTGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGTGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAgc                                                                                                            >  1:315092/1‑144 (MQ=255)
                                             aTGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAg                                                                                                             >  1:1114522/1‑143 (MQ=255)
                                             aTGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAg                                                                                                             <  2:1114522/143‑1 (MQ=255)
                                              tGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCt                                                                                                     <  1:140998/150‑1 (MQ=255)
                                                gCGTGAAACGACTGTGTAAAAAACGCCCGATCTGGCCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTgg                                                                                                   >  2:715053/1‑150 (MQ=255)
                                                       aCGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTg                                                                                            <  2:90702/150‑1 (MQ=255)
                                                          aCTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGa                                                                                                                        <  2:2110815/119‑1 (MQ=255)
                                                          aCTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGa                                                                                                                        >  1:2110815/1‑119 (MQ=255)
                                                             gtgtAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGtt                                                                                      <  2:1814388/150‑1 (MQ=255)
                                                              tgtAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGttt                                                                                     >  2:1001614/1‑150 (MQ=255)
                                                                    aaaCGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGcc                                                                               >  1:366068/1‑150 (MQ=255)
                                                                    aaaCGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGcc                                                                               <  1:1792602/150‑1 (MQ=255)
                                                                        gCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCctct                                                                                                                            <  2:184578/101‑1 (MQ=255)
                                                                        gCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCctct                                                                                                                            >  1:184578/1‑101 (MQ=255)
                                                                            gATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCg                                                                       >  1:1228649/1‑150 (MQ=255)
                                                                             aTCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGa                                                                      >  1:1658604/1‑150 (MQ=255)
                                                                              tCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAg                                                                     >  2:811854/1‑150 (MQ=255)
                                                                               cTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAgc                                                                    <  2:1735137/150‑1 (MQ=255)
                                                                                tGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACATTTATTGTGGATAACCttgtga                                                                                                                               >  1:1386567/1‑86 (MQ=25)
                                                                                tGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACATTTATTGTGGATAACCttgtga                                                                                                                               <  2:1386567/90‑5 (MQ=25)
                                                                                     gggCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGc                                                              >  1:1477088/1‑150 (MQ=255)
                                                                                      ggCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCt                                                             >  2:924258/1‑150 (MQ=255)
                                                                                       gCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTcc                                                                                                                  <  1:2211982/96‑1 (MQ=255)
                                                                                       gCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTcc                                                                                                                  >  2:2211982/1‑96 (MQ=255)
                                                                                          ccTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCtttt                                                                                                                                 >  2:2201102/1‑78 (MQ=255)
                                                                                          ccTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCtttt                                                                                                                                 <  1:2201102/78‑1 (MQ=255)
                                                                                           cTTGCCTGGCAGGCACCTGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCcg                                                       >  1:325062/1‑150 (MQ=255)
                                                                                             tGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCcgc                                                      <  1:300079/150‑1 (MQ=255)
                                                                                              gCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCcgcg                                                     >  2:844697/1‑150 (MQ=255)
                                                                                                           cTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGc                                        >  1:1318002/1‑150 (MQ=255)
                                                                                                            tGGGGGGGGAAACGGGGCCCTGTTCACAGGGTTATCCACAATACCTGTGAAGAGACGTTTCCTCTCGGAGGCTCCGTGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGcc                                       <  2:64918/150‑1 (MQ=255)
                                                                                                                    gAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTgttgtt                               >  2:1366368/1‑150 (MQ=255)
                                                                                                                           cGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTg                        <  1:2288741/150‑1 (MQ=255)
                                                                                                                            gCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTgg                       >  1:1724808/1‑150 (MQ=255)
                                                                                                                             cccTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGgc                      >  1:2173993/1‑150 (MQ=255)
                                                                                                                              ccTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCAGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGgcg                     >  2:1832093/1‑150 (MQ=255)
                                                                                                                               cTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAAtt                                                                                               <  2:1417182/75‑1 (MQ=255)
                                                                                                                               cTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAAtt                                                                                               >  1:1417182/1‑75 (MQ=255)
                                                                                                                                     acaAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGa                                                                      <  1:33572/94‑1 (MQ=255)
                                                                                                                                     acaAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGa                                                                      >  2:33572/1‑94 (MQ=255)
                                                                                                                                     acaAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAg              >  2:1994368/1‑150 (MQ=255)
                                                                                                                                     acaAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCAGGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAg              >  2:1687717/1‑150 (MQ=255)
                                                                                                                                       aaGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGc            >  2:1319863/1‑150 (MQ=255)
                                                                                                                                       aaGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGc            >  1:251709/1‑150 (MQ=255)
                                                                                                                                          gTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCATGCCGCTGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCcg         >  1:737294/1‑150 (MQ=255)
                                                                                                                                           ttATCCACATTTATTGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCcgc        >  2:915800/1‑150 (MQ=255)
                                                                                                                                             aTCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCgcgc      <  2:1685226/150‑1 (MQ=255)
                                                                                                                                               ccACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCt    >  1:407142/1‑150 (MQ=255)
                                                                                                                                                cacaATCAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGCCGTTTGCCGCCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTc   >  2:1954977/1‑150 (MQ=255)
                                                                                                                                                 acaATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTCg  >  2:583689/1‑150 (MQ=255)
                                                                                                                                                 acaATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTCg  >  2:1521079/1‑150 (MQ=255)
                                                                                                                                                 acaATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTGCCTGTCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCCGTCAAAGGCCCGCGCCTCg  >  1:1329365/1‑150 (MQ=255)
                                                                                                                                                 acaATAAACGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTCTTGTTCGCGGAGGCGCGCACATAGCCGCGCCTCg  >  2:441904/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAAGTAGGATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTCG  >  NC_002947/1915720‑1916014

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