Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 824,949 T→C 6.7% L107P (CTG→CCG)  ybhQ → putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913824,9490TC 6.7% 220.4 / 4.8 74L107P (CTG→CCG) ybhQputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (39/30);  new base C (3/2);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GTTTCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAG  >  NC_000913/824800‑825097
                                                                                                                                                     |                                                                                                                                                    
gtTTCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCt                                                                                                                                                      >  1:1189728/1‑150 (MQ=255)
gtTTCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCt                                                                                                                                                      >  1:689278/1‑150 (MQ=255)
  ttCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTgg                                                                                                                                                    >  2:1012199/1‑150 (MQ=255)
   tCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGc                                                                                                                                                   >  2:244323/1‑150 (MQ=255)
   tCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCGGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGc                                                                                                                                                   <  1:629317/150‑1 (MQ=255)
   acAGCGCCACACCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGCCCACCACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCAGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTgct                                                                                                                                                   >  2:756584/2‑148 (MQ=255)
    cAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCa                                                                                                                                                  >  1:151527/1‑150 (MQ=255)
     aGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCaa                                                                                                                                                 >  1:96009/1‑150 (MQ=255)
      gcgcCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAAt                                                                                                                                                <  2:564576/150‑1 (MQ=255)
          cACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCtt                                                                                                                                            >  2:1025840/1‑150 (MQ=255)
           aCCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTg                                                                                                                                           >  2:765609/1‑150 (MQ=255)
            ccccGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGc                                                                                                                                          >  1:1248113/1‑150 (MQ=255)
             cccGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGcc                                                                                                                                         <  2:223135/150‑1 (MQ=255)
                  gCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCCGGCAAGGATTGccagccc                                                                                                                                    <  1:369407/150‑1 (MQ=255)
                    aaCGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGccggcccgg                                                                                                                                  >  1:70043/1‑150 (MQ=255)
                        cTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGt                                                                                                                              >  2:292190/1‑150 (MQ=255)
                        cTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGt                                                                                                                              >  2:999307/1‑150 (MQ=255)
                                 ggTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCttt                                                                                                                     <  2:1211295/150‑1 (MQ=255)
                                    ggCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCt                                                                                                                  >  1:968285/1‑150 (MQ=255)
                                     gCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTg                                                                                                                 <  1:416664/150‑1 (MQ=255)
                                        gACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCg                                                                                                              >  2:1173263/1‑150 (MQ=255)
                                               tGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATa                                                                                                       >  1:825534/1‑150 (MQ=255)
                                                ggAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCCGGGCCGGCAAGGATTGccagccctgcaatt                                                                                                                             <  1:333954/127‑5 (MQ=38)
                                                ggAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCCGGGCCGGCAAGGATTGccagccctgcaatt                                                                                                                             >  2:333954/1‑123 (MQ=38)
                                                ggAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATaa                                                                                                      >  1:638475/1‑150 (MQ=255)
                                                ggAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATaa                                                                                                      >  1:1357951/1‑150 (MQ=255)
                                                  aaGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATaaga                                                                                                    <  1:404883/150‑1 (MQ=255)
                                                    gagCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATaagaag                                                                                                  >  2:1528147/1‑150 (MQ=255)
                                                      gCTGACCACGACCTGGTATTTTGGAGCAGCGCTGATTGTGGTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGtt                                                                                                <  1:395611/150‑1 (MQ=255)
                                                       cTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTa                                                                                               <  1:947647/150‑1 (MQ=255)
                                                        tGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTAc                                                                                              <  1:636637/150‑1 (MQ=255)
                                                        tGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTAc                                                                                              >  2:79624/1‑150 (MQ=255)
                                                         gACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACa                                                                                             >  1:1052527/1‑150 (MQ=255)
                                                            cACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCa                                                                                          >  2:987741/1‑150 (MQ=255)
                                                               gACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCt                                                                                       <  1:899233/150‑1 (MQ=255)
                                                                aCCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCtt                                                                                      <  1:1012199/150‑1 (MQ=255)
                                                                aCCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCtt                                                                                      <  2:363620/150‑1 (MQ=255)
                                                                   tGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACg                                                                                   <  1:1351034/150‑1 (MQ=255)
                                                                    ggTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGt                                                                                  >  2:666874/1‑150 (MQ=255)
                                                                    ggTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGt                                                                                  >  2:476516/1‑150 (MQ=255)
                                                                    ggTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGt                                                                                  >  2:1315950/1‑150 (MQ=255)
                                                                     gTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTc                                                                                 <  1:765609/150‑1 (MQ=255)
                                                                       aTTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTa                                                                               <  1:55439/150‑1 (MQ=255)
                                                                           tGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACa                                                                           >  2:24308/1‑150 (MQ=255)
                                                                              cgcAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATcg                                                                        <  1:999307/150‑1 (MQ=255)
                                                                                     cTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTAcgc                                                                 <  1:418205/150‑1 (MQ=255)
                                                                                       gATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTAcgccg                                                               >  1:287883/1‑150 (MQ=255)
                                                                                        aTTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTAcgccgc                                                              <  1:1456398/150‑1 (MQ=255)
                                                                                           gtgCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGa                                                           >  2:420957/1‑150 (MQ=255)
                                                                                               tgtgGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCt                                                       >  2:222132/1‑150 (MQ=255)
                                                                                                   gCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCg                                                   <  2:689278/150‑1 (MQ=255)
                                                                                                     tgttgtCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTg                                                 <  2:825534/150‑1 (MQ=255)
                                                                                                       ttgtCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGc                                               <  2:968285/150‑1 (MQ=255)
                                                                                                           cccGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGtt                                           <  1:244323/150‑1 (MQ=255)
                                                                                                           cccGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGtt                                           >  1:261397/1‑150 (MQ=255)
                                                                                                           cccGCGTTCTGGAAAACAACTTTTTGCTCGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACACTGTCGTGGCCGtt                                           >  1:999281/1‑150 (MQ=255)
                                                                                                                   cTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAg                                   >  2:1001799/1‑150 (MQ=255)
                                                                                                                   cTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAg                                   >  2:286234/1‑150 (MQ=255)
                                                                                                                   cTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAg                                   >  2:373229/1‑150 (MQ=255)
                                                                                                                     ggAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTa                                 >  1:3218/1‑150 (MQ=255)
                                                                                                                      gAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAg                                >  2:600034/1‑150 (MQ=255)
                                                                                                                                    gCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTAcgcc                                                                >  2:1447886/1‑104 (MQ=255)
                                                                                                                                    gCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTAcgcc                                                                <  1:1447886/104‑1 (MQ=255)
                                                                                                                                      tGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTc                >  1:21688/1‑150 (MQ=255)
                                                                                                                                      tGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGGCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCGGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTc                <  2:737548/150‑1 (MQ=255)
                                                                                                                                       ggCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTcc               >  1:124679/1‑150 (MQ=255)
                                                                                                                                       ggCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTcc               <  1:1510282/150‑1 (MQ=255)
                                                                                                                                        gCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCca              >  2:101883/1‑150 (MQ=255)
                                                                                                                                         cAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCcac             <  1:1173263/150‑1 (MQ=255)
                                                                                                                                          aaTTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCcaca            <  1:1421134/150‑1 (MQ=255)
                                                                                                                                               cAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAata       <  1:1458129/150‑1 (MQ=255)
                                                                                                                                                 gggCCGGCAAGGATTGCCAGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATc     >  2:1209341/1‑150 (MQ=255)
                                                                                                                                                 gggCCGGCAAGGATTGCCAGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATc     >  2:1214976/1‑150 (MQ=255)
                                                                                                                                                  ggCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATcc    <  2:638475/150‑1 (MQ=255)
                                                                                                                                                   gCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCa   <  1:476516/150‑1 (MQ=255)
                                                                                                                                                    cTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAg  >  2:982492/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTTCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAG  >  NC_000913/824800‑825097

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