Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,429,120 Δ1 bp 7.5% coding (186/621 nt) fklB → FKBP‑type peptidyl‑prolyl cis‑trans isomerase FklB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,429,1200G. 7.5% 529.2 / 34.4 159coding (186/621 nt)fklBFKBP‑type peptidyl‑prolyl cis‑trans isomerase FklB
Reads supporting (aligned to +/- strand):  ref base G (78/69);  new base . (4/8);  total (82/77)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 164 total aligned reads displayed.

GCAAGCTACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGG  >  NC_000913/4428971‑4429269
                                                                                                                                                      |                                                                                                                                                      
gcaagcTACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTg                                                                                                                                                        <  1:522928/150‑1 (MQ=255)
gcaagcTACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTg                                                                                                                                                        <  2:491719/150‑1 (MQ=255)
  aagcTACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCg                                                                                                                                                      <  1:940761/150‑1 (MQ=255)
   agcTACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGt                                                                                                                                                     <  2:1229187/150‑1 (MQ=255)
     cTACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGa                                                                                                                                                   >  2:598947/1‑150 (MQ=255)
       aCGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGaaa                                                                                                                                                 >  1:766892/1‑150 (MQ=255)
          gCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATcc                                                                                                                                              >  2:879660/1‑150 (MQ=255)
          gCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATcc                                                                                                                                              >  1:1273505/1‑150 (MQ=255)
           cATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCa                                                                                                                                             <  2:300608/150‑1 (MQ=255)
            aTTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCAc                                                                                                                                            <  1:100634/150‑1 (MQ=255)
               ggTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAg                                                                                                                                        >  2:12745/1‑150 (MQ=255)
               ggTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAg                                                                                                                                        >  1:802918/1‑150 (MQ=255)
                  ttGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTGGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAgcgc                                                                                                                                     <  1:553010/150‑1 (MQ=255)
                      aGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCg                                                                                                                                 >  1:1296550/1‑150 (MQ=255)
                       ggTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCTCGTGGATTTCACGCAgcgcgcga                                                                                                                                 <  2:597565/150‑3 (MQ=38)
                       ggTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGa                                                                                                                                >  1:1194852/1‑150 (MQ=255)
                         tAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATg                                                                                                                              >  2:1028356/1‑150 (MQ=255)
                          aGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTGGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGc                                                                                                                             <  2:262351/150‑1 (MQ=255)
                           gggCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCTCGTGGATTTCACGCAGCGCGCGATGca                                                                                                                            <  1:68489/150‑2 (MQ=37)
                           gggCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCTCGTGGATTTCACGCAGCGCGCGATGca                                                                                                                            <  1:916101/150‑2 (MQ=37)
                           gggCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCAGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGcc                                                                                                                            >  2:771652/1‑150 (MQ=255)
                            ggCAACAACTGATTGAATCTGGTCTGGAATGGCTTCTGCCATAAGCACTGTTTGCAGGTATTGCCGATGCGCTTGACGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATTGCGCGATCGTGGATTTCACGCAGCGCGCGATGcacc                                                                                                                          <  2:544633/150‑4 (MQ=255)
                            ggCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCa                                                                                                                                             >  2:390866/1‑133 (MQ=255)
                            ggCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCa                                                                                                                                             <  1:390866/133‑1 (MQ=255)
                               aacaacTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTtc                                                                                                                        >  1:335293/1‑150 (MQ=255)
                                 caacTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCTCGTGGATTTCACGCAGCGCGCGATGcaccacat                                                                                                                      <  2:1271046/150‑8 (MQ=25)
                                 caacTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCAGCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTtcgt                                                                                                                      <  2:1003402/150‑1 (MQ=255)
                                   acTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTtcgtcg                                                                                                                    >  1:1259622/1‑150 (MQ=255)
                                            aTCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCTCGTGGATTTCACGCAGCGCGCGATGcaccaca                                                                                                                       <  2:1120366/138‑7 (MQ=25)
                                             tCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTc                                                                                                          >  1:1381706/1‑150 (MQ=255)
                                              cTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTcc                                                                                                         <  1:771652/150‑1 (MQ=255)
                                              cTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTcc                                                                                                         >  1:872444/1‑150 (MQ=255)
                                              cTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTcc                                                                                                         <  2:1055856/150‑1 (MQ=255)
                                                  ccTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGc                                                                                                     <  1:1276741/150‑1 (MQ=255)
                                                  ccTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGc                                                                                                     >  1:1074614/1‑150 (MQ=255)
                                                    tGGAAGGGTTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGa                                                                                                   >  1:247762/1‑150 (MQ=255)
                                                           gctgctGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctgc                                                                                            <  2:1002579/150‑1 (MQ=255)
                                                           gctgctGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctgc                                                                                            >  2:484649/1‑150 (MQ=255)
                                                           gctgctGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGGCAGCGTTTCCAGGCGATGgctgc                                                                                            <  1:318071/150‑1 (MQ=255)
                                                            ctgctgCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctgct                                                                                           >  2:736225/1‑150 (MQ=255)
                                                               ctgcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGaa                                                                                        <  1:1122360/150‑1 (MQ=255)
                                                                     gAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGtgtg                                                                                  <  2:1074614/150‑1 (MQ=255)
                                                                     gAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGtgtg                                                                                  >  1:277610/1‑150 (MQ=255)
                                                                      aaGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGa                                                                                 <  1:1006965/150‑1 (MQ=255)
                                                                            tGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATAcc                                                                           <  2:766892/150‑1 (MQ=255)
                                                                                tGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGga                                                                       <  2:659106/150‑1 (MQ=255)
                                                                                   aGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGgaaga                                                                    <  1:183488/150‑1 (MQ=255)
                                                                                    ggTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGgaagaa                                                                   >  2:1169699/1‑150 (MQ=255)
                                                                                       aTTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAAc                                                                <  2:1396501/150‑1 (MQ=255)
                                                                                       aTTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAAc                                                                >  1:27961/1‑150 (MQ=255)
                                                                                          gCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGga                                                                       <  1:351511/140‑1 (MQ=255)
                                                                                          gCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGcc                                                             >  1:1050393/1‑150 (MQ=255)
                                                                                          gCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGcc                                                             >  2:946460/1‑150 (MQ=255)
                                                                                             gATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCaaa                                                          >  1:963043/1‑150 (MQ=255)
                                                                                               tGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCACGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCaaaaa                                                        >  1:498108/1‑150 (MQ=255)
                                                                                                gcgcTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCaaaaaa                                                       <  2:1195/150‑1 (MQ=255)
                                                                                                 cgcTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaag                                                      >  2:956606/1‑150 (MQ=255)
                                                                                                  gcTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaaga                                                     >  1:752454/1‑150 (MQ=255)
                                                                                                    tGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaagaag                                                   <  1:1023056/150‑1 (MQ=255)
                                                                                                    tGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaagaag                                                   <  2:835363/150‑1 (MQ=255)
                                                                                                        agtcAAACATCCGGCTTTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATTCCGTTCGTCGTCAGCGTTTCCAGTCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGgtgt                                               <  2:1355992/147‑1 (MQ=255)
                                                                                                         ggCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCATGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGtgtg                                              <  2:286249/150‑1 (MQ=255)
                                                                                                         ggCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGtgtg                                              >  1:1129227/1‑150 (MQ=255)
                                                                                                         ggCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGtgtg                                              >  1:824819/1‑150 (MQ=255)
                                                                                                          gCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTGCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAACGCCAAAAAAGAAGGTGTGaa                                            <  2:464646/150‑1 (MQ=255)
                                                                                                           cAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGaa                                            >  1:85156/1‑150 (MQ=255)
                                                                                                           cAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGaa                                            >  1:180103/1‑150 (MQ=255)
                                                                                                           cAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGACGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGaa                                            <  1:379151/150‑1 (MQ=255)
                                                                                                           cAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTAAGCGTTTCCAGGCGATGgc                                                                                               >  1:1374196/1‑99 (MQ=255)
                                                                                                               cATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGc                                        >  2:924268/1‑150 (MQ=255)
                                                                                                               cATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGc                                        >  2:431243/1‑150 (MQ=255)
                                                                                                                  ccGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAAt                                                                              <  2:529751/109‑1 (MQ=255)
                                                                                                                  ccGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCAcc                                     >  2:913141/1‑150 (MQ=255)
                                                                                                                   cGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCg                                    >  1:301920/1‑150 (MQ=255)
                                                                                                                    ggCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGa                                   <  2:1029487/150‑1 (MQ=255)
                                                                                                                        gTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCt                               <  2:1296550/150‑1 (MQ=255)
                                                                                                                        gTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCt                               >  2:1075334/1‑150 (MQ=255)
                                                                                                                           ccGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGc                            >  2:190148/1‑150 (MQ=255)
                                                                                                                             ggTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAATGTGTGAATAGCACCGAATCTGGCCt                          <  2:335293/150‑1 (MQ=255)
                                                                                                                               ttGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGATTAGCACCGAATCTGGCCTGc                        <  2:721376/150‑1 (MQ=255)
                                                                                                                                    gaacggCATCGGCGCGCTCGTGGATTTCACGCAGCGCGCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGCAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTc                   >  1:499140/6‑150 (MQ=25)
                                                                                                                                     tggtgCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTcc                  >  1:146351/1‑150 (MQ=255)
                                                                                                                                     tggtgCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTcc                  >  1:414781/1‑150 (MQ=255)
                                                                                                                                      ggtgCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCcg                 >  1:429932/1‑150 (MQ=255)
                                                                                                                                       gtgCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCcgc                <  2:851898/150‑1 (MQ=255)
                                                                                                                                       gtgCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCcgc                <  2:85156/150‑1 (MQ=255)
                                                                                                                                        tgCATCGCGCGCTGAGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCcgcg               >  2:956070/1‑150 (MQ=255)
                                                                                                                                         gCATCGGCGCGCTCGTGGATTTCACGCAGCGCGCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGcc                           <  2:38203/137‑1 (MQ=37)
                                                                                                                                         gCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAAt                                 >  1:744613/1‑131 (MQ=255)
                                                                                                                                         gCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGt              <  1:941656/150‑1 (MQ=255)
                                                                                                                                          cATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTg             <  1:1251687/150‑1 (MQ=255)
                                                                                                                                          cATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTg             <  2:988475/150‑1 (MQ=255)
                                                                                                                                          cATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTg             <  1:484649/150‑1 (MQ=255)
                                                                                                                                             cgcgcgCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCTTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATc          <  2:1259622/150‑1 (MQ=255)
                                                                                                                                             cgcgcgCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATc          >  2:1412570/1‑150 (MQ=255)
                                                                                                                                               gcgcgcTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCa         <  1:1028356/150‑1 (MQ=255)
                                                                                                                                                cgcgcTGCCTGAAATCCACGAGCGCGCCGATCCCGTTCGTCGTCAGCGTTTCCAGGCGATCGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCaa        >  2:521811/1‑150 (MQ=255)
                                                                                                                                                  cgcTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAAcc      <  1:1162895/150‑1 (MQ=255)
                                                                                                                                                      gCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAggg  <  2:247762/150‑1 (MQ=255)
                                                                                                                                                      gCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAggg  >  1:1092859/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GCAAGCTACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGG  >  NC_000913/4428971‑4429269

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